Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,642,358 T→A 5.8% intergenic (+27/+22) dtpB → / ← rsmJ dipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,642,3580TA 5.8% 368.1 / 9.1 104intergenic (+27/+22)dtpB/rsmJdipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (46/52);  new base A (3/3);  total (49/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

Only 100 of 105 total aligned reads displayed.

TGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAG  >  NC_000913/3642084‑3642625
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                           
tGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGt                                                                                                                                                                                                                                                                             <  2:362556/275‑1 (MQ=255)
          tttGATGGGCTTGATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTGGCCCGTCTATACCAACGTGGTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAcc                                                                                                                                                                                                                                                                   <  2:131382/275‑1 (MQ=255)
             gATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCCACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAccc                                                                                                                                                                                                                                                                  >  2:104563/1‑273 (MQ=255)
                gggTTTTTTTCTGGGGTGGTGTTTCCTCGNNCATGCTGCGCCGTTCTTGCTTGGAGGCTATNTTGCAAACTTTACTGCGGTGCCGGACAACATTACTGCTCCGCTTTTGACGTTGCCCGTCTATACCATGGTGTTTNTTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                                                                                                                                             <  3:63584/275‑1 (MQ=255)
                      tttgctcGGGGTTTGGTTCCTGACGCGGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTAACGATCCGCTTGAGACGTGGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtt                                                                                                                                                                                                                                                        <  4:262128/270‑1 (MQ=255)
                          atagtagacaggTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTac                                                                                                                                                                                                                                                    <  1:112220/264‑1 (MQ=255)
                            cGGGATGTGGTTCCTGACGCAGGCTGCCGCGTCCTTGCTGGGCGGCTATGTGGCAACATTTCCCGCGGTGCCGGACATCATTACCGACCCTCTTGAGACGTTGCCCGTCTATACCAACGTGCTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCCGTTTACGGTTCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTGTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                                                                                                                                                 >  2:283386/1‑275 (MQ=255)
                             gggATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTg                                                                                                                                                                                                                                                >  2:178726/1‑275 (MQ=255)
                               gATGGGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTgg                                                                                                                                                                                                                                               <  2:181597/274‑1 (MQ=255)
                                         ctgccgcgtgtggCCGCGTTTTTGCTGGGGGGGTAGTTGGCAACATTTACCGCGGTGCCGGACAACATTACCGAGCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCAt                                                                                                                                                                                                                                    <  4:104562/263‑1 (MQ=255)
                                         cTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCAt                                                                                                                                                                                                                                    <  2:53654/275‑1 (MQ=255)
                                           gtcggaGGCTGCCGCGTTTTTGTTGGCGGGTTATTTGGCAACATTTACCGGGGTGCCGGAAAACTTTACCGATCCGCGTGGGATGTGGTCCGTCTATACCAAGGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTGGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATa                                                                                                                                                                                                                                   <  4:191412/269‑1 (MQ=255)
                                              gCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGACAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGa                                                                                                                                                                                                                                 >  2:10671/1‑273 (MQ=255)
                                                            tCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGcc                                                                                                                                                                                                                  >  1:90658/1‑274 (MQ=255)
                                                                                   aCATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcg                                                                                                                                                                                          >  2:338136/1‑275 (MQ=255)
                                                                                               gTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCg                                                                                                                                                                              >  2:249068/1‑275 (MQ=255)
                                                                                                   aggACGACATGACCGATCCGCTTGAGAGGTTATCCGTCTATACCATCGTGTTTGTTAAGATTGGTCTGGGCCCGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAg                                                                                                                                                                          <  1:94217/274‑1 (MQ=255)
                                                                                                                    ccGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAg                                                                                                                                                                          >  2:190071/1‑258 (MQ=255)
                                                                                                                    ccGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCg                                                                                                                                                         >  2:176449/1‑275 (MQ=255)
                                                                                                                      ggttGAGACGTTGCCCGTTGATACCAACGTGTTTGGTAATATTGGTCTGGGAACGCTGTGCGTTGGAGTAGGGATGCTGTTGATGGGGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                                                                                                                                       <  4:384825/273‑1 (MQ=255)
                                                                                                                      gctatagactatGCCCGGCTATACCAAAGTGTTTTGTAAGATGGGTCTGGTGACGCGGGGCGTTGCAGTAGGGATGCTGTGGATGGGGCCGTGGCTGAAACGCAGGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                                                                                                                                       <  4:338132/264‑1 (MQ=255)
                                                                                                                        ttGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                                                                                                                                                      <  4:190069/258‑1 (MQ=255)
                                                                                                                                  gCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCAATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATATATATCAAAGCGATGCTCTTTAGTGACTACCGCGTTTGGCGTTGCGACATTCGCAAGCGGTGGCACGAAGTCCGGCCGCTTTaccac                                                                                                                                            >  4:292737/1‑274 (MQ=255)
                                                                                                                                    ccGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGc                                                                                                                                                                                                                                                                       >  6:377468/1‑149 (MQ=255)
                                                                                                                                           taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCt                                                                                                                                                                                                                                                                 >  6:93177/1‑148 (MQ=255)
                                                                                                                                           taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCt                                                                                                                                                                                                                                                                 <  5:93177/148‑1 (MQ=255)
                                                                                                                                           taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTGAATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTg                                                                                                                                                                                                                                                                >  8:464007/1‑149 (MQ=255)
                                                                                                                                             ccAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtt                                                                                                                                                                                                                                                              <  5:286052/149‑1 (MQ=255)
                                                                                                                                             ccAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtt                                                                                                                                                                                                                                                              <  7:25731/149‑1 (MQ=255)
                                                                                                                                               aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                                                                                                                                            >  5:78121/1‑149 (MQ=255)
                                                                                                                                               aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                                                                                                                                            >  8:456462/1‑149 (MQ=255)
                                                                                                                                                aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgc                                                                                                                                                                                                                                                               >  7:282330/1‑144 (MQ=255)
                                                                                                                                                aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgc                                                                                                                                                                                                                                                               <  8:282330/145‑2 (MQ=255)
                                                                                                                                                 cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGTCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGc                                                                                                                             >  2:252677/1‑274 (MQ=255)
                                                                                                                                                    gtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagttt                                                                                                                                                                                                                                                       <  6:78121/149‑1 (MQ=255)
                                                                                                                                                     tttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCt                                                                                                                                                                                                                                                                 >  2:84295/1‑138 (MQ=255)
                                                                                                                                                         gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                                                                                                                                             <  4:84294/138‑1 (MQ=255)
                                                                                                                                                           aaGATTGGTCTGGTCACGCTGGGCGTTGCAGGAGTGAGGCGGTTGATGGTGCCGTGGCTGAAACGCAGGATTGCGACGACGGAAAGCCATTAATTAGTCTTGCAGAAAGCAGGGTAGCGTGATCTCTACCCGGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCTGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTg                                                                                                                  <  2:282073/275‑1 (MQ=255)
                                                                                                                                                           aaGATGGGTCTGGTACCGTTGGGCGTTGCAGTAGGGAGGCTGTTGATGGTGCCGGGGCTGAAACGCATGATTGGAACACCGGAAAGAGATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACGCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCGCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTg                                                                                                                  <  4:275933/275‑1 (MQ=255)
                                                                                                                                                             gATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATTCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCCCCCCGCGTTTTGTCGCCAGCAAGccgg                                                                                                                  >  2:485261/1‑270 (MQ=255)
                                                                                                                                                             gATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCtt                                                                                                                                                                                                               >  2:521949/1‑180 (MQ=255)
                                                                                                                                                              aTTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAg                                                                                                               <  2:462880/275‑1 (MQ=255)
                                                                                                                                                                tggtttggtcAGGTGGGGGGTTGCAGTAGGGGTGCTGTTGATGGTGCCGTGGCTGAAACGGATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCTCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAgg                                                                                                              <  4:249064/274‑1 (MQ=255)
                                                                                                                                                                 ggtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGt                                                                                                                                                                                                           <  4:521944/180‑1 (MQ=255)
                                                                                                                                                                  gtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTAcc                                                                                                                                                                                                                                         >  7:163707/1‑149 (MQ=255)
                                                                                                                                                                       gtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCAt                                                                                                                                                                                                                                    >  7:248033/1‑149 (MQ=255)
                                                                                                                                                                             ttGGGCGTTGCAGTAGGTATGTTTTTGATGGGGCCGTGGCTGAAACGCTTGATTGCGACGCCGGAAAGCCATAAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCAt                                                                                                <  4:178724/274‑1 (MQ=255)
                                                                                                                                                                                   gTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGc                                                                                                                                                                                                                        >  7:183866/1‑149 (MQ=255)
                                                                                                                                                                                     tGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGa                                                                                                                                                                                                                      <  7:92493/149‑1 (MQ=255)
                                                                                                                                                                                       cAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATg                                                                                                                                                                                                                    >  6:112633/1‑149 (MQ=255)
                                                                                                                                                                                        agtagtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGc                                                                                                                                                                                                                   <  7:456462/149‑1 (MQ=255)
                                                                                                                                                                                                 gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGa                                                                                                                                                                                                                                 <  6:85852/126‑1 (MQ=255)
                                                                                                                                                                                                 gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGa                                                                                                                                                                                                                                 >  5:85852/1‑126 (MQ=255)
                                                                                                                                                                                                 gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTccc                                                                            >  2:209592/1‑275 (MQ=255)
                                                                                                                                                                                                    gTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagttt                                                                                                                                                                                                                                                       >  3:149300/1‑101 (MQ=255)
                                                                                                                                                                                                    gTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaaa                                                                                                                                                                                                                          >  6:134877/1‑130 (MQ=255)
                                                                                                                                                                                                    gTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaaa                                                                                                                                                                                                                          <  5:134877/130‑1 (MQ=255)
                                                                                                                                                                                                     ttGATGGTGGCGGGGCTGAAACGCATTATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTa                                                                                                                                                                                                      <  7:121110/149‑1 (MQ=255)
                                                                                                                                                                                                     ttGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGa                                                                                                                                                                                                                      <  7:333002/133‑1 (MQ=255)
                                                                                                                                                                                                     ttGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGa                                                                                                                                                                                                                      >  8:333002/1‑133 (MQ=255)
                                                                                                                                                                                                      tGATGGTGCCGTGGCTGAAACGCATGATTCCAACCCCGGAAGGCCATAAATAATTCTTGCGGAAGGAGGGGAAGAGTTATCGCTACCCTGTTTTATTTTTACACTGGCGTCCTTCAATAGATATCAAACGGAGCCCCTTTATTGACTAc                                                                                                                                                                                                     >  8:121110/1‑149 (MQ=255)
                                                                                                                                                                                                       tatgggtgcgaggCGGAAACGCAGGCTTGTGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGTTGGCGTTAGCGGTATCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACGGCGTTTGGAGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAg                                                                       <  4:23277/263‑1 (MQ=255)
                                                                                                                                                                                                        aTGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                   <  1:149300/101‑1 (MQ=255)
                                                                                                                                                                                                             gCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGttt                                                                                                                                                                                              <  7:100092/149‑1 (MQ=255)
                                                                                                                                                                                                              ccGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCATGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTTCCGCGTTtg                                                                                                                                                                                             <  8:183866/149‑1 (MQ=255)
                                                                                                                                                                                                               cGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaa                                                                                                                                                                                                                           >  2:145930/1‑118 (MQ=255)
                                                                                                                                                                                                                 tGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCcca                                                                           >  2:23277/1‑260 (MQ=255)
                                                                                                                                                                                                                 tGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGa                                                            >  2:184443/1‑275 (MQ=255)
                                                                                                                                                                                                                 gggCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGa                                                            <  2:391401/274‑1 (MQ=255)
                                                                                                                                                                                                                   gCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCg                                                                                                                                                                                                                       <  4:145929/118‑1 (MQ=255)
                                                                                                                                                                                                                                        gCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgc                                                                                                                                                                   <  8:163707/149‑1 (MQ=255)
                                                                                                                                                                                                                                           tgAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACGGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAGTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAgcgcg                                  <  1:97256/274‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                               >  8:147405/1‑128 (MQ=255)
                                                                                                                                                                                                                                                 ccATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                               <  7:147405/128‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTc                            <  2:292741/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  cATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCg                                                                                                                                                         <  5:112633/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        ttattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCt                                                                                                                                                   <  8:248033/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                   >  4:532943/1‑50 (MQ=255)
                                                                                                                                                                                                                                                           ttCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGa                                                                                                                                                                                                         >  4:548712/1‑92 (MQ=255)
                                                                                                                                                                                                                                                            tCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGt                                                                                                                                     >  4:532038/1‑159 (MQ=255)
                                                                                                                                                                                                                                                              ttGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTacca                                                                                                                                             >  5:438645/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               tGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTgg                                                                                                                                                                                                                                               <  2:532948/50‑1 (MQ=255)
                                                                                                                                                                                                                                                               tGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTAc                                                                                                                                                                                                     <  2:548717/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                gCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGGGCGTAGTCCGGGCGCTTGACCACCACGCGTTTGGTCGACAGCAAg                                                                                                                      >  4:454751/1‑170 (MQ=255)
                                                                                                                                                                                                                                                                gCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTgg                                                                                                                                 <  2:532043/159‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCAcgc                                                                                                                                       >  8:346295/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTg                                                                                                                  <  2:454756/170‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                                                                                                                                               <  5:381641/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaGCAGGGAAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                                                                                                                                               >  6:381641/1‑140 (MQ=255)
                                                                                                                                                                                                                                                                      aGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtgg                                                                                                                                                                                       >  6:462770/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                      aGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtgg                                                                                                                                                                                       <  5:462770/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCCCTAGCGCGCTTTTCTGCTTATGCGGGAACATCgg        >  1:97849/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                       gCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGt      <  2:168786/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       acAGGGTAGCGATAACGCTACCCTGCTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgc                                                                                                                                                                   <  5:406067/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                       acAGGGTAGCGATAACGCTACCCTGCTTTAGTTTTACACTGGCGTACATGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgc                                                                                                                                                                   >  6:406067/2‑118 (MQ=39)
                                                                                                                                                                                                                                                                         aGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTg                                                                                                                                  <  7:464007/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ggTAGCGTTATCGCTACCCTGTTTTAGTTTCACACTGGAGTACCTGCATAGATATCAACGCGATGCCCTTTAGTGCCTACCGCGTTTGGCGTGGCACCATTCCCCAGCGGTGGCGCGTAGTCCGGGCGCTTTAACACCACGCGTTTGGTCGCCAGCCAGCGGGCAGGGTGCAGTAATCCATCGGCAACAAGATCCGGTCCCACAAGCAACTGAAAGCAACCAATCTCTTTTTTCCCTAGCGCGCGTTTCTGCTTCTGGGGGGACAGCGGGtatgg  >  4:420377/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                             tAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGttt                                                                                                                                   <  8:158269/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                               >  6:233815/1‑100 (MQ=255)
                                                                                                                                                                                                                                                                             tAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                               <  5:233815/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                           
TGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAG  >  NC_000913/3642084‑3642625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: