Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 248,447 C→T pseudogene (1624/1713 nt) lfhA ← flagellar system protein, promoterless fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913248,4470CT96.5% 87.5 / ‑6.3 29pseudogene (1624/1713 nt)lfhAflagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/8);  total (21/8)
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

TCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATAT  >  NC_000913/248319‑248563
                                                                                                                                |                                                                                                                    
tCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATaa                                                                                                                  >  2:22894/1‑133 (MQ=255)
 cACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAACCACGGCGCACACAAGATTGtatacc                                                                                                                 >  2:313677/1‑131 (MQ=255)
    gccgcATCCGGCAGTTGTGGACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAAcgcg                                                                                                              <  1:516057/133‑1 (MQ=255)
      cacATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGcc                                                                                                            >  1:80805/3‑133 (MQ=255)
         aTCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAACCCCGGCGCAAACAACCTTGTATAACCCGCCAGc                                                                                                         >  2:640468/1‑133 (MQ=255)
          tacGGCAGTTGTGCAACAATGCTGCATACTTATATCAACGCTCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCGTGTATAACGCGCCAGCa                                                                                                        <  2:291664/131‑1 (MQ=255)
                            aTGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAACCAACCTTGCATAACGCGCCAGCAAATGGCGCAGCTGCGGTg                                                                                      >  2:216505/1‑133 (MQ=255)
                              gCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGc                                                                                    <  1:36973/133‑1 (MQ=255)
                              gCTGGATACTTACATCAACGCCCCCATAATATTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGc                                                                                    <  1:313677/133‑1 (MQ=255)
                                    tACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCTAGCAAAGAGCGTAGATGCGGTAGCACcagc                                                                              >  1:299128/1‑133 (MQ=255)
                                    tACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACcagc                                                                              <  2:405755/133‑1 (MQ=255)
                                      cTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCTAGCAACGGGCGCAGCTGCGGTAGCACcagcag                                                                            >  2:302295/1‑133 (MQ=255)
                                         aCATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCAc                                                                         <  2:453940/133‑1 (MQ=255)
                                              aaCGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACcag                                                                               <  1:379003/122‑1 (MQ=255)
                                              aaCGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACcag                                                                               >  2:379003/1‑122 (MQ=255)
                                                cGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGCTGGTACTAGCAGCTCCGGGTCt                                                                  >  1:296859/1‑133 (MQ=255)
                                                          aTCTTCAACTCCAGCCCACCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATACCGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCccgccg                                                        <  1:250717/133‑1 (MQ=255)
                                                                     cAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCGGCAGCGGGTCTTTCCCAGCCGCTTTCATCTGc                                             >  2:250717/1‑133 (MQ=255)
                                                                       gCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCAGCGGGTCTATCCCCGTCGCTTTAATCTGctc                                           >  1:211959/1‑133 (MQ=255)
                                                                            tCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTATCATCTGCTCTTTCa                                      >  1:313585/1‑133 (MQ=255)
                                                                                gCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATACCGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCCCCGGGTTTTCCCCCGCCGCTTTCATCTGCTCTTTCACCtg                                  <  2:324964/133‑1 (MQ=255)
                                                                                     tCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGAGGCa                             >  2:325287/1‑133 (MQ=255)
                                                                                        ttATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCg                          <  2:80805/133‑1 (MQ=255)
                                                                                         tataCGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCCCCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGt                         <  1:22894/133‑1 (MQ=255)
                                                                                              gACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTActct                    >  1:325744/1‑133 (MQ=255)
                                                                                                cAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTgg                  >  2:406268/1‑133 (MQ=255)
                                                                                                cAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTgg                  >  2:450004/1‑133 (MQ=255)
                                                                                                  gCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGaa                >  2:518299/1‑133 (MQ=255)
                                                                                                  gCAAATCCAGCCGCGGCCCAAACACCCTTTTAAAGCGCGCCAAAAAAGGGCGCAGCCGGGGCGCCACCGCCACCCCCGGTTCTTCCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGaa                <  2:477276/133‑1 (MQ=255)
                                                                                                   cacaTGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAAc               >  1:122449/1‑133 (MQ=255)
                                                                                                   cacaTGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAAc               >  1:56054/1‑133 (MQ=255)
                                                                                                       tGCAGCCCCGTCGCACACACCCTTCTATCACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCCCCGGGTCTTCCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGt           <  2:181029/133‑1 (MQ=255)
                                                                                                       tGCACCCCGGGCGCAAACACCCTTGTATACCCCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCACCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGt           <  2:325744/133‑1 (MQ=255)
                                                                                                        gCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGtt          >  2:451281/1‑133 (MQ=255)
                                                                                                        gCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGtt          >  2:590413/1‑133 (MQ=255)
                                                                                                          aGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGAATCTTATTCTGGAACTGGTTGa        >  1:213473/1‑133 (MQ=255)
                                                                                                                ggCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCatat  >  1:396169/1‑133 (MQ=255)
                                                                                                                                |                                                                                                                    
TCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATAT  >  NC_000913/248319‑248563

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