Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,6510CT51.7% 0.4 / 37.3 29intergenic (‑57/‑229)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
Reads supporting (aligned to +/- strand):  ref base C (6/8);  new base T (7/8);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCT  >  NC_000913/579425‑579900
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
tttATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGATACATTGTGACATTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCGTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCttt                                                                                                                                                                                                                                     <  1:3233/249‑1 (MQ=255)
         acGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTttgttg                                                                                                                                                                                                                          <  1:7802/251‑1 (MQ=255)
                                                   gCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGACGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGt                                                                                                                                                                                 >  1:138414/1‑250 (MQ=255)
                                                                 ataATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAAc                                                                                                                                                                   <  1:153016/250‑1 (MQ=255)
                                                                 ataATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTgg                                                                                                                                                                                                                                                 <  1:247868/172‑1 (MQ=255)
                                                                        gagaCGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTNc                                                                                                                                                            <  1:254294/250‑1 (MQ=255)
                                                                                             tCAGANTAAAACAATTCACAGTCTAAATCTTTTCGCANTNGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTNCCATGTGATANGAGNNCGCGNANTNTGCATTATCGTTTTTATCGTTTCAATCTGNTCTGACCTCCTTGTGTTTTGTTGNNGANNNATNNCAAat                                                                                                                                                                                                          <  2:242929/183‑1 (MQ=255)
                                                                                             tCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTACAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGCGCTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAgg                                                                                                                                      >  1:141256/1‑251 (MQ=255)
                                                                                             tCAGAATAAAACAATTCACAGTTTAAATCGTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGGGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTGTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAg                                                                                                                                       <  1:232151/250‑1 (MQ=255)
                                                                                             tCAGAATAAAACAATTCACAGTCTAAATCTTTTCGTAATTGATCGAATATTTCTTTAACAACAACAACCTGCGCCATTCGTAAAACCTTACATGAGATAAGAGGGCGCGTAGTTTGCATTATCGTTTTTCTCGTTTCAACGTGATCTGACATCCTTGTGTTTTGTTGATGATTTATGtaaact                                                                                                                                                                                                          >  1:242929/1‑178 (MQ=255)
                                                                                             tCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCGGCTGGCATTCTGGAg                                                                                                                                       <  1:110152/250‑1 (MQ=255)
                                                                                              cAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGtttt                                                                                                                                                                                            >  1:258300/1‑196 (MQ=255)
                                                                                                     aaaCAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGATCCATTGGTAAAACCTTCCATGTTATACGAGGGCGCGTAGTTTGCATTATCGTTTTGATCGCTTCAATCTGGTCGGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAg                                                                                                                                                                                                     <  1:209527/180‑1 (MQ=255)
                                                                                                       aCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAACCATGGCACCATGAGCCATTGGTAAA‑CCTTCCACGTGATAAGAGGGCGCGTAGTTTGCATTCTCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATAATAGGCATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGc                                                                                                                                                >  1:312077/1‑230 (MQ=255)
                                                                                                                  tctAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTAGCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCGAATATTACGAATGTTTTCACTTAATAGTATTGGTTGTGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATAc                                                                                                                               <  1:128063/235‑1 (MQ=255)
                                                                                                                   ctAAATCTTTTCACACTTGATCGCATATTTCTTTAAAGATCGCAACCTGCGCCATTGGTAAACACTTCCATGTGATACGAGGGCGAGAAGTTTGCCTTATCGTTTATATCGTTTCCATCTGGTCTGACCTCCTTGTGTTTTGTTGTTGTTTTATGTCAAATATTAGGACTGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGc                                                                                                                                                            >  1:318461/2‑207 (MQ=255)
                                                                                                                           tttCGCACTTGACCGAATATTTCTTTAAAAATGGCAACCTGAGCAATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTTTGTCAAATATTAGGAATTTTGGCAAGTAATAGTATTGTTTGCGTAACAAAGTGCGGTCCTGCTGGCATCCTGTAGGGAAATACAACCGACAGATGTATGTAAg                                                                                                           <  1:133270/248‑1 (MQ=255)
                                                                                                                               gCACTTGCTCGAATATTTCTTTACACATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCCTTATCGTTTTTATCGCTTCACTCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGTTTGCGTAACAAAGTGCGGTCctgctg                                                                                                                                                  >  1:261080/1‑205 (MQ=255)
                                                                                                                                        cGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATTTGAGACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGan                                                                                                                                                                                                                                           <  1:32607/107‑2 (MQ=255)
                                                                                                                                                   ttAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGt                                                                                                                                                                                               <  1:17523/140‑1 (MQ=255)
                                                                                                                                                            gCAACCTGAGCCATTGGTAAAACCTTNCATGTGNNACNAGGGCGCGTAGTTTGCATTNTCNNNTNNNTNGCTTCNATNNGGTCTGACCNCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATNGTATTGGTTNNGTAACAAAGTGCGGTCCTGCTNGCANNCNNGNGGGnaa                                                                                                                                  >  2:68187/1‑192 (MQ=255)
                                                                                                                                                            gCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAACTATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAg                                                                                                           >  1:228285/1‑215 (MQ=255)
                                                                                                                                                                     gCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGCAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTATGCCCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAAAAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCa                                                                                                       >  1:214775/1‑210 (MQ=255)
                                                                                                                                                                                   cTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGc                                                                                                         <  1:12497/194‑1 (MQ=255)
                                                                                                                                                                                     tCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGGCAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGa                                               <  1:133257/250‑1 (MQ=255)
                                                                                                                                                                                         tgtgATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTc                                                                                                                                            >  1:134912/1‑153 (MQ=255)
                                                                                                                                                                                          gtgATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTAAAATATTCGGAATGTTTTCAATTAATAGTCTTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCTTTCTGGAGGGACATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGc                                          >  1:324459/1‑250 (MQ=255)
                                                                                                                                                                                                                gCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCCAATATTAGGAATGTTTTCCCTTAATAGTATTGGTTGCGTAAAAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGAAAGCGCATGGAGCGACAAAATg                    >  1:232251/1‑250 (MQ=255)
                                                                                                                                                                                                                     aTCGTTTTTATCGTTTCAATCTGGTCTGACCGCCTTGTGTTTTGTTGATGATTTATGTCAAATATTan                                                                                                                                                                                                     <  1:84287/68‑2 (MQ=255)
                                                                                                                                                                                                                                 ncTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGATTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATctgct  >  1:196039/2‑251 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
TTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCT  >  NC_000913/579425‑579900

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