Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0380TC28.6% 26.5 / 10.5 21F273L (TTC→CTC) ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base T (3/12);  new base C (0/6);  total (3/18)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAA  >  NC_000913/1373973‑1374104
                                                                 |                                                                  
atCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtc                                                                  <  1:260274/68‑1 (MQ=255)
 tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGTGGGGGTCTtca                                                                 >  1:4141562/1‑68 (MQ=255)
 tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGGGGGGCTCTtct                                                                 >  1:2535576/1‑67 (MQ=255)
 tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGGTCTtcc                                                                 >  1:3904685/1‑67 (MQ=255)
 tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTcca                                                                 >  1:3608202/1‑68 (MQ=255)
 tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCTCttta                                                                 >  1:2595337/1‑66 (MQ=255)
 tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCCCTtca                                                                 >  1:1438333/1‑68 (MQ=255)
 tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGGGCTCTtca                                                                 >  1:3331519/1‑68 (MQ=255)
      aCAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCTCTTCATCatt                                                            >  1:3573142/1‑68 (MQ=255)
               cctttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCt                                                   <  1:387412/68‑1 (MQ=255)
                ctttgcctttAATGATTTAGGAAAAGCGGCGGCGAACTCGGGGGTGCTCTTCATCATTATTTTCCCTg                                                  >  1:3192424/1‑68 (MQ=255)
                  ttgcgtttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:81622/68‑1 (MQ=255)
                  ttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:3744485/68‑1 (MQ=255)
                       tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                           <  1:2117345/68‑1 (MQ=255)
                       tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                           <  1:567478/68‑1 (MQ=255)
                       tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                           <  1:998109/68‑1 (MQ=255)
                              aTTTAAGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGAcc                                    <  1:4174022/68‑1 (MQ=255)
                                ttAGGAAAAGCGGCGGCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:3473481/68‑1 (MQ=255)
                                ttAGAAAAAGGGGCGTCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:2656267/68‑1 (MQ=255)
                                 tAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:297562/68‑1 (MQ=255)
                                 tAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:2640781/68‑1 (MQ=255)
                                  aGGCAAGGCGGCGGCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGaa                               <  1:3380252/68‑1 (MQ=255)
                                    gAAAAGCGGCGGCGGTCTCGGGGGTGGTTTTTATCATTATTTTTGCTGTTACTCTTCTTACCACGaaa                              >  1:1681246/1‑68 (MQ=16)
                                      aaaGCGGCGGCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACgg                           <  1:3168278/68‑1 (MQ=255)
                                       aaGCGGCGGCGATCTCGGTGGTGCTCTTCTTCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACgg                           <  1:1858166/68‑1 (MQ=255)
                                        aGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACggg                          <  1:512505/68‑1 (MQ=255)
                                           ggcggcgATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa                      <  1:1999761/68‑1 (MQ=255)
                                           cgcgtcgATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa                      <  1:4181195/67‑1 (MQ=255)
                                            gcgccgATCTGGGGGGTGCTCTTCATCATTATTTTCGC‑GTCATTCTGCTGACCAGGAAACGGGTGAAc                     <  1:2073098/68‑1 (MQ=255)
                                             tgtcgATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:2980227/65‑1 (MQ=255)
                                               gcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCt                   <  1:2732091/68‑1 (MQ=255)
                                                      cGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGc            <  1:1678370/68‑1 (MQ=255)
                                                        gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaa          >  1:2525973/1‑68 (MQ=255)
                                                          ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTAAATGCaaaca        >  1:3534070/1‑68 (MQ=255)
                                                             gCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCACCAAAt     >  1:3292700/1‑68 (MQ=255)
                                                               tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaa   >  1:3677682/1‑68 (MQ=255)
                                                               tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaa   >  1:3446694/1‑68 (MQ=255)
                                                               tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaa   >  1:1522468/1‑68 (MQ=255)
                                                                ctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaaa  <  1:4046534/68‑1 (MQ=255)
                                                                ctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaaa  <  1:3230895/68‑1 (MQ=255)
                                                                 |                                                                  
ATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAA  >  NC_000913/1373973‑1374104

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