Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,454 T→A Q33L (CAG→CTG)  rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,4540TA97.4% 269.1 / NA 78Q33L (CAG→CTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (34/42);  minor base C (2/0);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGG  >  NC_000913/2867209‑2867699
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTct                                                                                                                                                                                                                                                  <  1:666062/251‑1 (MQ=255)
 cAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTcta                                                                                                                                                                                                                                                 <  1:282208/251‑1 (MQ=255)
 cAGACCACGATTGCCATAACGGCGGGCAATTTTTAACACCCGACGCAAGGTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCCCGCAGTGCGCGACGCGCCAAATAAACTTCTTCTTCGGCCGTTAACAGTGGGGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACAAGCTGTGTGGCTCCCTGCCATAACAGTTCCTCTTCTGCCCAAACGTTATTACTGGGGTCCCGTTct                                                                                                                                                                                                                                                  >  2:168841/1‑250 (MQ=255)
      cACGATTGCCATACCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAg                                                                                                                                                                                                                                            <  2:548061/251‑1 (MQ=255)
        cGATTGCCATAACGGCGGGCAATTTGTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGc                                                                                                                                                                                                                                          <  2:640842/251‑1 (MQ=255)
             gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTGCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCtttt                                                                                                                                                                                                                                     <  2:442432/251‑1 (MQ=255)
                      gcggGCTATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCaaaa                                                                                                                                                                                                                             >  1:555434/1‑250 (MQ=255)
                       cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAAc                                                                                                                                                                                                                           <  1:612687/251‑1 (MQ=255)
                          gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                                                                                                                                                                        <  1:3327/251‑1 (MQ=255)
                             aTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAAc                                                                                                                                                                                                                     <  1:173726/251‑1 (MQ=255)
                               ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTc                                                                                                                                                                                                                   >  1:670946/1‑251 (MQ=255)
                                       ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc                                                                                                                                                                                                            >  1:390041/1‑250 (MQ=255)
                                          gACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                         >  1:355127/1‑250 (MQ=255)
                                            cGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCCCGCAGTGCGCGACGCGCAAAATAACCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGAAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCAGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                                                                                                                                                                                                       >  1:709719/1‑250 (MQ=255)
                                            cGCAAGTTACTCTCGATCAGCCGGCAGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAGAAACCTCAACTCCGTTCTCATCaaa                                                                                                                                                                                                      <  2:174055/251‑1 (MQ=255)
                                             gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                                                                                                                                                                                                     <  2:615181/251‑1 (MQ=255)
                                              cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCAGCAAAtt                                                                                                                                                                                                    <  2:141998/251‑1 (MQ=255)
                                                     ctctCGATCATCCCGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCCAAATAAACTTCTTCTTCGGACGTTAACAGTGGTGAATAACCAATCTCACCCAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCACATCGTTATCACTGGGTTCCAGTTCTCCTAAGGCCTTTTCGTCAAAAACCTCAACTCCGATCTCATCAAATTCCGCATc                                                                                                                                                                                             >  2:246116/1‑251 (MQ=255)
                                                         cGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                                                                                                                                                                            <  1:192657/248‑1 (MQ=255)
                                                         cGATCAACCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGGTCGCATCAAATTCCGCAGCTTCa                                                                                                                                                                                         <  2:524529/251‑1 (MQ=255)
                                                            tcatcCCGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCACAATCAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAttt                                                                                                                                                                                      >  2:678334/1‑251 (MQ=255)
                                                                ccGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAGTTAAAt                                                                                                                                                                                  <  2:237455/251‑1 (MQ=255)
                                                                 cggcggCGAGATGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCGTCGGCCGTTAACAGTGGTTAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCGCAACTCCGTTCTCATCAAATTCTGCATCTTCATTGAAATc                                                                                                                                                                                 <  2:375543/251‑1 (MQ=255)
                                                                  ggcggcGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                                 >  2:560214/1‑250 (MQ=255)
                                                                     ggcgAGAGGCGACATCTCCTCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCGTCATTTAAATCATGa                                                                                                                                                                             <  2:279192/251‑1 (MQ=255)
                                                                     ggcgAGAGGCGACATCTCCACGCAGAGCGCGACGCGCCAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTAGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAATTCATg                                                                                                                                                                              >  2:105155/1‑250 (MQ=255)
                                                                       cgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAAc                                                                                                                                                                           >  2:260173/1‑251 (MQ=255)
                                                                        gagaGGCGACATCTCCACGCAGTGCGCGACGCGAAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACt                                                                                                                                                                          >  1:542306/1‑251 (MQ=255)
                                                                           agGCGACATCTCCACGCAGTGCGCGACGCGCCAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTACAGTTCTACTAAAGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCACCTTCATTTAAATCATGAACTTTc                                                                                                                                                                       >  2:649683/1‑251 (MQ=255)
                                                                               gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTGACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGCGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCCTCTTCATTTAAATCATGAACTTTCAGCg                                                                                                                                                                   >  2:64196/1‑251 (MQ=255)
                                                                               gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCg                                                                                                                                                                   >  2:153702/1‑251 (MQ=255)
                                                                               gACATCTCCACGACGGGCGCGACGCGCAAAATAAACTTCTTCTACCGCCGTTAACAGCGGAGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGACACCCTGCGATAACAGTTCCTCTTCGGCAAAATCGTTCTCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCCAATTCCGCATCTTCATTTAAATCATGAACATTCAGCg                                                                                                                                                                   >  2:433398/1‑251 (MQ=255)
                                                                                 cATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCg                                                                                                                                                                   >  2:239794/1‑249 (MQ=255)
                                                                                 cATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCg                                                                                                                                                                   <  1:239794/249‑1 (MQ=255)
                                                                                      ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                            <  2:547325/251‑1 (MQ=255)
                                                                                        ccGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATGCCGCATCTTCATTTAAATCATGAACTTTCAGCGTaggctgac                                                                                                                                                          <  2:696498/250‑8 (MQ=255)
                                                                                        aCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGAc                                                                                                                                                          <  1:619441/251‑1 (MQ=255)
                                                                                            aGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAt                                                                                                                                                      <  1:345526/251‑1 (MQ=255)
                                                                                                  cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat                                                                                                                                                                                                          >  1:572225/1‑193 (MQ=255)
                                                                                                  cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat                                                                                                                                                                                                          <  2:572225/193‑1 (MQ=255)
                                                                                                  cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg                                                                                                                                                <  1:249702/251‑1 (MQ=255)
                                                                                                   gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCAGCAAATGCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTgg                                                                                                                                               <  1:105155/251‑1 (MQ=255)
                                                                                                         cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                                 >  2:245332/1‑211 (MQ=255)
                                                                                                         cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                                 <  1:245332/211‑1 (MQ=255)
                                                                                                               aaaCTTCTTCTTCGGCCGTTAAAAGTGGTGAATAACCAATCTCACCTCGGCAAAGCTGAGTCGCGTACAACACACGATGTGAGTCTACCTGCGCTAACCGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGATCTACTAAGGCATTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCCTCTTCATTTAAATCATGAACTTTCAGCGAATTCTGCCTCTTAAGGTGGCTCCTACCCGTg                                                                                                                                   >  2:653249/1‑251 (MQ=255)
                                                                                                                       tcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTg                                                                                                                           >  1:724103/1‑251 (MQ=255)
                                                                                                                       tcttcGGCCGTGAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCAGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGTCTCATAAGGTGGCTCCTACCCGTGATCCCt                                                                                                                             <  2:708631/249‑1 (MQ=255)
                                                                                                                          tcGGCCGTTAACAGTTGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGGCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTGTCTTTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCTTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTTCTACCCGTGATCCCTTGa                                                                                                                          <  1:485162/249‑1 (MQ=255)
                                                                                                                          tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGTCTCCTACCCGTGATCCCTTg                                                                                                                           <  1:595338/248‑1 (MQ=255)
                                                                                                                             gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGaa                                                                                                                     <  2:272262/251‑1 (MQ=255)
                                                                                                                             gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGa                                                                                                                      >  1:191743/1‑250 (MQ=255)
                                                                                                                                  tAACAGTAGCGAATAACCACTCTAACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTCGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTCTTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTc                                                                                                                >  2:312922/1‑251 (MQ=255)
                                                                                                                                      aGTGGTGAATAACCAATCTCACCAAGTTAAAGATGAGTCGCGTCCAACACACGCTGAGAGGCTCCCTGCGATAACAGTTCCTCTTCGGCCCAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGAATCTTCATTTAAATCATGAACTTTCAGCGTATTATGACTCATAAGGTGGCTCCTACACGTGATCCCTTGACGGAACATTCCAGc                                                                                                            >  2:633786/1‑251 (MQ=255)
                                                                                                                                             aaTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGcc                                                                                                     >  2:238262/1‑251 (MQ=255)
                                                                                                                                                aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTgg                                                                                                  >  2:84702/1‑251 (MQ=255)
                                                                                                                                                aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTgg                                                                                                  <  2:463104/251‑1 (MQ=255)
                                                                                                                                                aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTgg                                                                                                  <  1:122732/251‑1 (MQ=255)
                                                                                                                                                  ccAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGACTGGtt                                                                                                >  2:417947/1‑251 (MQ=255)
                                                                                                                                                  acAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGAAAAAGCCTGGtt                                                                                                >  2:404841/2‑251 (MQ=255)
                                                                                                                                                             aaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                                                                                      >  1:471766/1‑250 (MQ=255)
                                                                                                                                                             aaGGTAAAGCGGAGCCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAGCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTa                                                                                     <  2:542306/251‑1 (MQ=255)
                                                                                                                                                                          gTCGCGTCCAACACACGCTGCGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCCGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGTTTTATCGCTGCGGCaa                                                                         >  2:514816/1‑250 (MQ=255)
                                                                                                                                                                             gcgTCCAACACACGCTGTGTGGCCCCCTCCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCACAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCTTATTCTTGCTCAAAATGTGGCTCCTACCCGTGATCCCTTGTCGGAACATTCAGGCAAAAGCCGGGTTCCGCCGATTTATCGCTGCGGTAAATaa                                                                     <  2:318411/251‑1 (MQ=255)
                                                                                                                                                                                  cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGc                                                                  >  1:27080/1‑249 (MQ=255)
                                                                                                                                                                                  cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAg                                                                <  2:23474/251‑1 (MQ=255)
                                                                                                                                                                                             gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACgg                                                     <  1:195227/251‑1 (MQ=255)
                                                                                                                                                                                             gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACg                                                      <  1:238262/250‑1 (MQ=255)
                                                                                                                                                                                               gtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAt                                                   >  1:284561/1‑251 (MQ=255)
                                                                                                                                                                                                   cacccTGCGATAACAGTTCCTCTGCGGCCAAATCGTTATCACGTGGTTCCAGTTCTGCTAAGGCCTTTTCGTCAAAACCCTCAACGCCGTTCTCATCAAATCCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACGCAAAAGGCGGCTCTGACCCGTGATCCCTTGACGGAGCATTCAAGCAACAGCCTGGTTCCGCCGCTTTATCGCTGCGGCAAATAACGCAGCGGGTTTCCGGATTTcc                                               <  2:22034/249‑1 (MQ=255)
                                                                                                                                                                                                       cggCGATAACAGTTCCTCGTCTGCCAAATCGCTATCACTGGTGTCCAGTTTTACTAAGTCCTTTTCGTCAAAAACCTCAACTCCGTTCTCAACAAATTCCGCATCTGCATTTAAATCACGAACTTTCAGCGTATTCTGACTCATAGGGTGGCACCTACCCCTGATCCCTTGGCGTTACATTCAAGCAAAAGCCTGGTTCCGCCGAATTAACGCTGTGGCAAATAACGCAGCGGGTTTTCGGATTTCCCCtt                                           <  2:734556/249‑1 (MQ=255)
                                                                                                                                                                                                       cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAACAACCTCAACTCAGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTCTTCTGACTCATAAGGCGGCTCCTACCCGTGACCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGCTTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCtt                                           >  1:736224/1‑251 (MQ=255)
                                                                                                                                                                                                                        ctTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAACAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATACCGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAAt                          >  1:128108/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCa                       >  2:271902/1‑250 (MQ=255)
                                                                                                                                                                                                                               aaaaTCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGGATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTGGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCg                   <  2:555434/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTGTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtg                 <  2:355127/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGaa            <  1:168841/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gTTCCAGTTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCgg  <  2:560202/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCgg  >  2:605892/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGG  >  NC_000913/2867209‑2867699

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