Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,024,6160GA21.4% 55.7 / 10.4 28intergenic (‑345/+323)PP_3547/emrBNAD(P)‑binding oxidoreductase/multidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base G (12/10);  new base A (2/4);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATC  >  NC_002947/4024493‑4024746
                                                                                                                           |                                                                                                                                   
aaaCAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGt                                                                                                                         >  2:381445/1‑136 (MQ=255)
        gAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAAcgcg                                                                                                                 >  1:317319/1‑136 (MQ=255)
        gAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGAGTCACTTTTTGTCATAcgcg                                                                                                                 >  1:17658/1‑136 (MQ=255)
        gAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGAGTCACTTTTTGTCAAAcgcg                                                                                                                 >  1:384205/1‑136 (MQ=255)
        gAGTTCGCATCTTGCTTCTAGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAAcgcc                                                                                                                 >  2:243846/1‑135 (MQ=255)
         aGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGa                                                                                                                >  2:353246/1‑136 (MQ=255)
            tCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACaa                                                                                                             >  2:45193/1‑136 (MQ=255)
                       ttcttgcttcttgcttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCaa                                                                                                  >  1:217409/1‑136 (MQ=40)
                       ttcttgcttcttgcttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCaa                                                                                                  >  1:186399/1‑136 (MQ=40)
                              ttcttgcttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGc                                                                                           <  1:381445/136‑1 (MQ=40)
                                cttgcttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTg                                                                                         <  1:353246/136‑1 (MQ=40)
                                    cttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTc                                                                                    >  2:167614/1‑136 (MQ=40)
                                    cttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTc                                                                                    >  2:169543/1‑136 (MQ=40)
                                                tgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAATAACGCTGCGCTCCCATCATCCgg                                                                         <  2:17658/136‑1 (MQ=14)
                                                tgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCgg                                                                         <  2:384205/136‑1 (MQ=255)
                                                tgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCgg                                                                         <  2:252366/136‑1 (MQ=255)
                                                         gCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACACTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCATGCGCt                                                               >  2:404546/1‑136 (MQ=255)
                                                         gCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACACTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCATGCGCt                                                               >  2:441356/1‑136 (MQ=255)
                                                                 aTCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACACTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCATGCGCTACGCTCCg                                                       <  1:404546/136‑1 (MQ=255)
                                                                 aTCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACACTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCATGCGCTACGCTCCg                                                       <  1:441356/136‑1 (MQ=255)
                                                                                  tcaCGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTTCCCTCACTCCg                                       <  2:48120/136‑1 (MQ=11)
                                                                                                 aTAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGggag                       <  2:317319/136‑1 (MQ=37)
                                                                                                 aTAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGggag                       <  2:217409/136‑1 (MQ=37)
                                                                                                 aTAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGggag                       <  1:169543/136‑1 (MQ=37)
                                                                                                 aTAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGggag                       <  1:167614/136‑1 (MQ=37)
                                                                                                  tAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGggagg                      <  1:45193/136‑1 (MQ=38)
                                                                                                            tCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGc            <  2:33020/136‑1 (MQ=255)
                                                                                                              cGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCAGCGCGACCCGCCGGGGGTGCCACGCGCCGGGCGTGCGGCCGGGAGGACCCCGCGGCac          >  1:328853/1‑135 (MQ=14)
                                                                                                                      cggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCGACGCTCCGGGGTTCCCTCGCGCCGGGCGTGCTCCCCGGGGGGCCGCGCCGAAGGcggtcta  >  2:296325/1‑130 (MQ=11)
                                                                                                                      cggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCCGCGCTACGCGCCGGGGGTGCCGCGCGCCCGGGGTGGTGCCGGGGGGACCGCGCGGCAGGCCCCCTc  >  1:188404/1‑136 (MQ=14)
                                                                                                                      cggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCCGCGCTACCCGCCGGGGGTGCCGCGCGCCGGGCGTGGTCCGGGGGGGGCCGCGCAGCCGGcccccac  >  1:430288/1‑133 (MQ=14)
                                                                                                                           |                                                                                                                                   
AAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATC  >  NC_002947/4024493‑4024746

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