Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,077,514 Δ1 bp intergenic (‑165/+149) cmtB ← / ← tktA PTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,077,5130C.98.6% 307.5 / ‑4.0 69intergenic (‑164/+150)cmtB/tktAPTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (37/31);  minor base G (0/1);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCT  >  NC_000913/3077438‑3077583
                                                                           |                                                                      
ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAc                                                                        <  2:5259040/76‑1 (MQ=255)
ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAc                                                                        >  1:3536156/1‑76 (MQ=255)
 cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcg                                                                       >  1:1055650/1‑75 (MQ=255)
 cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcg                                                                       >  2:2893578/1‑75 (MQ=255)
 cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcg                                                                       >  2:4213587/1‑75 (MQ=255)
 cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcg                                                                       >  1:2726503/1‑75 (MQ=255)
  ggTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgc                                                                      >  1:5396071/1‑74 (MQ=255)
      ttttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCatat                                                                 <  2:5278035/76‑1 (MQ=255)
      ttttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCatat                                                                 >  2:4654869/1‑76 (MQ=255)
       tttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                <  2:4538065/76‑1 (MQ=255)
       tttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                >  2:535134/1‑76 (MQ=255)
       tttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                >  1:3800552/1‑76 (MQ=255)
        ttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCg                                                               >  1:2880707/1‑76 (MQ=255)
         tGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGt                                                              <  2:490756/76‑1 (MQ=255)
            aaaCAGTGAATTATTAATTTCTTATATAATATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAt                                                           <  2:805495/76‑1 (MQ=255)
            aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAt                                                           >  2:916579/1‑76 (MQ=255)
            aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAt                                                           <  1:1434317/76‑1 (MQ=255)
              aCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                         >  1:3039663/1‑76 (MQ=255)
              aCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                         >  1:3672526/1‑76 (MQ=255)
                aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCat                                                       >  2:4472821/1‑76 (MQ=255)
                aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCat                                                       <  2:46811/76‑1 (MQ=255)
                aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCat                                                       >  1:3090602/1‑76 (MQ=255)
                aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCat                                                       <  1:1739677/76‑1 (MQ=255)
                   gAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATAtt                                                    <  1:4594764/76‑1 (MQ=255)
                     attattAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTAc                                                  >  1:3618415/1‑76 (MQ=255)
                     attattAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTAc                                                  <  1:1605778/76‑1 (MQ=255)
                     attattAATTTCTTATATAACATTGAG‑‑ATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACAt                                                <  1:3389500/76‑1 (MQ=255)
                       tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACAt                                                >  2:3143981/1‑76 (MQ=255)
                       tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACAt                                                >  1:3945995/1‑76 (MQ=255)
                        attaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATa                                               <  2:5705542/76‑1 (MQ=255)
                        attaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATa                                               <  2:3672526/76‑1 (MQ=255)
                           aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAtt                                            <  1:1482967/76‑1 (MQ=255)
                           aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAtt                                            <  1:1446700/76‑1 (MQ=255)
                           aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATGCGCATATCGTGATGCATATTACATAAtt                                            <  2:3689963/76‑1 (MQ=255)
                               tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGata                                        >  2:1078346/1‑76 (MQ=255)
                                 ttATGTAACATTGAGTTATATATATAACAACAACTGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                      <  1:4373416/76‑1 (MQ=255)
                                 ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                      <  1:2334336/76‑1 (MQ=255)
                                 ttATATAACATTGAGTTAGAGATATATCAACAAATGTCAATACGCATATCGTGTTGCATATTGCATAATTGatctg                                      <  2:689329/76‑4 (MQ=255)
                                   atataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattta                                    >  2:93860/1‑76 (MQ=255)
                                    tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttat                                   >  1:148584/1‑76 (MQ=255)
                                    tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttat                                   >  1:452101/1‑76 (MQ=255)
                                      taACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTg                                 >  1:3393060/1‑76 (MQ=255)
                                       aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                >  2:3149115/1‑76 (MQ=255)
                                       aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                <  2:1643512/76‑1 (MQ=255)
                                       aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                <  1:5290367/76‑1 (MQ=255)
                                          aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                             <  2:377119/76‑1 (MQ=255)
                                           ttGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAGGCATATTACATAATTGATATTTATTGCatat                            >  2:2647695/1‑76 (MQ=255)
                                             gAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaa                          >  2:4042366/1‑76 (MQ=255)
                                              aGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                         <  2:3393060/76‑1 (MQ=255)
                                                ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                       <  1:4387108/76‑1 (MQ=255)
                                                 tataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAatat                      >  2:953569/1‑76 (MQ=255)
                                                   taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAttt                    <  1:5298819/76‑1 (MQ=255)
                                                   taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAttt                    >  1:3694154/1‑76 (MQ=255)
                                                   taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAttt                    >  2:568599/1‑76 (MQ=255)
                                                   taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAttt                    <  2:3114416/76‑1 (MQ=255)
                                                    aGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                   >  2:1589057/1‑76 (MQ=255)
                                                    aGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                   <  2:5396070/76‑1 (MQ=255)
                                                      atatAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgtg                 >  2:5185782/1‑76 (MQ=255)
                                                        ataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt               <  2:3846669/76‑1 (MQ=255)
                                                         taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc              <  2:3090602/76‑1 (MQ=255)
                                                          aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCt             >  1:907286/1‑76 (MQ=255)
                                                          aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCt             <  2:3068182/76‑1 (MQ=255)
                                                          aacaacAAATGTCAATACGCATATAGTGATGCATATTACATAATTGATATTTATTGCATATAAATATATGTGATCt             >  1:1270641/1‑76 (MQ=255)
                                                            caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac           >  2:4043888/1‑76 (MQ=255)
                                                             aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTaca          >  1:2515528/1‑76 (MQ=255)
                                                                aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaaca       >  2:257821/1‑76 (MQ=255)
                                                                 aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac      >  1:4379824/1‑76 (MQ=255)
                                                                 aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac      >  2:5712128/1‑76 (MQ=255)
                                                                 aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac      <  1:365771/76‑1 (MQ=255)
                                                                  aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACg     <  2:2187641/76‑1 (MQ=255)
                                                                  aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACg     >  2:5114575/1‑76 (MQ=255)
                                                                  aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACg     <  2:1772899/76‑1 (MQ=255)
                                                                   tGTCAATACGCATATCGTGATGCATATTACATAATTGATATGTATTGCATATAAATATTTGTGATCTACAACACGc    >  1:1152076/1‑76 (MQ=255)
                                                                    gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc   >  2:1473760/1‑76 (MQ=255)
                                                                    gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc   >  1:1005834/1‑76 (MQ=255)
                                                                     tCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCt  >  1:492658/1‑76 (MQ=255)
                                                                           |                                                                      
CCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCT  >  NC_000913/3077438‑3077583

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