Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0340G.33.3% 33.9 / 22.3 18coding (813/882 nt)ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base G (0/12);  new base . (0/6);  total (0/18)
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.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAAT  >  NC_000913/1373973‑1374101
                                                             |                                                                   
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:998109/68‑1 (MQ=255)
                       tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                        <  1:567478/68‑1 (MQ=255)
                       tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                        <  1:2117345/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)
                                                             |                                                                   
ATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAAT  >  NC_000913/1373973‑1374101

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: ‑