Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,799 A→G N43N (AAT→AAC)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,7990AG95.4% 267.6 / 0.3 88N43N (AAT→AAC) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (35/48);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

ACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATA  >  NC_000913/72552‑73042
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACCGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTc                                                                                                                                                                                                                                                   >  2:274904/1‑251 (MQ=255)
aCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTc                                                                                                                                                                                                                                                   >  1:618288/1‑251 (MQ=255)
  cgcgCCATTAAATTATCAATCCCCATCTGGCGGGCGATTGCGTGCATTGTCGCATGCGGTTCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCCGCCCGATGTACTGCGCGACCGTCAGGCGACTGGACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCGGTCGTCAGCAACCCGGCGTTCAAGTTcag                                                                                                                                                                                                                                                 <  2:599803/251‑1 (MQ=255)
  cgcgCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTcag                                                                                                                                                                                                                                                 <  1:498405/251‑1 (MQ=255)
      ccATTAAATTCTCAATCCCCATCGGGCGGGCGATAGCGTGCATTTTCTCCTGCTTTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCACAGCCCGATGTTCTTCGCGACCTGCAGTTGACTGAACAGGTTGTTTTCCTGAAACAGCACCGACACCGGACGGCGTGCCGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTcagcatg                                                                                                                                                                                                                                             <  2:215203/251‑3 (MQ=255)
       cATTAAATTATCAATCCCCATCTGGCGTGCGATAGCGTTCACTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGTTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAggg                                                                                                                                                                                                                                            <  2:473237/251‑1 (MQ=255)
          tAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTa                                                                                                                                                                                                                                          <  1:505551/250‑1 (MQ=255)
            aaTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACtt                                                                                                                                                                                                                                       <  1:22573/251‑1 (MQ=255)
             aTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGCTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACttt                                                                                                                                                                                                                                      >  1:83352/1‑251 (MQ=255)
              ttATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACtttt                                                                                                                                                                                                                                     >  1:653675/1‑251 (MQ=255)
              ttATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAAAGCCGTCGATAGGCAGCGAACCGCTGGCTGGGTTCGGAAAAACAGCGAACGAGCTCAGCAGGGTACtttt                                                                                                                                                                                                                                     >  2:600580/1‑251 (MQ=255)
                  cAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAcc                                                                                                                                                                                                                                  <  1:229188/250‑1 (MQ=255)
                      ccccTTCTGGCGGGCCATAGCGTGCATTTTCTCCTGCTGGGCCGCGTTCAGTTTCAATCCCGGATTGAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACGGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGAGAGCAAGTTCAGCAGGGTACTTTTACCGg‑cgc                                                                                                                                                                                                                            <  2:244066/251‑1 (MQ=255)
                       cccATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGGCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑cgcc                                                                                                                                                                                                                           >  2:328623/1‑251 (MQ=255)
                          agcTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCt                                                                                                                                                                                                                        <  1:637950/249‑1 (MQ=255)
                             tGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTgg                                                                                                                                                                                                                     >  1:647474/1‑251 (MQ=255)
                             tGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAAAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTCCCTG‑CGCCGCTTgg                                                                                                                                                                                                                     >  2:664217/1‑251 (MQ=255)
                              ggcgggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTg                                                                                                                                                                                                                      <  1:87633/249‑1 (MQ=255)
                                     gATAGCGTGCAGTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGa                                                                                                                                                                                                             <  1:103855/251‑1 (MQ=255)
                                           gTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGcc                                                                                                                                                                                                        >  1:164242/1‑250 (MQ=255)
                                               aTTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTg                                                                                                                                                                                                   <  1:161604/251‑1 (MQ=255)
                                               aTTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTg                                                                                                                                                                                                   <  1:401387/251‑1 (MQ=255)
                                                     cccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTcgcc                                                                                                                                                                                             <  2:647474/251‑1 (MQ=255)
                                                      ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCCAGATCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCCCCTGCTcgccg                                                                                                                                                                                            >  2:73713/1‑251 (MQ=255)
                                                      ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTcgccg                                                                                                                                                                                            >  1:26406/1‑251 (MQ=255)
                                                            gtgccGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGCGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTct                                                                                                                                                                                      >  1:447914/4‑250 (MQ=255)
                                                                  cgTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGAAGGCGTAACGGCGGCGTCGTGGTGTGATAAACGCCGTCGAGAGTCAGCGGACCGCTGGCTGTCGTCAGAAAACCGACGATCAAATTCAGAAGGGTACTTTTACCGG‑CGCCTCTTAGCCCGAGGATCGCCACCTGCACGCCGCGTTCTACCGtt                                                                                                                                                                                >  2:467399/1‑251 (MQ=255)
                                                                   gTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCATGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCCAGATCAGAAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGTGCGCCACCTGCTCGCCGCGTTCTACCggt                                                                                                                                                                                >  2:532623/1‑248 (MQ=255)
                                                                     tCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTa                                                                                                                                                                               >  1:391171/1‑249 (MQ=255)
                                                                       cgttacaaccccGGATTTAGCCCCAGCCCGAGGGTCTGTGCGACCGTCAGGTGGCTGACCAGGTTGTCCTCCTGAAACATCATCGACACCTGACGGCGGTACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCTCGAGGATCGCCTCCTGCTCGCCGCGTTCCACCTTTAAGCt                                                                                                                                                                           <  2:55711/242‑1 (MQ=255)
                                                                              aTCCCGGATTTAGCCCCAGCCCGATTTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGc                                                                                                                                                                    <  2:618288/251‑1 (MQ=255)
                                                                                cccGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCAt                                                                                                                                                                  >  1:120129/1‑251 (MQ=255)
                                                                                 ccGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATc                                                                                                                                                                 <  1:344071/251‑1 (MQ=255)
                                                                                  cGGATTTAGCCCCATCCCGATGTTCCGCGCGCCCGTCAGGTGACGGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGTCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTG‑CGCCGCTGGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCTACCGGTAATCTAAAACGCATCg                                                                                                                                                                <  2:311618/251‑1 (MQ=255)
                                                                                   ggATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCgg                                                                                                                                                               <  1:607253/251‑1 (MQ=255)
                                                                                       ttAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCaa                                                                                                                                                            <  1:418436/250‑1 (MQ=255)
                                                                                         aGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCGGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtg                                                                                                                                                         <  2:140477/251‑1 (MQ=255)
                                                                                            cccATCCCGATGTTCTGCGCGACCGTCAGGTGAATGAACAGGTTTTTTTCCTGCAACAGCATCGACACCGGTCGGCGTGACGGCGGCGTCGTGTTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGGCAGAAAACCGGCGATCAATTTCAGCAGGGTACTTTTTCCTG‑CGCCGTTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttg                                                                                                                                                      <  2:566792/251‑4 (MQ=255)
                                                                                                gCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTaa                                                                                                                                                  <  1:328623/251‑1 (MQ=255)
                                                                                                    gATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGcc                                                                                                                                              >  1:322304/1‑251 (MQ=255)
                                                                                                         tCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGt                                                                                                                                          <  2:164242/250‑1 (MQ=255)
                                                                                                              gCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAgggggtatt                                                                                                                                    >  2:177651/1‑244 (MQ=255)
                                                                                                               cGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtca                                                                                                                                   <  1:274140/251‑1 (MQ=255)
                                                                                                               cGACCGTCAGGTGACTGAACAGGTTGTTTTCCGGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCTGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtca                                                                                                                                   <  1:664217/251‑1 (MQ=255)
                                                                                                                    gTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtcagtcag                                                                                                                              <  2:337658/251‑1 (MQ=255)
                                                                                                                         gtgacTGAACGGGTTGTTTTCCTGAATCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGTACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGGTCAGCAGGGTACTTTTACCTG‑CGCCGCTGGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTa                                                                                                                         <  2:588556/247‑1 (MQ=255)
                                                                                                                          tgacTGAACAGGTTGTTTTCCTGAAACAGCATCTACACCTGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                     <  2:320787/219‑1 (MQ=255)
                                                                                                                          tgacTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                     >  1:320787/4‑222 (MQ=255)
                                                                                                                            acTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAAc                                                                                                                       >  1:562290/2‑250 (MQ=255)
                                                                                                                             ctgtaCAGGTTGTTTTCCTGAAACAGCATCGACACAGGACGGCGTGAAGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGCACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAACAGGCTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTACCCGGGATCTCCCGTGAAACACAAACTCATCGCCAAATGTTGGTAAAGCCGGGTGATATACGCAAGTTTTAACAt                                                                                                                     >  2:511308/5‑251 (MQ=255)
                                                                                                                               gAACAGGTTTTTTTCCTGAAACAGCAGCGACACCGGACGGCGTGACGGCGGCGGCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAtc                                                                                                                    <  2:391171/250‑2 (MQ=255)
                                                                                                                                    ggTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGAGTGATCACCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCCTCGGCAAATGTTGGTAAAGCCCAGTGGTATCAGTCAGTTTTAACATTTCGccc                                                                                                              >  2:425865/1‑251 (MQ=255)
                                                                                                                                    ggTTGTTGTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGccc                                                                                                              <  1:291200/251‑1 (MQ=255)
                                                                                                                                     gttgttTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGcccc                                                                                                             <  1:44122/251‑1 (MQ=255)
                                                                                                                                          tttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTa                                                                                                         <  1:207164/248‑1 (MQ=255)
                                                                                                                                           ttcCTGAAACATCATCGACACCGGACGGAGTGACGGCGGGGTCGTGGTGTGAACAACGCCGTCGATAACCAGCGATCCGCTGGCTGGCGCCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGttg                                                                                                        <  2:524112/249‑3 (MQ=255)
                                                                                                                                           ttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAg                                                                                                        <  1:83965/249‑1 (MQ=255)
                                                                                                                                                   ccAGCATCGACACAGGACGGCGTGCCGGCTGCGTCGTGGTGTGGTCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGAACAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTTAGCTAAAACGCATCGGTCAATGTTGGTAAAGCCAGGAGCTATCAGTCAGTTTTAACCTTTCGCCCCGGTAATTTTTccac                                                                                               >  2:301675/2‑248 (MQ=255)
                                                                                                                                                   aCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAAt                                                                                               >  1:699840/1‑251 (MQ=255)
                                                                                                                                                    cAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGAAATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATc                                                                                              >  1:234547/1‑251 (MQ=255)
                                                                                                                                                     aGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAATAAGCGAACCGCTGGCTGGCGTCAGAAAACAAGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGATCGACGCGTTCTACCGTTAAGCTAAAACGCATCGGCAACTGTTGGTAAAGACAGGTGATATCAGTCAGGTTTAACATCTCGCCCCGGTAGTTTTTCAAta                                                                                              >  1:703972/1‑249 (MQ=255)
                                                                                                                                                        atcgaacccgtccGGCGTGACGGCGCGATCGTGGTGTGATCAAGGCCGTCGATAGTCAGCGAGCCGCGGGCTGGCGTCATAAACCCGTCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCGTGCCGCGTGCTACCGTTAAGCGAAAACGCAGCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTGTTCAATCACgg                                                                                          <  2:491846/239‑1 (MQ=255)
                                                                                                                                                              cccGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAca                                                                                    <  1:530579/250‑1 (MQ=255)
                                                                                                                                                               ccGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcag                                                                                   <  2:454654/251‑1 (MQ=255)
                                                                                                                                                               ccGGACGGCGTGACGGCGGCGTCGTGGTGTGAGCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTGCAATCACGGTAAAcag                                                                                   <  1:219768/251‑1 (MQ=255)
                                                                                                                                                                   aCGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagc                                                                                  >  1:17259/1‑248 (MQ=255)
                                                                                                                                                                                     cgtggtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGCTCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagcagc                                                                >  2:366371/5‑248 (MQ=255)
                                                                                                                                                                                         gtgtgATCAACGCCGTCGATAGTCAGCGCACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                          >  2:220108/1‑250 (MQ=255)
                                                                                                                                                                                         gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACCGAGTAGCAGCAGAAtt                                                          >  2:140136/1‑250 (MQ=255)
                                                                                                                                                                                         gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCCGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                          >  1:653924/1‑250 (MQ=255)
                                                                                                                                                                                          tgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATAGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTa                                                         >  1:469040/1‑250 (MQ=255)
                                                                                                                                                                                           gtgATCAACGACGTCGATAGTCAGCGAACCGATGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCGCGGTAGTGTTTCAATCACGGTAAACAGTAGAAAACAGAGTAGCAGCAGAATTAAc                                                       <  1:237695/251‑1 (MQ=255)
                                                                                                                                                                                                caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                  <  1:513750/249‑1 (MQ=255)
                                                                                                                                                                                                caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTGACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                  <  2:83218/249‑1 (MQ=255)
                                                                                                                                                                                                 aaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg                                                 >  1:702875/2‑251 (MQ=255)
                                                                                                                                                                                                  acgccggcGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGa                                                <  1:563089/244‑1 (MQ=255)
                                                                                                                                                                                                           aTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGt                                       >  1:30549/1‑251 (MQ=255)
                                                                                                                                                                                                             aGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                     >  1:209714/1‑139 (MQ=255)
                                                                                                                                                                                                             aGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                     <  2:209714/139‑1 (MQ=255)
                                                                                                                                                                                                              gTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCt                                    <  1:486554/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCt                                <  2:653924/251‑1 (MQ=255)
                                                                                                                                                                                                                        cGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGAt                           >  1:67130/1‑250 (MQ=255)
                                                                                                                                                                                                                           gggctggcGGCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTGCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGTCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCTGTCCTGGCTGCGATAGGa                       <  1:156338/250‑1 (MQ=255)
                                                                                                                                                                                                                            ggctggcGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAg                      <  1:73713/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCGCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAtt                <  1:309548/249‑1 (MQ=255)
                                                                                                                                                                                                                                       gAAAACCGGCTATCAAGTTCAGCAGGGTACTTG‑ACCCGGCGCCGCTTGGCCCGAGGATCGCTACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGTTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTg           <  2:26406/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aCCGGCGAGCAAGTTACGCAGGGTACTTTTCCCTG‑CGCCGCTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTag       <  2:93357/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gATCAAGTTCAGAAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCAACTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACAGCGCCGTACTGGCTGCGTTAGGAGCCAATTTGCTGGTAGAGATa  >  1:219478/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATA  >  NC_000913/72552‑73042

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