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.3% 260.6 / ‑2.8 59intergenic (‑164/+150)cmtB/tktAPTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (35/23);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTT  >  NC_000913/3077443‑3077584
                                                                      |                                                                       
gtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCata                                                                   >  2:3021946/1‑76 (MQ=255)
  tttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                 >  2:3379600/1‑76 (MQ=255)
    tGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGt                                                               >  1:1776088/1‑76 (MQ=255)
    tGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGt                                                               >  1:2539152/1‑76 (MQ=255)
     gAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTg                                                              <  2:4508576/76‑1 (MQ=255)
       aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAt                                                            >  1:943613/1‑76 (MQ=255)
       aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAACACGCATAGCGTGAt                                                            >  1:754832/1‑76 (MQ=255)
        aaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAta                                                           <  2:3657955/76‑2 (MQ=255)
        aaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATg                                                           <  1:4269122/76‑1 (MQ=255)
         aCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                          <  2:2906960/76‑1 (MQ=255)
          cAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCa                                                         >  1:2698755/1‑76 (MQ=255)
            gTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCata                                                       >  2:4666729/1‑76 (MQ=255)
              gAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATAtt                                                     >  1:1938202/1‑76 (MQ=255)
                attattAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTAc                                                   >  1:3019108/1‑76 (MQ=255)
                 ttattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                  >  2:4450873/1‑76 (MQ=255)
                 ttattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                  >  2:4487756/1‑76 (MQ=255)
                    ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATaa                                               <  1:2493007/76‑1 (MQ=255)
                    ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATaa                                               >  2:506135/1‑76 (MQ=255)
                     taATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAt                                              <  2:3426910/76‑1 (MQ=255)
                      aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATGTTACATAAtt                                             <  1:4444690/76‑1 (MQ=255)
                        tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGa                                           <  2:4100656/76‑1 (MQ=255)
                        tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGa                                           >  1:2543090/1‑76 (MQ=255)
                          tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGata                                         <  2:1661692/76‑1 (MQ=255)
                          tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGata                                         <  1:3021946/76‑1 (MQ=255)
                            ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                       >  1:2249962/1‑76 (MQ=255)
                              atataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattta                                     <  2:2321776/76‑1 (MQ=255)
                               tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttat                                    <  2:3502203/76‑1 (MQ=255)
                                ataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                   >  1:4150075/1‑76 (MQ=255)
                                 taACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTg                                  <  1:769823/76‑1 (MQ=255)
                                   aCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCa                                >  1:3016776/1‑76 (MQ=255)
                                        gAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaa                           >  2:2813726/1‑76 (MQ=255)
                                        gAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaa                           >  1:154918/1‑76 (MQ=255)
                                         aGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                          >  2:2710372/1‑76 (MQ=255)
                                           ttATAGATATAACAACAAGTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                        >  1:4540820/1‑76 (MQ=255)
                                           ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                        >  1:1925849/1‑76 (MQ=255)
                                           ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                        >  2:4488698/1‑76 (MQ=255)
                                           ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                        >  2:4483893/1‑76 (MQ=255)
                                               aGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                    >  1:3862010/1‑76 (MQ=255)
                                                gATATAACAACATATGTCAATCCGCATATCGTGATTCATATTACATAATTGATATTTATTGCATATAAATATTtgt                   <  1:2725295/76‑1 (MQ=255)
                                                  tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                 >  2:4079668/1‑76 (MQ=255)
                                                  tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                 <  1:3664228/76‑1 (MQ=255)
                                                  tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                 >  2:3251211/1‑76 (MQ=255)
                                                   ataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                >  2:1426839/1‑76 (MQ=255)
                                                   ataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                <  1:2813726/76‑1 (MQ=255)
                                                   ataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                >  2:882549/1‑76 (MQ=255)
                                                    taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc               >  2:103797/1‑76 (MQ=255)
                                                     aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATca              >  2:4690265/1‑75 (MQ=255)
                                                          caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaac         <  1:3404206/76‑1 (MQ=255)
                                                           aaaTGTCAATACGCATATCGTGATGCATATTACATCCTTGATATTGATTGCATATAAATATTTGTGCTCTaccaca        >  1:1553003/1‑76 (MQ=255)
                                                           aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaaca        <  2:1314396/76‑1 (MQ=255)
                                                           aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaaca        <  2:1287953/76‑1 (MQ=255)
                                                           aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTGATTGCATATAAATATTTGTGATCTacaaca        >  1:1514820/1‑76 (MQ=255)
                                                            aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac       <  1:3853681/76‑1 (MQ=255)
                                                            aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac       <  1:4337807/76‑1 (MQ=255)
                                                              tGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGc     >  1:2252054/1‑76 (MQ=255)
                                                               gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc    >  1:3840201/1‑76 (MQ=255)
                                                                 cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGTATATAAATATTTGTGATCTACAACACGCCtt  <  2:1072448/76‑1 (MQ=255)
                                                                 cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt  <  2:1047330/76‑1 (MQ=255)
                                                                 cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt  <  2:1042022/76‑1 (MQ=255)
                                                                      |                                                                       
GTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTT  >  NC_000913/3077443‑3077584

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