Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,866,200 T→C 5.5% intergenic (‑209/+21) mtlD ← / ← mtlA mannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding/fused mannitol‑specific PTS enzyme IIABC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,866,2000TC 5.5% 257.1 / 4.9 91intergenic (‑209/+21)mtlD/mtlAmannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding/fused mannitol‑specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base T (35/51);  new base C (3/2);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGGTGGTG  >  AP009048/3866131‑3866269
                                                                     |                                                                     
aacaaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAt                                                                       <  1:3588521/70‑1 (MQ=255)
  caaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                    <  1:1200245/71‑1 (MQ=255)
  caaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCACACCTTCTCCAtgtg                                                                    >  1:1312947/1‑71 (MQ=255)
     aTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                    >  1:4292489/1‑68 (MQ=255)
     aTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                    >  1:838799/1‑68 (MQ=255)
      ttGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGgaga                                                                <  1:3645101/71‑1 (MQ=255)
        gtgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                               <  1:1332377/70‑1 (MQ=255)
         tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                             <  1:2966907/71‑1 (MQ=255)
         tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                             >  1:3778970/1‑71 (MQ=255)
         tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAgg                                                              <  1:2653601/70‑1 (MQ=255)
         tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAgg                                                              <  1:1952358/70‑1 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                            <  1:2345324/70‑1 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                            <  1:5082437/70‑1 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTg                                                           >  1:2532348/1‑71 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTg                                                           >  1:1215358/1‑68 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGgga                                                        <  1:3630135/71‑1 (MQ=255)
                tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgt                                                                     <  1:1160235/56‑1 (MQ=255)
                tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggatt                                                      <  1:2979297/71‑1 (MQ=255)
                  cTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                                     <  1:144184/70‑1 (MQ=255)
                   tGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                                     <  1:4112772/69‑1 (MQ=255)
                        gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTgg                                                                   <  1:3275491/50‑1 (MQ=255)
                        gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                              <  1:3174406/71‑1 (MQ=255)
                        gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                              <  1:4312243/71‑1 (MQ=255)
                        gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                              <  1:5236795/71‑1 (MQ=255)
                         ccTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                              <  1:827125/70‑1 (MQ=255)
                           tGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                               <  1:3803544/51‑1 (MQ=255)
                           tGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttacttac                                           >  1:680712/1‑71 (MQ=255)
                           tGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAtta                                                >  1:2467193/1‑66 (MQ=255)
                            gATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTTCATGTGGAGAGGGTGGGATTGGAttacttac                                           >  1:4207061/1‑70 (MQ=255)
                            gATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAt                                                                       <  1:973503/42‑1 (MQ=255)
                            gATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGgaga                                                                >  1:1454356/1‑49 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACg                                          <  1:3487071/70‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                         >  1:2330691/1‑71 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                         >  1:1068446/1‑71 (MQ=255)
                              tAAGCGTAGCGCATCAGGCAATTAACCACACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGAc                                        >  1:4852193/1‑71 (MQ=255)
                                aGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGgaga                                                                >  1:5347440/1‑45 (MQ=255)
                                aGCGTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGgagaa                                                               <  1:1879473/46‑2 (MQ=21)
                                 gCGTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGGAGAagg                                                             <  1:347914/47‑1 (MQ=18)
                                  cGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATTTGGAGAGGGTGggattgga                                                   >  1:3617350/1‑56 (MQ=255)
                                  cGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGc                                    <  1:4021157/71‑1 (MQ=255)
                                   gTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                              <  1:2238873/60‑1 (MQ=255)
                                   gTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGGAGAAGGTg                                                           >  1:2077999/1‑47 (MQ=17)
                                      gcgcATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                                     >  1:1271706/1‑50 (MQ=255)
                                      gcgcATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcag                                 >  1:4445413/1‑70 (MQ=255)
                                       cgcATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcagc                                <  1:3026443/70‑1 (MQ=255)
                                         cATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                             >  1:4387983/1‑71 (MQ=255)
                                         cATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                             >  1:216850/1‑71 (MQ=255)
                                          aTCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                            <  1:2471884/39‑1 (MQ=255)
                                          aTCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGtt                            >  1:4544558/1‑71 (MQ=255)
                                          aTCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                             <  1:375539/70‑1 (MQ=255)
                                             aGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTcc                          <  1:3100337/70‑1 (MQ=255)
                                             aGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTc                           >  1:2038205/1‑69 (MQ=255)
                                              ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcagca                               >  1:4469543/1‑64 (MQ=255)
                                              ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                        <  1:1964857/71‑1 (MQ=255)
                                               gCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGATCTGCCAGCAGTTCCAg                        >  1:2510022/1‑70 (MQ=255)
                                               gCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                        >  1:4504366/1‑70 (MQ=255)
                                                cAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                      <  1:2791445/71‑1 (MQ=255)
                                                 aaTTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcagc                                <  1:101380/60‑1 (MQ=255)
                                                    tAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttactta                                            >  1:1445527/1‑45 (MQ=255)
                                                    tAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACtt                   <  1:4414191/70‑1 (MQ=255)
                                                     aaCCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCa                 >  1:141234/1‑71 (MQ=255)
                                                      aCCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCa                 <  1:591826/70‑1 (MQ=255)
                                                       tcccACCCTCTCCACATGGAGAAGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc               >  1:1312819/2‑71 (MQ=25)
                                                       ccccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc               >  1:2937224/1‑71 (MQ=255)
                                                       ccccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc               <  1:2218806/71‑1 (MQ=255)
                                                        cccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc               >  1:3572972/1‑70 (MQ=255)
                                                         ccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                             <  1:4366777/55‑1 (MQ=255)
                                                         ccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTc                  <  1:1991114/66‑1 (MQ=255)
                                                          cACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCAc            <  1:3502840/71‑1 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCAc            >  1:3682374/1‑70 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACg           >  1:1283496/1‑71 (MQ=255)
                                                            ccTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcag                                 <  1:2719237/48‑1 (MQ=255)
                                                            ccTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTcc                          >  1:1005395/1‑55 (MQ=255)
                                                            ccTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGc          <  1:1727949/71‑1 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGc                       <  1:2420796/57‑1 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACt                    <  1:1174665/60‑1 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt         >  1:1874980/1‑71 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt         <  1:4943847/71‑1 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt         <  1:4450701/71‑1 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt         <  1:1681338/71‑1 (MQ=255)
                                                              ttCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGAc                                        <  1:1412/39‑1 (MQ=255)
                                                              ttCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtg        >  1:4335628/1‑71 (MQ=255)
                                                                ctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCt                                      >  1:5070236/1‑39 (MQ=255)
                                                                ctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCa                         >  1:642705/1‑52 (MQ=255)
                                                                ctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACt                    >  1:4896151/1‑57 (MQ=255)
                                                                ctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtgg       <  1:104236/70‑1 (MQ=255)
                                                                ctcCACATGGAGAAGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtggt      >  1:5330517/1‑71 (MQ=39)
                                                                  ccATGTGGAGAGGGTGGGATTGGATTACTTACGAcc                                       <  1:4380401/36‑1 (MQ=255)
                                                                   cATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCATCACTTc                  <  1:691583/56‑1 (MQ=255)
                                                                   cATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtggtg     <  1:5333411/69‑1 (MQ=255)
                                                                   cATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtggtg     <  1:3677309/69‑1 (MQ=255)
                                                                   cATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtggtg     <  1:225612/69‑1 (MQ=255)
                                                                    aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcagca                               <  1:1019998/42‑1 (MQ=255)
                                                                    aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtggtgg    >  1:1285971/1‑69 (MQ=255)
                                                                    aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGgtggtg  <  1:1831368/71‑1 (MQ=255)
                                                                    aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGgtggtg  <  1:793763/71‑1 (MQ=255)
                                                                    aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGgtggtg  <  1:1783070/71‑1 (MQ=255)
                                                                    aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGgtggtg  >  1:5066698/1‑71 (MQ=255)
                                                                     tgtgGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTc                           <  1:2070336/45‑1 (MQ=255)
                                                                     |                                                                     
AACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGGTGGTG  >  AP009048/3866131‑3866269

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: