Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 250,303 A→G pseudogene (232/756 nt) lafU → flagellar system protein, promoterless fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913250,3030AG96.5% 90.3 / NA 29pseudogene (232/756 nt)lafUflagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (12/16);  total (12/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 = 9.98e-01

CCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCT  >  NC_000913/250174‑250435
                                                                                                                                 |                                                                                                                                    
ccGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCg                                                                                                                            >  1:717233/1‑140 (MQ=255)
       aaCCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATgcg                                                                                                                     >  2:574467/1‑140 (MQ=255)
           tGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCata                                                                                                                 <  1:574467/140‑1 (MQ=255)
                 cGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAg                                                                                                           <  1:445562/140‑1 (MQ=255)
                   tGGCGGGGCCGGAAGAGACGGAGTAAAAAGCGCGTGATGTGAATGAGAAAACGGGGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGATCTGGCAACCAGTTTCATCGCCATTGCCCGTGATGGGCATATGGAAGCg                                                                                                         >  2:276178/1‑140 (MQ=255)
                       ggTGCCGGAAGAGACTGAGAAATAACCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCTCCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATc                                                                                                     <  1:277174/140‑1 (MQ=255)
                       ggTGCCGGAAGAGACGGAGAAAAAACCGCTTGATGTGGATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATc                                                                                                     <  1:361573/140‑1 (MQ=255)
                         tGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCtg                                                                                                   >  2:774759/1‑140 (MQ=255)
                               aaGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCtggagatg                                                                                             <  1:774759/140‑1 (MQ=255)
                                             aaaGCGCGTGATGTGAATGAGACAACGGCGCTGCTCAAGGAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCa                                                                               <  2:209364/140‑1 (MQ=255)
                                             aaaGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCa                                                                               <  2:808891/140‑1 (MQ=255)
                                                gcgcGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAggg                                                                            <  2:13955/140‑1 (MQ=255)
                                                          tgAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTg                                                                  <  1:505280/140‑1 (MQ=255)
                                                          tgAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTg                                                                  <  2:717233/140‑1 (MQ=255)
                                                                  aaaaCGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTaa                                                          >  1:82781/1‑140 (MQ=255)
                                                                         cgcTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTTTTCGCCAGGGATAACGCTTGCTGATTAAAgacgac                                                   >  2:216555/1‑140 (MQ=255)
                                                                         cgcTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAgacgac                                                   <  1:559199/140‑1 (MQ=255)
                                                                         cgcTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAgacgac                                                   <  2:82781/140‑1 (MQ=255)
                                                                                      gaaATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAgacgac                                                   >  1:170677/1‑127 (MQ=255)
                                                                                      gaaATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAgacgac                                                   <  2:170677/127‑1 (MQ=255)
                                                                                      gaaATCCGCCACCGAGCTTGGTCAGGTGGCAACCAGTCTCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAata                                      <  1:276178/140‑1 (MQ=255)
                                                                                        aaTCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATg                                    >  1:682441/1‑140 (MQ=255)
                                                                                        aaTCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATg                                    >  2:51136/1‑140 (MQ=255)
                                                                                         aTCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGt                                   <  2:682441/140‑1 (MQ=255)
                                                                                                   gAGCTTGGTGAGCTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCgg                         <  2:461864/140‑1 (MQ=255)
                                                                                                           tGAGCTGGCAACCAGCATCAACACCCTTGACCGTGATACGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGa                 <  1:216555/140‑1 (MQ=255)
                                                                                                               cTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACGAGAACCGTAATGTGTTTGAACACGGCAGCGCGAAGAt             >  1:259662/1‑140 (MQ=255)
                                                                                                               cTGGCAACCAGTATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGAt             >  1:820200/1‑140 (MQ=255)
                                                                                                                          tATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCt  >  1:545383/1‑140 (MQ=255)
                                                                                                                          tATCAACGCCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCt  >  1:5213/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                    
CCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCT  >  NC_000913/250174‑250435

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