Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,875,8380AC40.0% 23.7 / 12.5 20T56P (ACC→CCC) nimT2‑nitroimidazole exporter
Reads supporting (aligned to +/- strand):  ref base A (11/1);  new base C (0/8);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGTTACTGGATACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGTAGCGCGACGTTTTGG  >  NC_000913/1875774‑1875902
                                                                |                                                                
aCCGTTACTGGATACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGacca                                                                  >  1:289663/1‑65 (MQ=255)
aCCGTTACTGGATACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGacca                                                                  >  1:386865/1‑65 (MQ=255)
aCCGTTACTGGATACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGacca                                                                  >  1:1618641/1‑65 (MQ=255)
      aCTGAATACGTTTCGTTCCGCTTCCTCGCTGGCGACAGCCCAACCCGGCTTATTGACCCCCCTGc                                                            <  1:3857585/65‑1 (MQ=255)
       cTGGATACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGACCACCCTggc                                                           >  1:3536662/1‑63 (MQ=255)
        tGGATACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCCGCTTATTGACCAACCTGCCa                                                          >  1:2613710/1‑65 (MQ=255)
        tGGATAAGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAAACGGCTTATTGACCACCCTGCCa                                                          >  1:2628117/1‑65 (MQ=255)
          taTACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAACCCGGCTTATTGCCCCCCCTGCCatt                                                        <  1:863604/64‑1 (MQ=255)
                aTTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGc                                                  >  1:2954951/1‑65 (MQ=255)
                     ttCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCttcg                                             >  1:1099098/1‑63 (MQ=255)
                      tCCGCTTACTCGCTGACGACAGCGCAACCCGGCTTATTCCCCCCCCTGCCATTATTGGCCTTTgc                                            <  1:2814081/65‑1 (MQ=255)
                      tCCGCTTACTCGCTGACGACAGCGAAAACCGGCTTATTGACCCCCCTGCCATTATTGGCCTTTgc                                            <  1:3274729/65‑1 (MQ=255)
                      tCCGCTTACCCGCTGCCGACAGCGCAAACCGGCTTATTGCCCCCCCTGCCATTATTGGCCTTTgc                                            <  1:1439679/65‑1 (MQ=255)
                      tCCGCTGACTCGCTGACGACAGCGCAAACCGGCTTATTGACCCCCCTGCCATTATTGGCCTTTgc                                            <  1:1640503/65‑1 (MQ=255)
                      cccGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGACCCCCCTGCCATTATTGGCCTTTgc                                            <  1:1626278/64‑1 (MQ=255)
                           ttACTCGCTGGCGACGCCGCAAACCGGCTTATTGACCGCCCTGCCATTATTGGCCTTTGCGCTaa                                       <  1:2046432/65‑1 (MQ=255)
                           ttACTCGCTGACGACAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTaa                                       <  1:4016141/65‑1 (MQ=255)
                           ttACTCGCTGACGACAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTaa                                       <  1:1837027/65‑1 (MQ=255)
                             ccTCGCTGACGACAGCGCAACCCGTCTTATTGCCCCCCCTGCCATTATTGGCCTTTGCGCTAAtc                                     <  1:2192326/64‑1 (MQ=255)
                                  cgcacgacAGCGCAAACCGGCTTATTCCCCCCCCTGCCATTATTGGCCTTTGCGCTAATCTCAcc                                <  1:1584362/62‑1 (MQ=255)
                                  cTGACGACAGCGCAACCCGGCTTATTGCCCCCCCTGCCATTATTGGCCTTTGCGCTAATCTCAcc                                <  1:3299899/65‑1 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGGCTGTGGGCTAATTTCACCTTTGGc                          >  1:1639859/1‑65 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTgcc                          >  1:3600128/1‑63 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGc                          >  1:320523/1‑65 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGc                          >  1:2159792/1‑65 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGc                          >  1:3632235/1‑65 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGc                          >  1:2039687/1‑65 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGc                          >  1:643355/1‑65 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGc                          >  1:1370890/1‑65 (MQ=255)
                                        acAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGc                          >  1:1152639/1‑65 (MQ=255)
                                          ggcgcAACCCGTCTTATTGCCCCCCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTg                        <  1:2544161/64‑1 (MQ=255)
                                           gcgcAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGc                       >  1:3532294/1‑65 (MQ=255)
                                            cccAACCCGGTTTATTGCCCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGcc                      <  1:2473238/63‑1 (MQ=255)
                                                  ccGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGTa                >  1:2520148/1‑65 (MQ=255)
                                                         tttGCCCCCCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGTAGCGCGAc         <  1:883047/64‑1 (MQ=255)
                                                         aTTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCCGTAGcgcgcc         >  1:2172572/1‑63 (MQ=255)
                                                           tGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGTAGCGCGACGt       >  1:1962812/1‑65 (MQ=255)
                                                                accCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGTAGCGCGACGTTTTgg  >  1:2149701/1‑65 (MQ=255)
                                                                accCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGTAGCGCGACGTTTTgg  >  1:1321682/1‑65 (MQ=255)
                                                                accCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGGAGCGCGACGTTTTgg  >  1:3029793/1‑65 (MQ=255)
                                                                |                                                                
ACCGTTACTGGATACGATTCGTTCCGCTTACTCGCTGACGACAGCGCAAACCGGCTTATTGACCACCCTGCCATTATTGGCCTTTGCGCTAATCTCACCTTTGGCTGCCCCGGTAGCGCGACGTTTTGG  >  NC_000913/1875774‑1875902

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