Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0400CT40.5% 13.9 / 42.8 37F273F (TTC→TTT) ‡ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base C (7/15);  new base T (0/15);  total (7/30)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAA  >  NC_000913/1373979‑1374104
                                                             |                                                                
aCAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtgatc                                                               >  1:1365084/1‑65 (MQ=255)
     aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGGTTTTCATCattat                                                          >  1:922526/1‑65 (MQ=255)
     aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCattat                                                          >  1:3643413/1‑65 (MQ=255)
     aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGAACTCGGGGGGGGTCTTCTTAattat                                                          >  1:3945959/1‑65 (MQ=255)
             tgcttttAGTGATTTCGGAAAAGGCGCGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGt                                                 <  1:1677685/65‑1 (MQ=255)
             tccctttAATGATTAGGAAAAAGGCGTGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGt                                                 <  1:778386/63‑1 (MQ=255)
              gcctttAATGGTTGAGAAAAAGCCGGGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTc                                                <  1:607765/65‑1 (MQ=255)
              gcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGt                                                 <  1:868140/65‑1 (MQ=255)
              gcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGt                                                 <  1:73475/65‑1 (MQ=255)
               cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:280525/65‑1 (MQ=255)
               cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:233270/65‑1 (MQ=255)
               cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:531633/65‑1 (MQ=255)
               cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:478792/65‑1 (MQ=255)
               cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:401551/65‑1 (MQ=255)
                   tAATGATTGAGAAAGAGGCGCGTCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTct                                           <  1:1322241/65‑1 (MQ=255)
                    aaTGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                           <  1:54565/65‑1 (MQ=255)
                    aaTGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                           <  1:371406/65‑1 (MQ=255)
                    aaTGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                           <  1:850773/65‑1 (MQ=255)
                     aTGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTGATCATTATTTTCGCTGTCATTctg                                          >  1:2656738/1‑65 (MQ=255)
                      tGATTTAGGAAAAGCGGCGGCGGTCTCGGTGGTGCTCTTCATCATTATTTTCGCCGTCATTctgc                                         >  1:3468165/1‑65 (MQ=255)
                      tGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTctgc                                         <  1:359470/65‑1 (MQ=255)
                         tttAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGa                                      <  1:21285/65‑1 (MQ=255)
                           tAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGAcc                                    <  1:2005200/65‑1 (MQ=255)
                             gggaaaGCGGTGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:314810/62‑1 (MQ=255)
                              gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:2435602/65‑1 (MQ=255)
                                aaaGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGaa                               <  1:2886206/65‑1 (MQ=255)
                                   gcggcgtcgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACgg                           <  1:1047590/65‑1 (MQ=255)
                                        gtcgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa                      <  1:153796/63‑1 (MQ=255)
                                        gtagATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa                      <  1:1056791/62‑1 (MQ=255)
                                          cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:3534090/65‑1 (MQ=255)
                                          cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:825436/65‑1 (MQ=255)
                                          cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:1690612/65‑1 (MQ=255)
                                          agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:519013/64‑1 (MQ=255)
                                          agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:3089781/64‑1 (MQ=255)
                                          agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAAAGGGTGaaaa                    <  1:1440816/64‑3 (MQ=255)
                                          agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGAAAAGGAAAAGGGTGaaaa                    <  1:983827/64‑3 (MQ=255)
                                                  gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATgg             >  1:121684/1‑65 (MQ=255)
                                                    ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca           >  1:304428/1‑65 (MQ=255)
                                                       gCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaaca        >  1:706779/1‑65 (MQ=255)
                                                       gCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaaca        >  1:2096714/1‑65 (MQ=255)
                                                         tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaacaag      >  1:3643754/1‑64 (MQ=255)
                                                         tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa      >  1:3520500/1‑65 (MQ=255)
                                                            tcatcaTTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaa   >  1:2640347/1‑65 (MQ=255)
                                                             catcatTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaaa  <  1:3097116/65‑1 (MQ=255)
                                                             |                                                                
ACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAA  >  NC_000913/1373979‑1374104

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