Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,1190AG46.2% ‑6.7 / 13.0 14intergenic (+33/+132)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (3/3);  total (8/6)
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.12e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCC  >  NC_000913/579874‑580327
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                
gACAAAATGAATAAAGAACAATCTGCTGATGACCCCTCCGAGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATCTTGTAACTGCATGTATAGAAAATAAGGAGTCTCTGGAAGCATTCCGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTACCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCtgt                                                                                                                                                                                                             >  1:119639/1‑251 (MQ=255)
         aaTAAAGAACAATCTGCTGTTGATCCCGCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCAAGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAATCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGGCTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGcc                                                                                                                                                                                                    <  1:174109/251‑1 (MQ=255)
                                                      gtAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGan                                                                                                                                                                                                          <  1:103855/200‑2 (MQ=255)
                                                       tAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTc                                                                                                                                                       >  1:69141/1‑250 (MQ=255)
                                                                 gCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtc                                                                                                                                                                                     >  1:74622/1‑210 (MQ=255)
                                                                                                                                                         gAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCATGGGGGTTGACTGATGACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGt                                                    >  1:9262/1‑251 (MQ=255)
                                                                                                                                                                     gTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCCTAACTGATAATCCTTCAAAATACTTACCCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGtag                                                                                                                                                                     >  1:85263/1‑124 (MQ=255)
                                                                                                                                                                        aaGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAAGTGATGTATATGCTCTCTTTTCt                                    <  1:144481/251‑1 (MQ=255)
                                                                                                                                                                            aCGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGTAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGAc                                 <  1:167561/251‑1 (MQ=255)
                                                                                                                                                                                            cgtgtgctctccgatctttGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGTAAGAACGCGGGATGTTCATGCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACTTTAGTCTCCGACGGc                 <  1:52560/234‑1 (MQ=255)
                                                                                                                                                                                                     ccTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCa         >  1:64862/1‑250 (MQ=255)
                                                                                                                                                                                                      cTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAAt       >  1:7511/1‑251 (MQ=255)
                                                                                                                                                                                                           ggtTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGACAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAAATGCAACTCAATTACTGCAATGCCCTAGTAATTAAGTGAATTTACAATATCGTAATGTTCGTAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGCTGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGAccc  >  1:129946/1‑251 (MQ=255)
                                                                                                                                                                                                                ncNGANCACCANAACTNCTAATCATTCAAACTNCTTAACCTTTGACAGAGCCNNCNCGCAGNctgtc                                                                                                                                                                                     <  2:74622/66‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                
GACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCC  >  NC_000913/579874‑580327

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: