Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,136,1140TA31.2% 12.1 / 10.1 16V184D (GTT→GAT) flgHflagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base A (2/3);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAGCGCAAAATATGGG  >  NC_000913/1135978‑1136222
                                                                                                                                        |                                                                                                            
aaTGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGATTAACCACGcca                                                                                                  <  4:19523/149‑3 (MQ=255)
              cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGATTAACCACGccaga                                                                                                <  1:108328/137‑5 (MQ=255)
              cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGATTAACCACGccaga                                                                                                >  2:108328/1‑133 (MQ=255)
                          gTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATaccgtaccgt                                                                        >  3:22605/1‑149 (MQ=255)
                              acggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTAc                                                                    <  3:164862/149‑1 (MQ=255)
                                                      gTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCgg                                                         >  1:79856/1‑136 (MQ=255)
                                                      gTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCgg                                                         <  2:79856/136‑1 (MQ=255)
                                                      gTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCAtt                                              <  4:22605/147‑1 (MQ=255)
                                                           cGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCaa                                                                                   <  3:121740/105‑1 (MQ=255)
                                                           cGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCaa                                                                                   >  4:121740/1‑105 (MQ=255)
                                                                  cTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCaa                                <  7:6679/149‑1 (MQ=255)
                                                                                 gAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTc                                                                       >  7:123045/1‑95 (MQ=255)
                                                                                 gAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTc                                                                       <  8:123045/95‑1 (MQ=255)
                                                                                                aTTAATCAGGGTACCGAATTTATTCGCATCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCAACATTAACGAAGCGCAAAATATggc  >  1:17338/1‑148 (MQ=255)
                                                                                                                                 tgCGTGGATTAACCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAGCGCAAAat       <  1:73469/110‑1 (MQ=255)
                                                                                                                                 tgCGTGGATTAACCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAGCGCAAAat       >  2:73469/2‑111 (MQ=255)
                                                                                                                                        |                                                                                                            
AATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAGCGCAAAATATGGG  >  NC_000913/1135978‑1136222

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