Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 250,268 C→T 33.3% pseudogene (197/756 nt) mbhA → ECK0231:JW5812:b0230; flagellar system protein, promoterless fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb250,2680CT33.3% 33.2 / 24.1 30pseudogene (197/756 nt)mbhAECK0231:JW5812:b0230; flagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base C (9/11);  new base T (4/6);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

AGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCA  >  W3110S.gb/250208‑250322
                                                            |                                                      
agaCGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCtt                                              <  1:3410576/71‑1 (MQ=255)
      agaAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAg                                        >  1:1904215/1‑71 (MQ=255)
      agaAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAg                                        >  1:2245690/1‑71 (MQ=255)
      agaAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAg                                        >  1:2241017/1‑71 (MQ=255)
        aaaaaaGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAg                                        >  1:594987/1‑69 (MQ=255)
        aaaaaaGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAg                                        >  1:1373605/1‑69 (MQ=255)
          aaaaGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCTTGGTGAGCTgg                                    >  1:2260377/1‑71 (MQ=255)
          aaaaGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTgg                                    >  1:2301094/1‑71 (MQ=255)
           aaaGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTgg                                    >  1:358724/1‑70 (MQ=255)
              gcgcGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAAc                                <  1:3722486/71‑1 (MQ=255)
               cgcgTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCTTGGTGAGCTGGCAAcc                               <  1:3748834/71‑1 (MQ=255)
               cgcgTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAAcc                               <  1:3344960/71‑1 (MQ=255)
                     aTGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATc                         <  1:576308/71‑1 (MQ=255)
                        tgAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCaa                       >  1:458795/1‑70 (MQ=255)
                            tGAGAAAACGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCt                                               <  1:2436305/42‑1 (MQ=255)
                            tGAGAAAACGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCt                                               <  1:2607089/42‑1 (MQ=255)
                            tGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAAc                                <  1:3859295/57‑1 (MQ=255)
                            tGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAAc                                <  1:1345141/57‑1 (MQ=255)
                              agaAAACGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCTTg                                             >  1:2188926/1‑42 (MQ=255)
                              agaAAACGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCTTg                                             >  1:3483311/1‑42 (MQ=255)
                              agaAAACGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCTTg                                             >  1:1002481/1‑42 (MQ=255)
                              agaAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTTAGCTGGCAACCAGTATCaa                       >  1:3121365/1‑64 (MQ=255)
                                aaaaCGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGc              <  1:1615733/71‑1 (MQ=255)
                                    cGGCGCTGCTGAAGAAGAAATCCGTCACCGAGCTTGGTGAGCTAGCAACCAGTATCAACACCATTGCCCGt          <  1:1626328/71‑1 (MQ=255)
                                     ggCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGt          <  1:1776862/70‑1 (MQ=255)
                                            ctgAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCa  <  1:959739/71‑1 (MQ=255)
                                            ctgAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCa  <  1:3574588/71‑1 (MQ=255)
                                            ctgAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCa  <  1:1313129/71‑1 (MQ=255)
                                            ctgAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCa  <  1:1256758/71‑1 (MQ=255)
                                                 gaagaaATCCGTCACCGAGCTTGGTGAGCTGGCAACCAGTATc                         <  1:2694341/43‑1 (MQ=255)
                                                            |                                                      
AGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCA  >  W3110S.gb/250208‑250322

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