Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913292,4330GA54.8% ‑5.8 / 58.6 42F172F (TTC→TTT) yagKuncharacterized protein YagK
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base A (13/10);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGG  >  NC_000913/292191‑292676
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
aaaCGTTAGTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATGTTTTGAAACAAGACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTGACGGTCTCCTTCACCAAACACTTTGCGTTCAACTGTGGGCAGGAAGTCTAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGt                                                                                                                                                                                                                                             <  1:191403/251‑1 (MQ=255)
  aCGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTAc                                                                                                                                                                                                                                           <  1:63292/251‑1 (MQ=255)
    gTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACtt                                                                                                                                                                                                                                         >  2:94101/1‑250 (MQ=255)
     ttATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTAt                                                                                                                                                                                                                                        >  1:100469/1‑251 (MQ=255)
               gCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAg                                                                                                                                                                                                                              <  2:275042/251‑1 (MQ=255)
                          gatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTa                                                                                                                                                                                                                    <  1:155404/250‑1 (MQ=255)
                                   cAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCa                                                                                                                                                                                                           >  2:77774/1‑250 (MQ=255)
                                   cAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCa                                                                                                                                                                                                           <  1:77774/250‑1 (MQ=255)
                                   cAGTAATCTATAGAAAATTGAGTAACGAAACTGCATTTTAGTAATTTTTTGACACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTc                                                                                                                                                                                                            >  2:302988/1‑249 (MQ=255)
                                                          aaCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTg                                                                                                                                                                                    >  2:292373/1‑250 (MQ=255)
                                                             gAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCg                                                                                                                                                                                <  1:129651/251‑1 (MQ=255)
                                                              aaaCTGTATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCAATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTAAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCCGGGTGATCGt                                                                                                                                                                               >  2:50555/1‑251 (MQ=255)
                                                               aaatgaaatttAGTAATATTTTGAAACAAAACAAATAAGGCCTTCTCACCGAGAAGGCAATTACCGACTTACAGTTCGATTAGGCTACAGACGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTAAACTTTGGTCAGGTAGTACAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCAAATACTTACAGTTTTCAGGAAAATGGACAAGACCAGCGTGATCgac                                                                                                                                                                              >  2:274438/9‑249 (MQ=255)
                                                                   gCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACGTTGGTCAGGTATGCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGtctctc                                                                                                                                                                          <  1:7773/251‑1 (MQ=255)
                                                                               ttttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAgg                                                                                                                                                              >  1:190416/1‑251 (MQ=255)
                                                                                  tttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGgc                                                                                                                                                            <  2:138815/250‑1 (MQ=255)
                                                                                      aaaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACGTACAGTTTTCAGGAAAATGGACCAGACCAGGGGGATCGTCTCTCTCCAGTCTCAGGGCGCTGt                                                                                                                                                       <  2:183441/251‑1 (MQ=255)
                                                                                        aCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTAc                                                                                                                                                     >  1:85677/1‑251 (MQ=255)
                                                                                          aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTAcc                                                                                                                                                    >  1:36553/1‑250 (MQ=255)
                                                                                          aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCa                                                                                                                                                   <  1:79980/251‑1 (MQ=255)
                                                                                                           ctcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATcc                                                                                                                                  >  2:141176/1‑251 (MQ=255)
                                                                                                             cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATccc                                                                                                                                 <  1:21907/250‑1 (MQ=255)
                                                                                                                   gaAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAgg                                                                                                                          >  1:235861/1‑251 (MQ=255)
                                                                                                                        ccATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATc                                                                                                                     >  2:198603/1‑251 (MQ=255)
                                                                                                                             cccGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCt                                                                                                                <  2:78752/250‑1 (MQ=255)
                                                                                                                                 aCTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATc                                                                                                            <  1:282583/251‑1 (MQ=255)
                                                                                                                                      cAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAAt                                                                                                       <  2:300867/251‑1 (MQ=255)
                                                                                                                                      cAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGGTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCGCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAAt                                                                                                       <  2:300846/251‑1 (MQ=255)
                                                                                                                                                  gCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtggt                                                                                           >  1:216261/1‑251 (MQ=255)
                                                                                                                                                          cGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCg                                                                                   >  1:38247/1‑251 (MQ=255)
                                                                                                                                                          cGAAGATACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGACAGGAAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGCGGAAGTAAGCg                                                                                   >  1:282163/1‑251 (MQ=255)
                                                                                                                                                                                            ttCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTgg                                                  <  2:235861/250‑1 (MQ=255)
                                                                                                                                                                                                     ggTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCATTCTCAGGGCTCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                        <  1:50555/251‑1 (MQ=255)
                                                                                                                                                                                                        cAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAt                                      <  1:61317/250‑1 (MQ=255)
                                                                                                                                                                                                              gTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAGTCAGAAt                               <  2:30867/251‑1 (MQ=255)
                                                                                                                                                                                                                 cAGTCGGGTGAACAGCGGCTGGCAGTTCTGCTGAAAGTGAGTTCTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGCACCAGACCAGGGTGGTCGTCTCCCTCCAGGCTCAGGGCGCTCTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGGCTTTGTTGAACTACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAAtgcg                            <  2:162602/251‑4 (MQ=255)
                                                                                                                                                                                                                               gcgcCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCatat               >  2:128066/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gcCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat            <  1:252218/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTAGTTCTGCTGAAAGTCAGTACTATTGGTATCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAACTGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaat        >  2:254137/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aGTTCTGCTGGAAGGCAGTACTATTTGTGTCTAGCACATACTTACAGTTTGCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGTCTCCTGTGATCAGCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTTATTCAGAATACTCTTGCGCCCATATAacaaataa      <  2:111466/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATActaaataa      >  2:45250/1‑249 (MQ=255)
                                                                                                                                                                                                                                           tCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAggg  >  1:215740/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
AAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGG  >  NC_000913/292191‑292676

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