Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,882 T→G 5.3% I163S (ATC→AGC)  yagG → CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913285,8820TG 5.3% 920.8 / 37.3 269I163S (ATC→AGC) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base T (113/140);  new base G (8/6);  total (123/146)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-02

Only 100 of 270 total aligned reads displayed.

GCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGT  >  NC_000913/285613‑286155
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 
   cGGCCGGGTCCGCCCTGTCCTGGTGTTGTGGTCATGTCCGTTGAGCACGGTCGGGGTGCGAACGTTATACCCGCCGGACTGCTCGGCGCGTGGCAAGTTCGTCTACGCCGGCGTAAACTACATTCGCACGACCCTGGTCTGCACGTGCGTTACCGTGCCGTACTGCGCCATGCCGGGCGTCCTCGCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTg                                                                                                                                                                                                                                                                            <  4:332671/274‑1 (MQ=255)
               cccGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGTGTCATCACCGCCGACCCGAAAGAGCGTCACTCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGATcctca                                                                                                                                                                                                                                                            >  1:73848/1‑278 (MQ=255)
                    tCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCCCCCGCTGGTGAGCATCATCggaa                                                                                                                                                                                                                                                      >  1:112688/1‑277 (MQ=255)
                                                                                aCTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCNGTCCTGCCGCCCCTCCCTGGCCGAGGCGGGCTCGCTCGCCCTCCGCGGCACCGCGCTGCCGCGGGTGAGCAGCACCCGCANTGCGGAGGACCAGGTGTGCCACATTNGCGGCACGTGCGTCCTcggggggccgggcggg                                                                                                                                                                                           >  3:182970/1‑266 (MQ=255)
                                                                                              ggCAAGATCATCTACGCCTGCGTGCCCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCCCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGCCTCCCTTTGCCTGCTGGGGCTGAGCGACGTGTTTCTGctcct                                                                                                                                                                                >  2:126535/1‑273 (MQ=255)
                                                                                              ggCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGCCTACACCTTCGTTACCGTGCCGTACTGCGCCATGCCGGGCGTCCTCACCGCCGACCCGAAAGAGCGTCACGCCCTCNCCTCCTGGCGCCTCTTCCCGGCCNCGGCNGGCTCGCTCGCCATCAGCGGCATCCCCCTGCCGCGCGTGCGCCCCATCCGCaa                                                                                                                                                                                                                                                     >  4:696294/1‑206 (MQ=255)
                                                                                                  ggATCATCTACGCCTGCGTGACCTACATTCTGCGGACGCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcggcg                                                                                                                                                                                             <  2:696299/257‑1 (MQ=255)
                                                                                                             gCCTGCGTGACCTACATTCTCCTGACCCGGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCa                                                                                                                                                                  <  2:390766/274‑1 (MQ=255)
                                                                                                               cggCGTGACCTACCTTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcggcg                                                                                                                                                                                             <  1:27059/243‑1 (MQ=255)
                                                                                                                  cGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCg                                                                                                                                                                                                                                                         >  4:229132/1‑182 (MQ=255)
                                                                                                                  cGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCTCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCTCGAcc                                                                                                                                                             >  2:176994/1‑274 (MQ=255)
                                                                                                                      aCCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaa                                                                                                                                                                                                                                                     <  2:229135/182‑1 (MQ=255)
                                                                                                                      aCCTACATTCTCCGGACCCTGGGTGAACCCTTCGTTAACGGGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaaa                                                                                                                                                          <  4:422902/273‑1 (MQ=255)
                                                                                                                         tACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATc                                                                                                                                                                                                                                                                                  >  6:953258/1‑149 (MQ=255)
                                                                                                                          aCATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATc                                                                                                                                                                                                                                                                                  >  8:371337/1‑149 (MQ=255)
                                                                                                                            aTTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATcgc                                                                                                                                                                                                                                                                                >  7:1642656/1‑149 (MQ=255)
                                                                                                                            aTTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTa                                                                                                                                                  <  2:481886/275‑1 (MQ=255)
                                                                                                                             ttCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGTGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATcgcg                                                                                                                                                                                                                                                                               >  8:1213049/1‑149 (MQ=255)
                                                                                                                             ttCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATcgcg                                                                                                                                                                                                                                                                               >  6:1551852/1‑149 (MQ=255)
                                                                                                                               ctcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCt                                                                                                                                                                                                                                                                             <  6:1398431/149‑1 (MQ=255)
                                                                                                                                     ccccTGCTCTACACCTTCGTTAACGTGCCGTACGGCNCCATGCCGNGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACg                                                                                                                                                                                                                                             <  3:51272/174‑1 (MQ=255)
                                                                                                                                      cccTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTg                                                                                                                                                                                                                                                                      >  8:913614/1‑149 (MQ=255)
                                                                                                                                       ccTGATCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCGCTga                                                                                                                                                                                                                                                                     <  6:1020252/149‑2 (MQ=255)
                                                                                                                                         tGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTg                                                                                                                                                                                                                                                                   >  8:167017/1‑149 (MQ=255)
                                                                                                                                         tGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCa                                                                                                                                                            <  2:617837/252‑1 (MQ=255)
                                                                                                                                         tGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCGCTgata                                                                                                                                                                                                                                                                   <  5:1519330/149‑4 (MQ=255)
                                                                                                                                          ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                                                                                                                                                  <  5:453129/149‑1 (MQ=255)
                                                                                                                                          ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                                                                                                                                                  <  5:1258950/149‑1 (MQ=255)
                                                                                                                                           gTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATc                                                                                                                                                                                                                                                                                  >  5:1176488/1‑132 (MQ=255)
                                                                                                                                            tCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGc                                                                                                                                                                                                                                                                <  7:1562032/149‑1 (MQ=255)
                                                                                                                                              tACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGcat                                                                                                                                                                                                                                                              >  5:699882/1‑149 (MQ=255)
                                                                                                                                                  ccTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCAACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCatcatc                                                                                                                                                                                                                                                          <  7:301401/149‑1 (MQ=255)
                                                                                                                                                   cTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCg                                                                                                                                                                                                                                                         >  5:970077/1‑149 (MQ=255)
                                                                                                                                                    ttCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCgg                                                                                                                                                                                                                                                        <  7:1072277/149‑1 (MQ=255)
                                                                                                                                                    ttCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaa                                                                                                                                                                                                                                                     <  4:460424/152‑1 (MQ=255)
                                                                                                                                                        ttAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaaa                                                                                                                                                                                                                                                    <  6:707828/149‑1 (MQ=255)
                                                                                                                                                          aaCGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAgg                                                                                                                                                                                                                                                  <  8:1007415/149‑1 (MQ=255)
                                                                                                                                                           ctttGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAggg                                                                                                                                                                                                                                                 <  7:926443/146‑1 (MQ=255)
                                                                                                                                                                ccGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGa                                                                                                                                                                                                                                            <  7:316572/149‑1 (MQ=255)
                                                                                                                                                                ccGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGa                                                                                                                                                                                                                                            <  6:200037/149‑1 (MQ=255)
                                                                                                                                                                ccGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgc                                                                                                                                                                                                 >  2:599059/1‑192 (MQ=255)
                                                                                                                                                                 cGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGCGCAGGTGGGCTACTTCGGCGCCATGTGCGTCCTGGGGCTGAGCGGCGTGGTGCTTCTCTACCGCTGCttcttcccc                                                                                                                                                                >  4:177694/1‑221 (MQ=255)
                                                                                                                                                                     cTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCa                                                                                                                                                            <  2:177697/224‑1 (MQ=255)
                                                                                                                                                                       gcgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTg                                                                                                                                                                                                                                     <  8:1476022/149‑1 (MQ=255)
                                                                                                                                                                        cgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTgg                                                                                                                                                                                                                                    <  7:1421682/149‑1 (MQ=255)
                                                                                                                                                                        cgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTa                                                                                                      <  2:477133/275‑1 (MQ=255)
                                                                                                                                                                          ccATGCCGGGCGTCATCACCGCCGACCCGAAGGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTCCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                                                                                                                                                  <  8:1358778/149‑1 (MQ=255)
                                                                                                                                                                          ccATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGcatcat                                                                                                                                                                                                                                                           >  5:253376/1‑124 (MQ=255)
                                                                                                                                                                          ccATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                                                                                                                                                  <  8:438132/149‑1 (MQ=255)
                                                                                                                                                                          ccATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCa                                                                                                                                                                                                                     >  2:413370/1‑162 (MQ=255)
                                                                                                                                                                            aTGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcggcg                                                                                                                                                                                             <  4:598082/184‑1 (MQ=255)
                                                                                                                                                                             tGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCCCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAgct                                                                                                  >  2:455140/1‑274 (MQ=255)
                                                                                                                                                                              gCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgtg                                                                                                                                                                                                                 <  4:413366/162‑1 (MQ=255)
                                                                                                                                                                                cGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTc                                                                                                                                                                                                                            <  6:995644/149‑1 (MQ=255)
                                                                                                                                                                                   gCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGgc                                                                                                                                                                                                                         <  7:757084/149‑1 (MQ=255)
                                                                                                                                                                                     gTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCCTCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGgcgc                                                                                                                                                                                                                       >  8:692182/1‑149 (MQ=255)
                                                                                                                                                                                      tcatcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCGTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGcc                                                                                                                                                                                                                      >  6:340303/1‑149 (MQ=255)
                                                                                                                                                                                      tcatcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGCCGGCGGGCTCGATCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGcc                                                                                                                                                                                                                      <  6:1225379/149‑1 (MQ=255)
                                                                                                                                                                                       catcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCCGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCa                                                                                                                                                                                                                     <  8:78433/149‑1 (MQ=255)
                                                                                                                                                                                            ccgccgACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgc                                                                                                                                                                                                                <  7:607601/149‑1 (MQ=255)
                                                                                                                                                                                            ccgccgACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgc                                                                                                                                                                                                                >  8:1635616/1‑149 (MQ=255)
                                                                                                                                                                                             cgccgACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcg                                                                                                                                                                                                               >  6:573970/1‑149 (MQ=255)
                                                                                                                                                                                               ccgACCCGAAAGAGCGTCACGCCCTTCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcgtg                                                                                                                                                                                                             <  8:195973/149‑1 (MQ=255)
                                                                                                                                                                                                 gACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCt                                                                                                                                                                                                           >  5:104951/1‑149 (MQ=255)
                                                                                                                                                                                                   cccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgg                                                                                                                                                                                                         <  5:1352803/149‑1 (MQ=255)
                                                                                                                                                                                                      gAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGc                                                                                                                                                                                                      >  7:552734/1‑149 (MQ=255)
                                                                                                                                                                                                       aaaGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCt                                                                                                                                                                                                     <  8:468984/149‑1 (MQ=255)
                                                                                                                                                                                                             cGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCAtgtg                                                                 <  2:546321/275‑1 (MQ=255)
                                                                                                                                                                                                               tCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcggcg                                                                                                                                                                                             >  8:601959/1‑149 (MQ=255)
                                                                                                                                                                                                                 aCGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGc                                                                                                                                                                                                                                          <  7:827638/102‑1 (MQ=255)
                                                                                                                                                                                                                    cccTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggtg                                                                                                                                                                                        <  5:340303/149‑1 (MQ=255)
                                                                                                                                                                                                                        gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctgc                                                                                                                                                                                    <  5:1544563/149‑1 (MQ=255)
                                                                                                                                                                                                                             ccTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGGCCTTCAGCTTCTCCTGGGCCACAGCCAGTGCCGCAGAATTGGCGAGTTCAAGatgctg                                                 >  4:264140/1‑275 (MQ=255)
                                                                                                                                                                                                                              cTGGCGCATCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTggg                                                                                                                                                                                                                                   <  8:1161744/96‑1 (MQ=255)
                                                                                                                                                                                                                               tGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGt                                                                                                                                                                             <  7:966180/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cgcgTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGAc                                            <  2:264143/275‑1 (MQ=255)
                                                                                                                                                                                                                                      tcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGcttc                                                                                                                                                                      <  6:77280/149‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCa                                                                                                                                                                  <  7:335061/149‑1 (MQ=255)
                                                                                                                                                                                                                                             ggcggcggcGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGa                                                                                                                                                               <  6:36063/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     cggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAagag                                                                                                                                                       <  7:427775/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         ctcgctcgCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgg                                                                                                                                                                                                         >  6:1146104/1‑95 (MQ=255)
                                                                                                                                                                                                                                                         ctcgctcgCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgg                                                                                                                                                                                                         <  5:1146104/95‑1 (MQ=255)
                                                                                                                                                                                                                                                          tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAagag                                                                                                                                                       <  6:1653829/144‑1 (MQ=255)
                                                                                                                                                                                                                                                            gctcgctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTaca                                                                                                                                                <  6:481307/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             ctcgctATCAGCGGCATCGCGCTGCCGCTGGTGCGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCgggg                  <  2:7811/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                gctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCtt                                                                                                                                            >  8:957420/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                gcgaTCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCtt                                                                                                                                            >  6:1177202/4‑149 (MQ=255)
                                                                                                                                                                                                                                                                 ctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCCTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCttt                                                                                                                                           <  7:1175312/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTg                                                                                                                                          >  6:1297054/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aTCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGa                                                                                                                                         <  6:801770/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg                                                                                                                                        <  7:167017/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg                                                                                                                                        <  5:573970/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaaa                                                                                                                                                          >  5:1699667/2‑131 (MQ=255)
                                                                                                                                                                                                                                                                     cAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg                                                                                                                                       >  7:675525/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg                                                                                                                                       >  5:757636/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           cATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                                                                                                                                                  <  1:142478/52‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGcc                                                                                                                                >  7:1213613/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGcc                                                                                                                                >  5:374760/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGGCGCGCTACTTCGt  >  2:109249/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                             tCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCg                                                                                                                               >  8:306709/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 
GCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGT  >  NC_000913/285613‑286155

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