Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,866,201 G→A 5.5% intergenic (‑210/+20) 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,2010GA 5.5% 300.3 / 5.3 91intergenic (‑210/+20)mtlD/mtlAmannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding/fused mannitol‑specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base G (34/51);  new base A (3/2);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGGTGGTGTG  >  AP009048/3866133‑3866271
                                                                    |                                                                      
caaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                      <  1:1200245/71‑1 (MQ=255)
caaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCACACCTTCTCCAtgtg                                                                      >  1:1312947/1‑71 (MQ=255)
   aTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                      >  1:838799/1‑68 (MQ=255)
   aTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                      >  1:4292489/1‑68 (MQ=255)
    ttGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGgaga                                                                  <  1:3645101/71‑1 (MQ=255)
      gtgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                                 <  1:1332377/70‑1 (MQ=255)
       tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                               >  1:3778970/1‑71 (MQ=255)
       tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                               <  1:2966907/71‑1 (MQ=255)
       tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAgg                                                                <  1:1952358/70‑1 (MQ=255)
       tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAgg                                                                <  1:2653601/70‑1 (MQ=255)
         aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                              <  1:5082437/70‑1 (MQ=255)
         aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                              <  1:2345324/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:4312243/71‑1 (MQ=255)
                      gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                                <  1:5236795/71‑1 (MQ=255)
                      gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                                <  1:3174406/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)
                          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:216850/1‑71 (MQ=255)
                                       cATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                               >  1:4387983/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:2218806/71‑1 (MQ=255)
                                                     ccccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc                 >  1:2937224/1‑71 (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:1681338/71‑1 (MQ=255)
                                                           cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt           <  1:4450701/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:3677309/69‑1 (MQ=255)
                                                                 cATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtggtg       <  1:5333411/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:5066698/1‑71 (MQ=255)
                                                                  aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGgtggtg    <  1:1783070/71‑1 (MQ=255)
                                                                  aTGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGgtggtg    <  1:793763/71‑1 (MQ=255)
                                                                   tgtgGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTc                             <  1:2070336/45‑1 (MQ=255)
                                                                    gtgGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGGTGgtgtg  >  1:2604227/1‑71 (MQ=255)
                                                                    |                                                                      
CAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGGTGGTGTG  >  AP009048/3866133‑3866271

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: