Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,8310CG21.1% 84.4 / 19.6 38intergenic (+30/+109)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (17/13);  new base G (4/4);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCG  >  NC_000913/2282685‑2282976
                                                                                                                                                  |                                                                                                                                                 
tACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg                                                                                                                                                 >  5:77289/1‑149 (MQ=255)
      gAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCgggg                                                                                                                                           >  2:48958/1‑149 (MQ=255)
          ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGtttt                                                                                                                                       <  6:77289/149‑1 (MQ=255)
          ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGtttt                                                                                                                                       >  6:153956/1‑149 (MQ=255)
               aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACCAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGAGCGGGGttttt                                                                                                                                      >  1:115277/1‑145 (MQ=255)
               aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGtttttt                                                                                                                                     <  5:153956/146‑1 (MQ=255)
               aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGAGCGGGGttttt                                                                                                                                      >  5:125904/1‑145 (MQ=255)
               aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGAGCGGGGttttt                                                                                                                                      <  2:115277/145‑1 (MQ=255)
               aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGAGCGGGGttttt                                                                                                                                      <  6:125904/145‑1 (MQ=255)
                         ttgACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGAtt                                                                                                                        <  1:48958/149‑1 (MQ=255)
                                                aaaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGc                                                                                                                             <  4:151695/121‑1 (MQ=255)
                                                aaaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGc                                                                                                                             >  3:151695/1‑121 (MQ=255)
                                                aaaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                                >  2:55257/1‑134 (MQ=255)
                                                aaaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                                <  1:55257/134‑1 (MQ=255)
                                                    cTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGAGCGGGGtttttta                                                                                                                                    <  3:117143/110‑2 (MQ=39)
                                                    cTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGAGCGGGGtttttta                                                                                                                                    >  4:117143/1‑109 (MQ=39)
                                                     tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTg                                                                                                                              >  2:193868/1‑115 (MQ=255)
                                                     tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTg                                                                                                                              <  1:193868/115‑1 (MQ=255)
                                                                       cccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                                >  7:82213/1‑111 (MQ=255)
                                                                       cccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                                <  8:82213/111‑1 (MQ=255)
                                                                        ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGttttt                                                                                                                                      >  6:21470/1‑88 (MQ=255)
                                                                                           agcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                   >  8:71559/1‑136 (MQ=255)
                                                                                           agcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                   <  7:71559/136‑1 (MQ=255)
                                                                                            gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgc                                                                    >  5:89014/1‑134 (MQ=255)
                                                                                            gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgc                                                                    <  6:89014/134‑1 (MQ=255)
                                                                                                aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTg                                                                                           >  2:96430/1‑107 (MQ=255)
                                                                                                aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTg                                                                                           <  1:96430/107‑1 (MQ=255)
                                                                                                   gACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAacac                                               >  2:130095/1‑148 (MQ=255)
                                                                                                            cgcgCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAt                                      <  1:130095/148‑1 (MQ=255)
                                                                                                              cgcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                   >  2:144353/1‑149 (MQ=255)
                                                                                                                cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                                 <  1:144353/149‑1 (MQ=255)
                                                                                                                     ttGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgccgc                            >  6:98988/1‑149 (MQ=255)
                                                                                                                           aaTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACg                                                                                                           >  3:136906/1‑64 (MQ=255)
                                                                                                                           aaTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACg                                                                                                           <  4:136906/64‑1 (MQ=255)
                                                                                                                                     aaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAt                                      <  7:36751/123‑1 (MQ=255)
                                                                                                                                     aaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAt                                      >  8:36751/1‑123 (MQ=255)
                                                                                                                                       aaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGCCGTCCGGCGCTGc           >  8:43439/1‑148 (MQ=255)
                                                                                                                                               gCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTcg  >  5:185698/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCG  >  NC_000913/2282685‑2282976

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