Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,887 C→A 14.8% L165M (CTG→ATG)  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,8870CA14.8% 255.4 / 33.4 88L165M (CTG→ATG) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base A (6/7);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-02

GCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAAT  >  NC_000913/285618‑286160
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                    
gCCGGTTCGGCCCGTTCCTGCTGTGGGGGGCGATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACCCGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACGTTCTCCTGACCCTGGTCTACACCGTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCg                                                                                                                                                                                                                                                                                 <  4:337375/275‑1 (MQ=255)
              ttCCTGCTGTGGGGGGCCATCCCGGTCGGAATCGTCTGCGGGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCAGCTCCGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGcatcat                                                                                                                                                                                                                                                                <  4:392691/275‑1 (MQ=255)
                    gtgtgGGGGGCGCTACCGTTCGGCATCTTCTGCGTGCTGACCTTCTACCCGCCGGACTTCTCCGCACAGGGCAAGTTCATCTACGCCTGCGTGACCTACATTCTCCGGACCCTGGTCTACACCTGCGTTATCGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCaa                                                                                                                                                                                                                                                          <  4:498793/274‑1 (MQ=255)
                         gggggCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCCGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAACgggg                                                                                                                                                                                                                                                     >  2:224222/1‑275 (MQ=255)
                            ggCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTCCGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCCCGCCCTGCAGTCCTGGAGCTTCTTCCTGGCGGCGGCGGGCTCGCGCGCTATCCGCGGCATCGCGCTGCCGCTG‑‑‑GTGCGCATCATAGGCCAAGGGGa                                                                                                                                                                                                                                                    >  1:91473/1‑273 (MQ=255)
                            gcccgtCCCGTTCGGCCTCGTCTGCGGGCTGCCCTTCTACACGCCGTACTTCTCCGCACAGGGCAATATCATCTACGNGTGCGTGACCTACATTCTCCTGCCCCTGGTCTACACCTTCGTTAACGGGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACg                                                                                                                                                                                                                                                  <  3:129533/270‑1 (MQ=255)
                                                     gTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGcc                                                                                                                                                                                                                           >  2:124888/1‑273 (MQ=255)
                                                         tGACCGTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGAACCTGGTCTACACCTTCGTTCACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgtg                                                                                                                                                                                                                      <  1:109409/274‑1 (MQ=255)
                                                                        cGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATAATCGTCAAAGGGGCCGAGCAGCTGGGCTACTTCGGCGCCTTGTGCGTGCTGGGGCTggc                                                                                                                                                                                                       >  1:19527/1‑272 (MQ=255)
                                                                            ctcatcCGCACAGGGCAGGTTCATCAACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTACGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATTACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcggcg                                                                                                                                                                                                  <  2:511364/271‑1 (MQ=255)
                                                                            cTTCTCCGGCCGGGGCAAGTTCATCGACGCCGGCGGGACCTACATCCTCCTGACCCTGGTCTACACCTGCGTTAACGGGCCGTACTGCGCCATGCCGGGCGTCAGCACCGCCGAGCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCGCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcggcg                                                                                                                                                                                                  <  4:492460/275‑1 (MQ=255)
                                                                               ctcCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCCCGCCCGCCATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCCTCCTCGGCAAAGGGGACGAGACGTTGGGCTCCTTCGGCGCCATGTGCGTGCTGGGGCTGCggggcggg                                                                                                                                                                                                >  4:255953/1‑272 (MQ=255)
                                                                                  cGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACAGTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgc                                                                                                                                                                                            <  2:487975/275‑1 (MQ=255)
                                                                                     acaGGGCAAGATCATCTACGGCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctgc                                                                                                                                                                                         <  2:100682/275‑1 (MQ=255)
                                                                                         ggCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGccct                                                                                                                                                                                      >  2:159371/1‑274 (MQ=255)
                                                                                             aGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgtg                                                                                                                                                                                                                      <  2:178936/238‑1 (MQ=255)
                                                                                                   tcTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGCGCATCATCGGCCAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCCTGTGCGTGCTGGGGCTGAGCGGCGTGGTTCTGCTCTACCTcgtagtc                                                                                                                                                                           >  1:12028/1‑269 (MQ=255)
                                                                                                      aCGCCTGTGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGGGCTGCTCTACGTCTGcttctt                                                                                                                                                                         >  2:44596/1‑274 (MQ=255)
                                                                                                        gCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCa                                                                                                                                                                       <  2:369401/274‑1 (MQ=255)
                                                                                                             cGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCCCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGccc                                                                                                                                                                  >  2:68772/1‑274 (MQ=255)
                                                                                                                    tACCTTCTCCTGCCCCTGGTCTACACCTTCGTTGACGTGCCGTACTGCGCCAGGCCGAGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTCAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgcg                                                                                                                                                          <  4:205676/275‑1 (MQ=255)
                                                                                                                       atttcccTGACCCTTGTTCACCCCTTCGTTACCGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTa                                                                                                                                                       <  4:44595/270‑1 (MQ=255)
                                                                                                                               gACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCg                                                                                                                                                                                                                                                                                 <  8:897245/148‑1 (MQ=255)
                                                                                                                                       tCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTTTTCCTGGCGGCGGCGGGCTCGCTCGCTATCTGCGGCATCGCGCTGCCTCTG‑‑‑GTGAGCTTCATCGGCACAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTgacg                                                                                                                                       >  2:375651/1‑272 (MQ=255)
                                                                                                                                                tCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCCCCGCCGACCCCAAAGAGCGTCCCGCCCTGCAGTCCTGGCCCTTCTTCCTTGCGCCCGCGGGCTCGCCCGCGATCACCCGCCTCGCTCCGCCCCCG‑‑‑TCGGCCATCCTCGGCAACNCTCCCGCCCAGGTGTGCCACTTCTGTCCCATGTGTGTGCTtggtgcccnncgc                                                                                                                                                                                                   >  4:544624/1‑196 (MQ=255)
                                                                                                                                                tCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCgcgc                                                                                                                                                                                                                                                                                      >  7:1075341/1‑126 (MQ=255)
                                                                                                                                                tCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCgcgc                                                                                                                                                                                                                                                                                      <  8:1075341/126‑1 (MQ=255)
                                                                                                                                                tCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTTCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCCCGACCAAAGACCGCTACACCTTTGCGGTGCACCCGGGCtc                                                                                                                               >  2:275773/1‑274 (MQ=255)
                                                                                                                                                tCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCCCTACGTCTGCTTCTTCACGACCAAAGAGCGCCAAACCTTTGGGGGGCAGCCGGGCtgg                                                                                                                              >  4:411009/1‑273 (MQ=255)
                                                                                                                                                                        tGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGCGCGGCGTGGTTCTGCTCCACGTCTGCTTCTTCACGTCAAAAGAGCGCCACCCCTTTGAGGTGCCGCCGGGCTCGCCGGTGGCGAAAGACCTGAGgct                                                                                                       >  4:430070/1‑274 (MQ=255)
                                                                                                                                                                        tGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAgct                                                                                                       >  2:344465/1‑274 (MQ=255)
                                                                                                                                                                          ccGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCCCGCTCGCGATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCCTCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCCTGTTCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAACGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGCTGGTGGCGTAAGACCTTAAgctgc                                                                                                     >  4:444136/1‑274 (MQ=255)
                                                                                                                                                                                  catcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTggg                                                                                                                                                                                                                                        >  6:354337/1‑135 (MQ=255)
                                                                                                                                                                                  catcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTggg                                                                                                                                                                                                                                        <  5:354337/135‑1 (MQ=255)
                                                                                                                                                                                   atcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCa                                                                                            >  2:417037/1‑274 (MQ=255)
                                                                                                                                                                                   acagccgccgACCCGAAAGGGCGTAACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCGCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCa                                                                                            <  2:521051/270‑1 (MQ=255)
                                                                                                                                                                                      aCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgt                                                                                                                                                                                                                       >  7:926430/1‑148 (MQ=255)
                                                                                                                                                                                                  aaaGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGCGCATCATCGGCaaa                                                                                                                                                                                                                                                         >  6:170249/1‑102 (MQ=255)
                                                                                                                                                                                                  aaaGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCaaa                                                                                                                                                                                                                                                         <  5:170249/102‑1 (MQ=255)
                                                                                                                                                                                                  aaaGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTTCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTTAGGTGCAGCCGGGCTAGTCGGTGGCGAAAGACCTTACACTGTCCAGGGGCACCAACCAGTGgct                                                                               >  4:521050/1‑271 (MQ=255)
                                                                                                                                                                                                        cGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTTCAGCCGGGCTCGTCGGTGGCGCAAGCACTTAAGCCGCTTCTGGGCCACAGCCCGTTgcgccacctgtt                                                                      >  4:67131/1‑267 (MQ=255)
                                                                                                                                                                                                        cGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATGAtgtg                                                                      >  2:327593/1‑275 (MQ=255)
                                                                                                                                                                                                           cACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑Gt                                                                                                                                                                                                                                                                         >  7:188702/1‑77 (MQ=255)
                                                                                                                                                                                                           cACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑Gt                                                                                                                                                                                                                                                                         <  8:188702/77‑1 (MQ=255)
                                                                                                                                                                                                                  tGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGTGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACCGCCAGTGGCGCATCATGTTCGCGTTCaa                                                             >  1:116263/1‑274 (MQ=255)
                                                                                                                                                                                                                   gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTg                                                                                                                                                                                                                                                                        <  6:1029415/70‑1 (MQ=255)
                                                                                                                                                                                                                   gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTg                                                                                                                                                                                                                                                                        >  5:1029415/1‑70 (MQ=255)
                                                                                                                                                                                                                    cAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctgc                                                                                                                                                                                         <  8:926430/148‑1 (MQ=255)
                                                                                                                                                                                                                     aGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCaa                                                                                                                                                                                                                                                          <  6:124626/82‑1 (MQ=255)
                                                                                                                                                                                                                     aGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCaa                                                                                                                                                                                                                                                          >  5:124626/1‑82 (MQ=255)
                                                                                                                                                                                                                           ggCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTg                                                                                                                                                                                                                                          <  8:646381/92‑1 (MQ=255)
                                                                                                                                                                                                                           ggCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTg                                                                                                                                                                                                                                          >  7:646381/1‑92 (MQ=255)
                                                                                                                                                                                                                                ttcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                                                                                                                                                       <  7:344661/90‑1 (MQ=255)
                                                                                                                                                                                                                                ttcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                                                                                                                                                       >  8:344661/1‑90 (MQ=255)
                                                                                                                                                                                                                                ttcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTg                                              <  2:67132/275‑1 (MQ=255)
                                                                                                                                                                                                                                ttcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTg                                              <  2:255954/275‑1 (MQ=255)
                                                                                                                                                                                                                                   ttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGCCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCCCGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGGCGCGCAAGCCCTTAAGCTGCTGCTGGGCAACAGCCAGTGCCGCACACTGTGCCCGTTCACGATGATGTGGACCtc                                              >  4:151436/1‑271 (MQ=255)
                                                                                                                                                                                                                                        ggcggcggcGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGcatca                                                                                                                                                                                                                                                                 <  7:235788/56‑1 (MQ=255)
                                                                                                                                                                                                                                        ggcggcggcGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGcatca                                                                                                                                                                                                                                                                 >  8:235788/1‑56 (MQ=255)
                                                                                                                                                                                                                                              ggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgtg                                                                                                                                                                                                                      >  8:682003/1‑93 (MQ=255)
                                                                                                                                                                                                                                              ggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgtg                                                                                                                                                                                                                      <  7:682003/93‑1 (MQ=255)
                                                                                                                                                                                                                                                         tcgctATCATCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCg                     <  2:328748/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          cgctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTg                                                                                                                                                                                                         <  6:332689/94‑1 (MQ=255)
                                                                                                                                                                                                                                                          cgctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTg                                                                                                                                                                                                         >  5:332689/1‑94 (MQ=255)
                                                                                                                                                                                                                                                           gctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCt                                                                                                                                                  >  7:575220/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           gctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCt                                                                                                                                                  <  8:575220/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            ctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTCCAGCCTGGCTCGTCGGTGGCCAAAGCCCTTAAGCTGCTGCTGGGCAACACCCAGTGCCGCATCATGTGCTCGTCCAACATGATTGCGACCTTCTCCACCGTGGTTCGCGGCGGGGCGCCg                  >  4:290070/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             tATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACg                  <  2:127502/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              accAGCGGCAGCGCGATGCCGCTG‑‑‑ATGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcg                                                                                                                                                                                                                    <  7:1025917/77‑1 (MQ=38)
                                                                                                                                                                                                                                                              accAGCGGCAGCGCGATGCCGCTG‑‑‑ATGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcg                                                                                                                                                                                                                    >  8:1025917/3‑79 (MQ=37)
                                                                                                                                                                                                                                                               tcagggggtgggcggTGCGGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTCCTTCGGCGCCATGTGCGGGCTGGGGCTGAGCGGCGTGGGGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGct                <  4:124887/263‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccAGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTc                                                                                                                                                                                 <  8:193549/112‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccAGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTc                                                                                                                                                                                 >  7:193549/2‑113 (MQ=255)
                                                                                                                                                                                                                                                                 aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggtg                                                                                                                                                                                             >  5:744373/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                 aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggtg                                                                                                                                                                                             <  6:744373/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg                                                                                                                                            >  7:171459/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  gCGGCAGCGCGATGCCGCTGATAGCGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCt                                                                                                                                                                                                          >  8:633498/1‑88 (MQ=25)
                                                                                                                                                                                                                                                                  gCGGCAGCGCGATGCCGCTGATAGCGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCt                                                                                                                                                                                                          <  7:633498/88‑1 (MQ=25)
                                                                                                                                                                                                                                                                  gCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGt                                                                                                                                           <  8:171459/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCt                                                                                                                                                        <  8:154627/134‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCt                                                                                                                                                        >  7:154627/1‑134 (MQ=255)
                                                                                                                                                                                                                                                                   cGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaa                                                                                                                                                                <  8:412329/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaa                                                                                                                                                                >  7:412329/1‑126 (MQ=255)
                                                                                                                                                                                                                                                                    ggCATCGCGCTGCCGCTG‑‑‑GGGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACtt          <  2:474003/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCg                                                                                                                                                                                                                                <  8:273177/59‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCg                                                                                                                                                                                                                                >  7:273177/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                       aTCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACAcc                                                                                                                                                   >  5:789384/1‑135 (MQ=255)
                                                                                                                                                                                                                                                                       aTCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACCACCAAAGAGCGCTAcccc                                                                                                                                                   <  6:789384/135‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cgcgCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGGGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTg                                                                                                                                                                               <  8:76931/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cgcgCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGCGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTTCGCGTTCAAGCTGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTg      >  2:2091/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                         cgcgCTGCCGCTG‑‑‑ATGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTg                                                                                                                                                                               >  7:76931/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                           cgcTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATTTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGaaa   <  2:411010/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cgcTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCgggg                                                                                                                                  >  5:873290/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                            gcTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGGTGGCTACTTCGGGGCCATGTGCGTGCTGGGGCTGAGCTGCGTGGTGCTGCTCTACGTCTGCTTCGTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAAt  <  2:518711/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                    
GCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAAT  >  NC_000913/285618‑286160

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