Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,642,361 T→A 5.8% intergenic (+30/+19) 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,3610TA 5.8% 358.4 / 7.8 104intergenic (+30/+19)dtpB/rsmJdipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (46/51);  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 = 8.65e-01

Only 100 of 106 total aligned reads displayed.

TTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACC  >  NC_000913/3642094‑3642635
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                  
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)
                                                          ggcggcTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGa                                                                                                                                                                                                                   >  2:275937/1‑275 (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)
                                                                                                                                 taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTGAATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTg                                                                                                                                                                                                                                                                          >  8:464007/1‑149 (MQ=255)
                                                                                                                                   ccAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtt                                                                                                                                                                                                                                                                        <  5:286052/149‑1 (MQ=255)
                                                                                                                                     aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                                                                                                                                                      >  5:78121/1‑149 (MQ=255)
                                                                                                                                     aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                                                                                                                                                      >  8:456462/1‑149 (MQ=255)
                                                                                                                                      aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgc                                                                                                                                                                                                                                                                         <  8:282330/145‑2 (MQ=255)
                                                                                                                                      aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgc                                                                                                                                                                                                                                                                         >  7:282330/1‑144 (MQ=255)
                                                                                                                                       cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGTCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGc                                                                                                                                       >  2:252677/1‑274 (MQ=255)
                                                                                                                                        gtgtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagt                                                                                                                                                                                                                                                                   <  5:377468/149‑1 (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                                                                                                                                                                                                                                           >  5:85852/1‑126 (MQ=255)
                                                                                                                                                                                       gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGa                                                                                                                                                                                                                                           <  6:85852/126‑1 (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                                                                                                                                                                                                                                >  8:333002/1‑133 (MQ=255)
                                                                                                                                                                                           ttGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGa                                                                                                                                                                                                                                <  7:333002/133‑1 (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)
                                                                                                                                                                                                                                            aattaTTCTTGCAGAAACGTGGTTACCGTTATCGCTACCCTGTTTTAGTTTTTCACTGGCTTACCTGAATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTGGGCGTGGCGACATTCGCCAGCAGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAAGCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCt                                  <  4:184441/274‑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)
                                                                                                                                                                                                                                                      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)
                                                                                                                                                                                                                                                                   tAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGttt                                                                                                                                             <  8:158269/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGttt                                                                                                                                             >  7:158269/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                   tAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                                         >  6:233815/1‑100 (MQ=255)
                                                                                                                                                                                                                                                                   tAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                                         <  5:233815/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tCGCTACCCTGTTTTAGTTTTACACTGGCATACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAaccacc  >  2:158751/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                  
TTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACC  >  NC_000913/3642094‑3642635

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