Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 847,214 T→A 5.4% intergenic (‑95/+44) glnP ← / ← glnH glutamine transporter subunit/glutamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913847,2140TA 5.4% 367.5 / 7.6 103intergenic (‑95/+44)glnP/glnHglutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (51/46);  new base A (2/4);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

Only 100 of 109 total aligned reads displayed.

CAATAAAGACCAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGC  >  NC_000913/846942‑847488
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
cAATAAAGACCAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                                                                                                                                                   >  2:500669/1‑274 (MQ=255)
  aTAAAGATCAGCGCTACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATACACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATGCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgg                                                                                                                                                                                                                                                                                 <  2:158207/274‑1 (MQ=255)
        aCCAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgag                                                                                                                                                                                                                                                                             >  2:363381/1‑272 (MQ=255)
         ccAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAggg                                                                                                                                                                                                                                                                         <  4:363380/275‑1 (MQ=255)
                    gattcgtcgtgccACCCGCCGAAGGTGCAGGCAAACCCTGACAGCAATCCGCTTACCAGGGCGCCTTGCAGGCGGCGGACGGAAATCCACAGGGTCATTTTGGGACCGTCAAACAGAAGCGGAATGGCAGGACAGATGGCACTCCATGCACACTGCATCTGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAGTTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCGCTCCCCATAGGGGAGAGGGAACTGCCAGTg                                                                                                                                                                                                                                                              <  4:459728/264‑1 (MQ=255)
                       ttGGCTATCCAACCTCCGAAGGGGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCCATAGGGGAGAGGGAACTGCCAGTg                                                                                                                                                                                                                                                              <  2:126960/273‑1 (MQ=255)
                              tCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCCGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCCTATGTTGTTCCTGTTACCGTGGTGTAACTAATCCAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCCACCCATAGGGGAGAGGGAACTGCCATTTCGGTTTTCAg                                                                                                                                                                                                                                                    >  4:532014/1‑275 (MQ=255)
                                               cttgaaaaaCCGGCCAGCAAGCCGAGTACCAAACCGCCTGCCAGACCGAGGACCGACATCCACAGGGTCACTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATAGGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTGGCGCCCCTCACCATAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAttt                                                                                                                                                                                                                                   <  4:50951/270‑1 (MQ=255)
                                               cgtgcAAAACCTGCCAGCAATACGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGACGCGGAATGGCAGGCCAGATGTCCCTACAGTCAAACTGCATATGTTGTTCCTGGTACCGTGGTGTAAATAATCAAATTAATGGCAATTAATTGCCCCGGCGCGGGTCTGCGCACCTCACCCTAAGCCACTCACCATCGGGGAAAGGGAACTTGCGGTGCGCTTTACCGGTGTTGCGTCATTAttt                                                                                                                                                                                                                                   >  4:182302/1‑275 (MQ=255)
                                                                                   ccTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGCGGGAACTGCCAGTGCGTTTTACCGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCCTTTTTTGTAGTTTTCGt                                                                                                                                                                                                >  2:412616/1‑274 (MQ=255)
                                                                                    cTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                                                                                                                                               >  2:162497/1‑274 (MQ=255)
                                                                                                    cGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCCTTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTcg                                                                                                                                                                               >  2:450468/1‑274 (MQ=255)
                                                                                                        aTCCACAGGGTCATTTTGGACCCTCAAATCAGAAGCGGAATGGATGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGGGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAg                                                                                                                                                                          <  4:352718/275‑1 (MQ=255)
                                                                                                         gccACAGGGTAATTTTGGAACGTTAAATCGGAAGCGGAATGGCAGGCCAGATGGCACTCCAGGAAAACTGCATAGGTTGTTCCTGTGACCGGGGTGTAAATAATCAAAGGCAGGGTAATTAATGGCACCGGAGCGGGTTTGCGCCCCTCACCCGAAGCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTATCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAg                                                                                                                                                                          <  4:212295/273‑1 (MQ=255)
                                                                                                          ccACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTAc                                                                                                                                                                                                                                                      <  4:319532/197‑1 (MQ=255)
                                                                                                                             ccTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAt                                                                                                                                                                                                                                                                                    >  8:16018/1‑148 (MQ=255)
                                                                                                                              cTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                                                                                                                                                    >  2:311453/1‑275 (MQ=255)
                                                                                                                                    tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgag                                                                                                                                                                                                                                                                             <  7:16018/148‑1 (MQ=255)
                                                                                                                                           cGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCt                                                                                                                                       <  4:162496/275‑1 (MQ=255)
                                                                                                                                             gAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGAGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCTTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCAAGCCGCtct                                                                                                                                     >  2:210192/1‑274 (MQ=255)
                                                                                                                                              aaTGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGgag                                                                                                                                                                                                                                                                             >  7:43344/1‑138 (MQ=255)
                                                                                                                                                   cAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgggg                                                                                                                                                                                                                                                                               <  7:215505/131‑1 (MQ=255)
                                                                                                                                                   cAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgggg                                                                                                                                                                                                                                                                               >  8:215505/1‑131 (MQ=255)
                                                                                                                                                                   ccAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGgaga                                                                                                                                                                                                                                                                            <  5:681332/118‑1 (MQ=255)
                                                                                                                                                                   ccAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAg                                                                                                                                                                                                                                              >  8:333957/1‑148 (MQ=255)
                                                                                                                                                                    cacgcANACTTCATATGTTTTTCCTGTTACCGTGGTGTAAAATATACACATCGGGGGACTTAAGGGCCCCGGAGCGGCTGTGTGCCCCGTGCCCTCATCCCCCCCCCATAGGGGGGAGGGAACGGCCAGGGCGTTTTACTGGTGTAGCGTTGTTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCaa                                                                                                               <  4:511808/270‑1 (MQ=255)
                                                                                                                                                                    cAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGc                                                                                                                                                                                                                                             <  7:333957/148‑1 (MQ=255)
                                                                                                                                                                    cAGTCAAACTGCAAATGTTTTTCCTGTTACCGGGGTGTAAAAAATCAAATTAAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTAACCCTAATCCTCACCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCaa                                                                                                               <  4:450467/274‑1 (MQ=255)
                                                                                                                                                                           aCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                                                                                                                                                   <  6:636236/103‑1 (MQ=255)
                                                                                                                                                                           aCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                                                                                                                                                   >  5:636236/1‑103 (MQ=255)
                                                                                                                                                                              gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTTCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGTTTTCGTTGTAATTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCATCTCGTCGCTCCCTTTCGGGAACGCAAACCCGTTTtgc                                                                                                    >  4:123020/1‑275 (MQ=255)
                                                                                                                                                                                         tCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGa                                                                                                                                                                                                                                                                        >  7:645189/1‑100 (MQ=255)
                                                                                                                                                                                         tCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGa                                                                                                                                                                                                                                                                        <  8:645189/100‑1 (MQ=255)
                                                                                                                                                                                         tCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTg                                                                                                                                                                                                                                                                    <  7:579092/104‑1 (MQ=255)
                                                                                                                                                                                         tCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTg                                                                                                                                                                                                                                                                    >  8:579092/1‑104 (MQ=255)
                                                                                                                                                                                          ccTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCTGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCag                                                                                         >  2:93725/1‑274 (MQ=255)
                                                                                                                                                                                           cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                                                                                                                                                   <  5:1117705/87‑1 (MQ=255)
                                                                                                                                                                                           cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                                                                                                                                                   >  6:1117705/1‑87 (MQ=255)
                                                                                                                                                                                            tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGtta                                                                                                                                                                                                                                         <  7:200583/128‑1 (MQ=255)
                                                                                                                                                                                            tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGtta                                                                                                                                                                                                                                         >  8:200583/1‑128 (MQ=255)
                                                                                                                                                                                              ttACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATCCCGTATTGCTGCGCTTCCAGAGAGt                                                                                    >  2:439469/1‑275 (MQ=255)
                                                                                                                                                                                              ttACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCAACCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGtttt                                                                                                                                                                                                                                                        <  6:323757/111‑1 (MQ=255)
                                                                                                                                                                                              ttACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCAACCTAATCCTCTCCCCATAGGGGAGAGCGAACTGCCAGTGCGtttt                                                                                                                                                                                                                                                        >  5:323757/1‑111 (MQ=255)
                                                                                                                                                                                               tACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGaa                                                                                                                                                                                                                                                                       >  6:969510/1‑95 (MQ=255)
                                                                                                                                                                                               tACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGCaa                                                                                                                                                                                                                                                                       <  5:969510/95‑1 (MQ=255)
                                                                                                                                                                                                        gtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTTCGCTTCCAGAGAGTCACCTCCCGc                                                                          >  2:302678/1‑275 (MQ=255)
                                                                                                                                                                                                         tAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCCCCCCTCACCCTAATCCGCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCt                                                                         <  1:44488/275‑1 (MQ=255)
                                                                                                                                                                                                           aataatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAt                                                                                                                                                                                                                                                                                    >  8:656798/1‑70 (MQ=255)
                                                                                                                                                                                                           aataatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAt                                                                                                                                                                                                                                                                                    <  7:656798/70‑1 (MQ=255)
                                                                                                                                                                                                               atCAAATTCTGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGt                                                                                                                                                                                                                               >  7:68792/1‑119 (MQ=255)
                                                                                                                                                                                                               atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGt                                                                                                                                                                                                                               <  8:68792/119‑1 (MQ=255)
                                                                                                                                                                                                                 cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTAc                                                                                                                                                                                                                                                      >  6:900162/1‑94 (MQ=255)
                                                                                                                                                                                                                 cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACa                                                                                                                                                                                                                                                     >  7:173472/1‑95 (MQ=255)
                                                                                                                                                                                                                            taattaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAg                                                                                                                                                                                                                                                                                  >  5:671658/1‑55 (MQ=255)
                                                                                                                                                                                                                            taattaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAg                                                                                                                                                                                                                                                                                  <  6:671658/55‑1 (MQ=255)
                                                                                                                                                                                                                             aattaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACa                                                                                                                                                                                                                                                     <  5:342514/83‑1 (MQ=255)
                                                                                                                                                                                                                             aattaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACa                                                                                                                                                                                                                                                     >  6:342514/1‑83 (MQ=255)
                                                                                                                                                                                                                                      ccccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTTCGCTTCCAGAAAGTCACCTACCGCTTTGAACTGACGGTTACCGGCGGTTTtgtt                                            >  2:475430/1‑273 (MQ=255)
                                                                                                                                                                                                                                        ccGGGGCGGGGTTGGGCCCCTCACCCAAACCCGCTCCCCAAAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTtgatg                                           <  4:302677/274‑1 (MQ=255)
                                                                                                                                                                                                                                        ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCGATCCCGTATTGCTGCGCTTCCAGGGAGTCACCTACCGATTTGAACTGACCGTTCCCGGCGGTTTtgagga                                          >  4:184726/1‑275 (MQ=255)
                                                                                                                                                                                                                                         ctgggCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTtgatga                                          <  4:93724/270‑1 (MQ=255)
                                                                                                                                                                                                                                         cGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATATGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTtgatg                                           <  2:134109/273‑1 (MQ=255)
                                                                                                                                                                                                                                          ggAGCGGGTTTGCGCCTCTCACCCTAATCCTCTCCCCCTAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCCCGCAGCTCGTCGCTACCTTTCGGGAACGCCATACCGTATTGCTGCGCTTCCCTAGAGCCCCCTCCCCCTTTGAACTTACCGTTACCGGCGGTTTGTATGaaa                                        >  4:255593/1‑274 (MQ=255)
                                                                                                                                                                                                                                                 gTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAgg                                                                                                                                                                                                                                                                          >  8:1092105/1‑42 (MQ=255)
                                                                                                                                                                                                                                                 gTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAgg                                                                                                                                                                                                                                                                          <  7:1092105/42‑1 (MQ=255)
                                                                                                                                                                                                                                                  tttGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGt                                                                                                                                                                                                                                           <  5:185509/72‑1 (MQ=255)
                                                                                                                                                                                                                                                  tttGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGt                                                                                                                                                                                                                                           >  6:185509/1‑72 (MQ=255)
                                                                                                                                                                                                                                                    tGCGCCCCTCACGCTAATCCGCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGTGTCAGTATCGAACCAGTTTTTGTAGATTTCGTTTTAAGTTCCGTTCTCGCGCACGCTTATAAACGCGCAGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCATCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGt                              <  4:255710/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAttttt                                                                                                                                                                                                            <  6:34382/96‑1 (MQ=255)
                                                                                                                                                                                                                                                         cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAttttt                                                                                                                                                                                                            >  5:34382/1‑96 (MQ=255)
                                                                                                                                                                                                                                                          ccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGc                         <  2:215859/274‑1 (MQ=255)
                                                                                                                                                                                                                                                          ccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGc                         <  2:532015/274‑1 (MQ=255)
                                                                                                                                                                                                                                                               cccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGgtgt                                                                                                                                                                                                                                                <  8:660286/54‑1 (MQ=255)
                                                                                                                                                                                                                                                               cccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGgtgt                                                                                                                                                                                                                                                >  7:660286/1‑54 (MQ=255)
                                                                                                                                                                                                                                                                  tAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTTTGATTTCGGTTCAGTACCGATCCATTTGTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGGTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGc                <  2:182303/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ctctCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGtt                                                                                                                                                                                   <  6:329133/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tctcCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGAc                                                                                                                                                     >  5:719499/1‑136 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttt                                                                                                                                                                                                             >  4:555640/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCCCGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCCGATAGTCACCTCCCGCTTTGAACTGACCGTTACCGGCGGTTTTGTTGAAGTACCGCATGCTTTGCGTATCGTGCAGAACGGc        >  1:49542/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATtgc                                                                                                    >  3:137255/1‑183 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCCCCTACCGCTTTGAACTGCCCGTTACCGGCGGTTTTGATGAAGTCCAGAATGTTTGGCGTATCGTGCAGAACGGc        >  4:313426/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTTACCTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGc        >  4:478735/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGTTGAATTACAGAATGTTTGGCGTATCGTGCAGAACGGc        >  2:362931/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGc        >  2:322385/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                          tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGNCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGANTGTTACCGGGGGTTTTGATGAAGTANAGAAGTTTTGGCGTTTTGTGCAGAACggg        >  3:94214/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                            ccccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGGCGCGTTGaaa                                                                                                                                                    >  7:178818/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                             cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtct                                                                                                                                                                                 <  7:454440/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtct                                                                                                                                                                                 >  8:454440/1‑104 (MQ=255)
                                                                                                                                                                                                                                                                               cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGt                                                                                                                                                                                                         <  2:555641/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAg                                                                                                                                                                                                       >  6:925039/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                               cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAg                                                                                                                                                                                                       <  5:925039/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAg                                                                                                                                                                                                       <  8:340657/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                                                                                                                                                      >  8:1008831/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                               cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                                                                                                                                                      <  7:1008831/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTgcg                                                                                                <  1:137255/183‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aTAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGc                                                                                                                                           <  8:625707/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aTAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGc                                                                                                                                           >  7:625707/1‑139 (MQ=255)
                                                                                                                                                                                                                                                                                 tAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg                                                                                                                                  >  5:240120/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                                 tAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg                                                                                                                                  <  6:240120/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTTCAGAACGGCGTTTgg  >  2:63057/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
CAATAAAGACCAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGC  >  NC_000913/846942‑847488

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