Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,6510CT52.2% 5.7 / 60.7 46intergenic (‑57/‑229)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base T (15/9);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAAT  >  NC_000913/579405‑579895
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ttACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGAATTATCGTTTTTATCGTTTc                                                                                                                                                                                                                                                   >  1:23475/1‑250 (MQ=255)
         ttttAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCt                                                                                                                                                                                                                                         >  1:93286/1‑251 (MQ=255)
           ttAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGGGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTg                                                                                                                                                                                                                                        <  1:304300/250‑1 (MQ=255)
             aaGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGAcc                                                                                                                                                                                                                                     <  1:59861/251‑1 (MQ=255)
             aaGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTGTTTATCGCTTCAATCTGGTCTGAcc                                                                                                                                                                                                                                     <  2:271039/251‑1 (MQ=255)
                          aaCACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttg                                                                                                                                                                                                                        <  1:24654/251‑1 (MQ=255)
                          aaCACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttg                                                                                                                                                                                                                        >  1:432837/1‑251 (MQ=255)
                            cacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttgtt                                                                                                                                                                                                                      <  1:96330/251‑1 (MQ=255)
                                gtgtTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATCCGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTtgatg                                                                                                                                                                                                                  <  1:327093/251‑1 (MQ=255)
                                    tATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGAttt                                                                                                                                                                                                              >  1:43193/1‑251 (MQ=255)
                                          tttCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTGTTGTTGATGATTTATGTc                                                                                                                                                                                                         <  1:458559/250‑1 (MQ=255)
                                                          aaTATTATCTGTTGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTTTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGTAATGttt                                                                                                                                                                                        <  2:521666/251‑1 (MQ=255)
                                                                  ctgtgcttatctAGTAACTATAATGCGAGACATTGTGACGTTTTAGATCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTATAACCTTCCATGTGTTACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTa                                                                                                                                                                                 >  2:81437/12‑250 (MQ=255)
                                                                               gTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATGAATAGTATTGGTTGc                                                                                                                                                                   <  1:463662/251‑1 (MQ=255)
                                                                                    tataATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAAc                                                                                                                                                              >  1:132307/1‑251 (MQ=255)
                                                                                     ataATGTGAGAAGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAAc                                                                                                                                                              >  1:500496/1‑250 (MQ=255)
                                                                                           tgAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGTATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGc                                                                                                                                                       >  2:230994/1‑251 (MQ=255)
                                                                                                cgttGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTcc                                                                                                                                                  <  2:23475/249‑1 (MQ=255)
                                                                                                 aTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCct                                                                                                                                                 <  1:54427/251‑1 (MQ=255)
                                                                                                               gTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATTTGATACGAGGGCGCTTAGTTTGCATTATCTTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGAGTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGa                                                                                                                                   <  1:586560/251‑1 (MQ=255)
                                                                                                                 tCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTACTAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAgg                                                                                                                                 >  1:301563/1‑251 (MQ=255)
                                                                                                                        aaaaCAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATAc                                                                                                                          >  1:164131/1‑251 (MQ=255)
                                                                                                                           acaattccccgactAAATCTTCTCGCACTTGATCGAATATTTCTTTAAAAATCGCAACCTGCGCCATTGGTAATACCTTCCATGTGATACAAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGTTCTGTCCTCATTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAGTGTTTTCACTTAATAGTATTGGTCGCGCAACAAAGGTCGGTCCTGCTGGCATTCTGGAGGGAAATACAAc                                                                                                                       >  1:497491/14‑251 (MQ=255)
                                                                                                                                tCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGCATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGAc                                                                                                                   >  1:175866/1‑250 (MQ=255)
                                                                                                                                 cacaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAACTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGAc                                                                                                                   >  2:82296/1‑249 (MQ=255)
                                                                                                                                                ttCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGcc                                                                                                   >  1:8636/1‑250 (MQ=255)
                                                                                                                                                     aCTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTg                                                                                             >  1:562664/1‑251 (MQ=255)
                                                                                                                                                         gagcGAATATTTCATTAAAAATTGCAACCTGAGCCATGGGTAAAACCTTCCATTTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTGCAATCTGGTCTGACCTCGTTGTGTTTTGTGGATGATGTATGTCGAATATTAGGAATGTTGTCAATGAATAGTACTGGGTGCGTAACAAAGTGCGGGCCTGCTGGCATTCTGTAGGGAAAGACAACCGGCATATTTATGTAAGTCCAATGTGCTCa                                                                                         <  2:309800/248‑1 (MQ=255)
                                                                                                                                                                tatTTCTTTCAAAATGGGAACATGAGACATTGATAAAACCTTCCCTGTGATACGCGGGCGAGAAGTTTGCATTATCGTTTTTCTCACTTCAACAAGGACTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAAGAAAGTGCGGTCCTGATGGCATTCTGGAGGGAAATACAACCGACATATGTATGTAAGGACACAGTGCTCagaactac                                                                                  >  2:560355/1‑244 (MQ=255)
                                                                                                                                                                 atTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTTGTTTTCATTATCTTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGATGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTc                                                                                  <  2:456767/250‑1 (MQ=255)
                                                                                                                                                                      tttAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATAc                                                                              >  1:314726/1‑249 (MQ=255)
                                                                                                                                                                          aaaaTGGCAACCTGAGCCCTTGGTAAAACCTTCCATGTCATAAGAGTGTGCGCAGTTTGCATTATCGATTTTATCGTTTCAATCTGGTCTGTCCCCCTTGTGTTTTGTTGATGTTTTATGTCAAATATTAGTAATGTTTTCACTTAATAGTATTGTTTGATTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGGTGTATGTAAGGCCACCGTGCTCaattaatcatacagaaa                                                                         >  2:542323/1‑235 (MQ=255)
                                                                                                                                                                          aaaaTGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGaa                                                                          >  1:348470/1‑249 (MQ=255)
                                                                                                                                                                            aaTGGCAACCTGCGCCATTGGTACAACCTTCCATGTGATACGAGGGCGCGTCGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATCTTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACAGCCAGATGTATGTAAGGCCAACg                                                                                               >  2:366320/1‑226 (MQ=255)
                                                                                                                                                                            aaTGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTGTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATAGTAGGAATGTTTTCACTTAATAGTATGGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACg                                                                                               <  1:366320/226‑1 (MQ=255)
                                                                                                                                                                                 cAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGa                                                                  >  1:36464/1‑250 (MQ=255)
                                                                                                                                                                                                    aaCCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAg                                              <  2:43193/251‑1 (MQ=255)
                                                                                                                                                                                                            aTGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAag                                        >  2:197657/1‑249 (MQ=255)
                                                                                                                                                                                                                 aTACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAgc                                 >  1:310548/1‑251 (MQ=255)
                                                                                                                                                                                                                               agtttgaagtgttgTTTTGATCGTTTCAATCTGGTCTGACCTCATTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATTTTTACAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCCTTCTGGAGGGAAATACAACCGACAGATGTATGTAAGACCAACGTGCTCAAACCTGCATGCAGCAAGATTTGAAGCAGTATTTTAAGCGCTAGATGAAGAGCAAGCGCATGGAGCGACa                    <  2:213875/237‑1 (MQ=255)
                                                                                                                                                                                                                                   tGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATg               >  1:2077/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ttATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAACATGAAt            >  1:145954/1‑250 (MQ=255)
                                                                                                                                                                                                                                             tttttATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATTTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGa                          <  2:522972/230‑1 (MQ=255)
                                                                                                                                                                                                                                             tttttATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGa                          >  1:522972/1‑230 (MQ=255)
                                                                                                                                                                                                                                              ttttATCGTTTCAATCTTGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAAc     <  2:432837/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ttATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAAt  >  1:202863/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
TTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAAT  >  NC_000913/579405‑579895

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