Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,022,5960AC36.1% 31.1 / 27.2 36intergenic (‑291/+377)PP_3547/emrBNAD(P)‑binding oxidoreductase/multidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base A (12/11);  new base C (6/7);  total (18/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 = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGC  >  NC_002947/4022449‑4022739
                                                                                                                                                   |                                                                                                                                               
tttGACAAAAAGTGACTCGCCGTATGGGCGAAAAGGTGAGTCCGCGTCGCCATCGCAAACGGACGGTTCGCGGCGTCAAAACCGACGAGTTCGCATCTTGCGTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGa                                                                                                                                                 >  5:70928/1‑148 (MQ=18)
tttGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCGAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGa                                                                                                                                                 <  8:2646/148‑1 (MQ=255)
tttGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGa                                                                                                                                                 <  7:125525/148‑1 (MQ=255)
   gACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTa                                                                                                                                              >  3:14399/1‑148 (MQ=255)
    aCAAAAAGTGACTCGCCGTAAGGGCGAAATGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGCCGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGGTTCTTGGATGTGGGCATTCACATCGAGTAg                                                                                                                                             >  3:67510/1‑148 (MQ=255)
               cTCGCCGTAAGGGCGAAACGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTattcattc                                                                                                                                   <  3:111028/147‑1 (MQ=255)
                tCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCAc                                                                                                                                 >  7:23681/1‑148 (MQ=255)
                  gCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCAc                                                                                                                                 >  8:31682/1‑146 (MQ=255)
                    cGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAg                                                                                                                             >  2:52478/1‑148 (MQ=255)
                       acgggCGAAACGGGGAGGCCCCGTCGCCATCGCAACCGGACGGTGCGCGCAGTCAACCCAGACGAGTGCGCATCGTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGCTGGGGGCATTCACATCGAGTAGGTATTCATTCACGCAgggg                                                                                                                          <  8:157931/146‑3 (MQ=17)
                         gggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTgg                                                                                                                         <  3:155174/147‑1 (MQ=255)
                              aaaaGGTGAGTCATCGTCGCCATCGGATACGGACGGTTGGGTGAGTCAAATCAGACGAGTTCTCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCAGATCGAGTAGGTATTCATTCACGGAGGGGGCGCATa                                                                                                                   <  8:142913/148‑1 (MQ=17)
                              aaaaGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATa                                                                                                                   <  3:146487/148‑1 (MQ=255)
                                   gTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCAc                                                                                                              >  7:40689/1‑148 (MQ=255)
                                   gTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCAc                                                                                                              <  4:128973/148‑1 (MQ=255)
                                   gTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCAc                                                                                                              <  1:98252/148‑1 (MQ=255)
                                     gAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCAcct                                                                                                            <  6:53988/148‑2 (MQ=255)
                                         cAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGTATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTcc                                                                                                        <  1:145008/148‑1 (MQ=255)
                                                     cGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTAcggcg                                                                                            >  4:96590/1‑148 (MQ=255)
                                                       cAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTAcggcgg                                                                                           <  2:87751/147‑1 (MQ=255)
                                                       cAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTAcggcgg                                                                                           <  3:45995/147‑1 (MQ=255)
                                                       cAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGACACCCTTCCGCCCTTACGGCggg                                                                                          >  1:152408/1‑148 (MQ=255)
                                                       cAAACGGACGGTTCGCTGAGGCAAAACAGACGAGTGCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGGGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTAcggcgg                                                                                           <  7:80714/147‑1 (MQ=255)
                                                             gacggTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACtt                                                                                    >  3:88521/1‑148 (MQ=255)
                                                                ggTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACttttt                                                                                 <  6:176079/148‑1 (MQ=255)
                                                                ggTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACttttt                                                                                 <  8:146246/148‑1 (MQ=255)
                                                                     gCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCaa                                                                            >  6:132009/1‑148 (MQ=255)
                                                                          gTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAAcgc                                                                        <  7:7252/147‑1 (MQ=255)
                                                                           tCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGGATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCAACTTTACGCCCTTAAGGCGTGTCAATTTTTGTCAAAcgcg                                                                       >  7:147215/1‑147 (MQ=14)
                                                                           tCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAAcgcg                                                                       >  5:130963/1‑147 (MQ=255)
                                                                           tCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGa                                                                      <  4:88521/148‑1 (MQ=255)
                                                                               aaCAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACaa                                                                   >  3:76871/1‑147 (MQ=255)
                                                                                 cAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACaaaaa                                                                <  7:134029/148‑1 (MQ=255)
                                                                                                                   ttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGcc                              >  3:148146/1‑148 (MQ=39)
                                                                                                                      ttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCt                           >  7:32231/1‑148 (MQ=39)
                                                                                                                             ttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTAc                    >  5:9079/1‑148 (MQ=39)
                                                                                                                                     gggCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCggg            <  6:70928/148‑1 (MQ=39)
                                                                                                                                        cATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGtt         >  7:83806/1‑148 (MQ=39)
                                                                                                                                               aTCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTACCTCGc  >  5:106246/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGC  >  NC_002947/4022449‑4022739

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