Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,2710GA60.0% 5.7 / 15.5 20H108H (CAC→CAT) ybcVDUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (7/5);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTG  >  NC_000913/579075‑579505
                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
nggAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTAtt                                                                                                                                                                                     >  1:132148/2‑251 (MQ=37)
  gAGTGTCCGGGAATAATTAAAAGCGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGAGTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATTCTCTCAGGAAATAGTAATGAATTATCAAAGGTATAGTAATATCTGTTATGTTCGTGGATAGGTGTAATCCAGCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTAtcg                                                                                                                                                                                    <  1:209056/251‑3 (MQ=38)
                                            cgcgTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGCAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAACTAGTAATGAGTTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGCTTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGt                                                                                                                                          >  1:164546/1‑251 (MQ=25)
                                                                             gTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATATTAATATCGTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGAGTTTCCCGGTATTGCTTAAATGTGATTTTTCTTTATTGCAACCTATCAGAGGACGTTTCTATAAGATGCGTATTTCGTGAGAATTTAAc                                                                                                         <  1:121190/251‑1 (MQ=25)
                                                                              tGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGACAGATTTGTGTAGTGTTATGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATCTTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAAc                                                                                                         >  1:103238/1‑250 (MQ=255)
                                                                                  gACACGGAAGAAACCGGACGTTATGATCTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGACAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTg                                                                                                                                                                >  1:112559/1‑191 (MQ=255)
                                                                                       ggAAGAAACCGGACGTTATGATTTAGNGTGGAANNATNTGTGTAGTGTTCTGAATGCNCTNNNTNNNTNGTAATNAANNATCAAAGGTNTAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCNTTAGCAAGANNTTCCCTGTATTGCTGAAATGTNATTNNTNNTNATTNCa                                                                                                                                                          >  2:47267/1‑192 (MQ=21)
                                                                                       ggAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAgg                                                                                                                                              <  1:47267/204‑1 (MQ=255)
                                                                                                      ttATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGTATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTGCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAgg                                                                                                                                              <  1:138101/189‑1 (MQ=255)
                                                                                                      ttATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTaacc                                                                                                                                                                                                                             <  1:91797/110‑3 (MQ=255)
                                                                                                                    ggAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTACATAGTAATGAATTATCAAAGGGATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGTTTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGTTGAGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGtt                                                                  >  1:10846/1‑251 (MQ=255)
                                                                                                                                          gaaggCTCTCAGTAAATAGTAAATAATTATCATAGGTATAGTAATACCTTTTATGTTCGTTGAGATTTAGAATCCATCGGAAATCTCCTGATTTAGCAAGATTTTCCCCGTATTGCTGAAATGTGATTTCTCTTg                                                                                                                                                                <  1:149050/131‑1 (MQ=255)
                                                                                                                                          gAATGCTGNCNNNANNTAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTNGATGTTTGTATTCCANCGNAAAANTCANGCTTTAGCAAGATTTTCCCTGTATTGCTGACATGTGATTTCTCTTg                                                                                                                                                                >  2:149050/1‑135 (MQ=18)
                                                                                                                                                    aGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCtg                                  >  1:69631/1‑251 (MQ=255)
                                                                                                                                                       aaaTAGTAATGAATTATCAAAGCTATAGTAATATCGTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGc                                                                                                                                                                                   <  1:125019/103‑1 (MQ=255)
                                                                                                                                                          nnnnnntGNNNNNTCANAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAg                                                                                                                                 <  2:132148/144‑1 (MQ=21)
                                                                                                                                                                            ggTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAAtgtg          <  1:45331/251‑1 (MQ=255)
                                                                                                                                                                                  ntAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGAcgt    >  1:131074/2‑249 (MQ=255)
                                                                                                                                                                                   tAATANNTTNNNTCNNCATGGATATTTGTAACCCATCGGNAAACTCCTGCTTTAGCAAGATTTTACCTGNANTNCTg                                                                                                                                                                                 >  2:162816/1‑77 (MQ=17)
                                                                                                                                                                                    aaTATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTtg  <  1:47705/251‑1 (MQ=255)
                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
GGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTG  >  NC_000913/579075‑579505

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑