Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,283,2300AG23.3% 51.3 / 13.6 30intergenic (‑155/+23)yjcF/actPpentapeptide repeats protein/acetate transporter
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base G (3/4);  total (13/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 = 9.22e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACC  >  NC_000913/4283087‑4283368
                                                                                                                                               |                                                                                                                                          
ataATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAgacga                                                                                                                                       >  2:105817/1‑149 (MQ=255)
 taATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAgacgac                                                                                                                                      <  6:175135/149‑1 (MQ=255)
    ttAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCgg                                                                                                                                   <  1:105817/149‑1 (MQ=255)
         aagaaagaAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTATTATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGaa                                                                                                                               >  3:4581/1‑148 (MQ=255)
         aagaaagaAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAg                                                                                                                              >  4:143343/1‑149 (MQ=255)
          agaaagaaAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTATTATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGaata                                                                                                                             <  4:4581/149‑3 (MQ=255)
           gaaagaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGa                                                                                                                                >  3:10271/1‑145 (MQ=255)
           gaaagaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGa                                                                                                                                <  4:10271/145‑1 (MQ=255)
           gaaagaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAggg                                                                                                                            <  6:169655/149‑1 (MQ=255)
               gaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGaatactctt                                                                                                                        <  7:17692/149‑8 (MQ=255)
                     gCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTaa                                                                                                                  >  4:132572/1‑149 (MQ=255)
                      cATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAAt                                                                                                                 <  6:6941/149‑1 (MQ=255)
                      cATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAAt                                                                                                                 <  3:143343/149‑1 (MQ=255)
                       atatTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAggg                                                                                                                            <  8:43853/137‑1 (MQ=255)
                       atatTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAggg                                                                                                                            >  7:43853/1‑137 (MQ=255)
                                                   cTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCg                                                                                    >  3:60866/1‑149 (MQ=255)
                                                                                       gTTTAATATTATGGGAAAAGGTGCTGCATTTGGGAGAGGAAGAGTATGCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTtc                                                 <  7:86739/148‑1 (MQ=255)
                                                                                           aaTATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGGCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTtcacgc                                             <  3:132572/148‑1 (MQ=255)
                                                                                             tattatGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTtcacgctca                                          >  6:65804/1‑149 (MQ=255)
                                                                                                    ggAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACgcgcgc                                   <  5:65804/149‑1 (MQ=255)
                                                                                                                           agGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            <  4:60866/149‑1 (MQ=255)
                                                                                                                             gAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGc                                                                          >  2:148201/1‑85 (MQ=255)
                                                                                                                             gAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGc                                                                          <  1:148201/85‑1 (MQ=255)
                                                                                                                             gAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            >  4:113224/1‑147 (MQ=255)
                                                                                                                              aaGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            >  6:15582/1‑146 (MQ=255)
                                                                                                                                     ttCCCCGGTCGTCTGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGcc                                  >  2:26927/1‑117 (MQ=255)
                                                                                                                                     ttCCCCGGTCGTCTGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGcc                                  <  1:26927/117‑1 (MQ=255)
                                                                                                                                     ttCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAAcc  <  5:15582/149‑1 (MQ=255)
                                                                                                                                      tCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            <  2:150970/138‑1 (MQ=255)
                                                                                                                                      tCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            >  1:150970/1‑138 (MQ=255)
                                                                                                                                               |                                                                                                                                          
ATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACC  >  NC_000913/4283087‑4283368

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