Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 286,013 A→C 12.5% T207P (ACC→CCC)  yagG → CP4‑6 prophage, putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913286,0130AC12.5% 245.9 / 3.0 98T207P (ACC→CCC) yagGCP4‑6 prophage, putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base A (46/34);  new base C (5/12);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 5.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCC  >  NC_000913/285786‑286247
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
tGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCACAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GGGCGTGCTGGGGCGGAGGGGGG‑GGGTGCTGCTCTCCGGCTGCTTCTTTCCGACCAAAAAGCGCCCCAccccgctcctt                                                                                                                                                                                                                                  >  2:723325/1‑230 (MQ=255)
    gggCGTCATCCCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGc                                                                                                                                                                                                                              <  1:1044822/239‑1 (MQ=255)
     ggCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGGGGGCTACTTCGGCGCCCT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCCAAGAGCGCTCCACCTTTGAGGTGCAGcc                                                                                                                                                                                                                             >  1:841302/1‑239 (MQ=255)
          catcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGCCGAGCCGGTGGGCTACTTCCGCGCCCT‑GTGCGTGCTGGGCCTGGGCGGCG‑GGGGGCGGCCCTACGCTGGCTTCTTCCCGACCAAAGCGCGCTACCCCTTTGGGGGTCCGCCGGGc                                                                                                                                                                                                                         >  2:945943/1‑238 (MQ=255)
           atcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCTGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGcgc                                                                                                                                                                                                                       >  1:162294/1‑237 (MQ=255)
                  ccgACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAACGGGGACGCGCAGGTGGGCTACTTCGGCGCCAT‑TTGCGTGCTGGGGCTGAGCGGCG‑GGGTGCTGCCCTACGTCTGCTTCTTCACGACCAAAAGACGCTCCACCTTTGAGGTGCAGCCGGGCCcggcgggg                                                                                                                                                                                                                >  2:907559/1‑237 (MQ=255)
                       ccGGAAGGGCGTCACGCACTGCATTCCTGGCGCTTTTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGaa                                                                                                                                                                                                           <  1:1136916/239‑1 (MQ=255)
                       ccGAAAGGGCGTCCCGCCCTGCAGTCCGGGCCCTTCGTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGaa                                                                                                                                                                                                           <  1:945943/239‑1 (MQ=255)
                       ccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGaa                                                                                                                                                                                                           <  1:387004/239‑1 (MQ=255)
                       ccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGaa                                                                                                                                                                                                           <  1:401113/239‑1 (MQ=255)
                               gCGTCACGCCCGGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCt                                                                                                                                                                                                     <  1:225409/237‑1 (MQ=255)
                                cGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGCCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGGGCTGCTCTACGTCTGCTTCTTCACGCCCAAAGAGCGCTACCCCCTTGAGGTGCAGCCCGGCCCGTCCGTGGGGAAAAACCCTaa                                                                                                                                                                                                  >  2:1101022/1‑239 (MQ=255)
                                 gTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCCCTACCCCTTTGAGGGGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTaa                                                                                                                                                                                                  >  1:809899/1‑238 (MQ=255)
                                          tGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGACCGCTACACCTTTGAGGTGCAGCCGGGCTCGGCGGTGGCGAAAGACCTTAAGctactgctg                                                                                                                                                                                        >  1:1154890/1‑239 (MQ=255)
                                            cAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGGGCTGCTCTACGTCTGCTTCTTCACGGCCAAAGAGCGCCACCCCTTTGAGGGGCAGCCGGGCTCCTCGGGGGGGAAAAACCTTAAactgcttggggg                                                                                                                                                                                      >  2:677239/1‑233 (MQ=255)
                                                 cTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAAc                                                                                                                                                                                  >  1:74105/1‑238 (MQ=255)
                                                  tGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑GGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTCCACCTTTGCGGTGCAGCCGGGCTCGTCGGTGGGGAAAAACCCTAAACTGCTGGTggggcaaa                                                                                                                                                                                 >  2:296050/1‑233 (MQ=255)
                                                      gcTTTTTCCTGGCGGGGGGTGGGTCGCCTGGCTATCGGCGGCATCGGCGTGCTGCTGGTGAGCACCATCGGGAAGGGGGAGGAGCAGGTGGGCTTCTTCGGCGCCAT‑GTGCGTGCGGGGGCTGGTCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGGGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCa                                                                                                                                                                             <  2:1102991/239‑1 (MQ=255)
                                                       cttctcCCTGGCGGGGGGGGGCTCG‑CCCGCTATCAGCGGCCTCGCTCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGt                                                                                                                                                                           <  2:587879/239‑1 (MQ=255)
                                                         tcttcCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTCCGTGCTGGGGCTGAGCGGCG‑TGGTCCTGCCCTACGCCTGCTTCTTCACGACCAAAGGGCCCTACACCCTTTAGGGGCGACCGGGCCCGCGGGTGGCGAAAGACCTTAAACCGCGGGCGGGCaaaccaccggtg                                                                                                                                                                         >  2:1111654/1‑229 (MQ=255)
                                                         tcttcCCGGCGGCGGCGGGCTCG‑CTCGCTAGCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTgg                                                                                                                                                                         <  1:725503/239‑1 (MQ=255)
                                                                      gcggggCTCG‑CTCGCTTTCGGCGGCCTCGCGCTGCCGCTGGTGGGCCTCATGGGCAAAGGGGACGGGCAGGTGGGCTCCTTCGGCGCCAT‑GTGCGTGCTGGGGGTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTgcg                                                                                                                                                            <  2:727345/236‑1 (MQ=255)
                                                                       gcggGCTAG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAAAGCCCTACCACCTTGGGGTGCAGCCGGGCTCCCCCGGGGCGGAAGACCCTCACCGGCTTCTGGGCAACCGcccgggcggcacacgtggccg                                                                                                                                                           >  2:1079585/1‑220 (MQ=255)
                                                                                        taaGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCg                                                                                                                                          <  1:826483/237‑1 (MQ=255)
                                                                                          gggggcacggggccgCCGCTGGGGGGCATCGACGGCAAAGGGGAGGGGCAGGTGGGCTACTTCGGCGCCTG‑GGGCGTGCTGGGGCTGGGCGGCG‑GGGTGCTGCTCGACGCCGGCTTTTTCACGACCAAAGAGCGATACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGAc                                                                                                                                        <  2:1243221/225‑1 (MQ=255)
                                                                                          cggggctccgcgtcgCCGGTGGTGGGCTGCTTCGGCAAGGGGGACGGGGAGGGGGGCTACTTCGGCGCCAT‑GTGGGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGAc                                                                                                                                        <  2:1211279/225‑1 (MQ=255)
                                                                                             ggcctggcgcTCCGGGTGGGGGGCATCGGCGGCAAGGGGGAGGGGGGGGGGGGCTTCTTCGGGCCCAT‑GGGGGTCCTGGGGCTGGGCGGGG‑GGGTGCTGCTCTACGTCTGCTTGTTCACGACCAAAGAGCGCTACCCCTTTGAGGTGCAGCGGGGCGCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTg                                                                                                                                     <  2:1344849/233‑1 (MQ=255)
                                                                                             ggCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGGAGGTGGGCTACTTCGGCGCCAT‑GGGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGAc                                                                                                                                        <  1:901278/236‑1 (MQ=255)
                                                                                                aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGCCCTGCTc                                                                                                                                  >  1:63869/1‑239 (MQ=255)
                                                                                                aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGcc                                                                                                                                   >  1:901299/1‑237 (MQ=255)
                                                                                                aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCt                                                                                                                                   <  1:1309582/238‑1 (MQ=255)
                                                                                                    cgcTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGGCACACCCAGTGGCGCATCAT‑TGCGCCGTTAAAAGTGATGGGGAccgccccaaag                                                                                                                              >  2:103828/1‑230 (MQ=255)
                                                                                                       tGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACgtg                                                                                                                           <  1:1006275/239‑1 (MQ=255)
                                                                                                               gtgcgctccacgggCAAAGGGGACAAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCAACGTCTGCTTCTTCACGACCAAAGAGCGCTACCCCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGcgcg                                                                                                                    <  1:887734/226‑1 (MQ=255)
                                                                                                                     atcatgGGCAAGGGGGGCGGGGTGGTGGGCTACCTGGGCGCCAT‑GTGCGTGCGGGGGCTGGGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACCCCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCgggg                                                                                                              <  2:809899/238‑1 (MQ=255)
                                                                                                                     atcatcGGCAAAGGGGACGGGCAGGGGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACCCCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCgggg                                                                                                              <  1:1113972/238‑1 (MQ=255)
                                                                                                                     atcatcGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAG‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCgggg                                                                                                              <  1:1108400/238‑1 (MQ=255)
                                                                                                                         tcGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGCCCTTAAGCTGCTGCTGGGCAACAGCACGGGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGCACTTCCCCACCGGGGGGCGCGGCGGGGCGAc                                                                                                          >  2:1111519/1‑238 (MQ=255)
                                                                                                                         tcGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCCGGGCGAc                                                                                                          >  1:522978/1‑238 (MQ=255)
                                                                                                                         tcGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGAGGATGCCGACCTGCTCCAACGGGGTGCGCGGCGGGGCGAc                                                                                                          >  1:14262/1‑238 (MQ=255)
                                                                                                                         tcGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGCCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGAc                                                                                                          >  1:1005836/1‑238 (MQ=255)
                                                                                                                              aaaGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCCACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAAGTGATGGGGCCCCGCCCCAACGTGGTGCCgggggggggcacggccc                                                                                                     >  2:599567/1‑229 (MQ=255)
                                                                                                                               aaGGGGCCGAGCAGGTGGGCTCCGTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCCCGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTAc                                                                                                   <  1:1087542/239‑1 (MQ=255)
                                                                                                                                aGGGGACGAGCAGGTGGGCTTCTTCGGCGCCAT‑GTGCGGGCTGGGGCTGGGCGGCG‑TGGTGCTGCTCTACGTTTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACt                                                                                                  <  1:296050/239‑1 (MQ=255)
                                                                                                                                  gggCCGGGCGGGGGGGCTCTTGCGGCGCCAT‑GTGCGGGTTGGGGGTGGGCGGGGTTGGTG‑TGCTCTTCGTCGGCTTCTTCACTCCCAAAGAGCGCTCCACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACt                                                                                                  <  2:541155/237‑1 (MQ=255)
                                                                                                                                        gggggggggggttcctCGGCGCCATGGTGGGTG‑TGGGGGGGAGGGGGG‑GGGTGCTGCTCTGAGTCTTCTTCTTCCCGACCAAAGGGGGCTACACCTTTGAGGGGCCGCCGGGCTCGGCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGaa                                                                                           <  2:43299/223‑1 (MQ=255)
                                                                                                                                         gctggTGGGCTATTTCGGCGCCTT‑TTGCGTGCTGGGGTGGGGCGGGG‑TGGTGCTGCTTGCCGTCGGCTTCTGCACGACCAAAGGGCGCTACGCCTGTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGaaa                                                                                          <  2:14262/235‑1 (MQ=255)
                                                                                                                                         gCAGGTGGGCTACTTCGGCGCCAT‑GGGGGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGTTACCCCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGaaa                                                                                          <  1:1026459/238‑1 (MQ=255)
                                                                                                                                          cgggTGGGGTCCTTCGGCGCCAT‑GTGCGTGCTGGGGGTGGGCGGGG‑TGGTGCTGCTCTACGTTTGCTTCTTCACGACCAAAGAGCGCTACCCCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAAt                                                                                         <  2:901299/236‑1 (MQ=255)
                                                                                                                                           aGGTGGGCTACGTCGGCGCCAT‑GGGGGTGCTGGGGCTGGGCGGCG‑GGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACCCCTTTGAGGTGCAGCCGGGCTCGGCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATa                                                                                        <  1:850596/238‑1 (MQ=255)
                                                                                                                                            ggTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGTTACGCGGCGGGGCGCCGCTCTACTTccggagatccg                                                                                      >  1:717186/1‑229 (MQ=255)
                                                                                                                                             gTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAACCTGCTGCGTGGCAACAGCCAGTGGCGCATCATGTGCG‑CGGTCAAGATGGTGGGGACCTGCCCCACCTTGGTGGCCGGCGGGGCGCCGCCCTACTTCGGGAAAAACGt                                                                                     >  2:779636/1‑239 (MQ=255)
                                                                                                                                                 gCTCCTTCGGCGCCAT‑GGGGGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGAt                                                                                  <  1:797522/238‑1 (MQ=255)
                                                                                                                                                     tccgtggcgcaattgtctggctggggctgggCGGCG‑TGGTGCTTCTATGCGTCGGCTTCTTCCGAACCAAAGAGCCTCACCCCTTGGGGGTCACGCCGGGCTCGTCGGTGGCGAAAGGCCTTAAGCTGCTGCTGGGCACCAGCCAGTTGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGAt                                                                              <  2:294036/220‑1 (MQ=255)
                                                                                                                                                          gcgcCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGGGCG‑CGTTCAAGATGATGGCGACCTGCCCCCAGGTGGGGCGTGCGGGGGCGACCCCCTACTTCGGGAAATACGTGAGGGAGCcccca                                                                        >  2:1237326/1‑238 (MQ=255)
                                                                                                                                                             ccAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCACCGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAg                                                                     >  1:869893/1‑239 (MQ=255)
                                                                                                                                                                t‑gtgCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAGAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGtt                                                                                                                                                       <  1:133521/154‑1 (MQ=255)
                                                                                                                                                                t‑gtgCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAGAGACCTTAACCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGtt                                                                                                                                                       >  2:133521/1‑154 (MQ=255)
                                                                                                                                                                       tgtgggggTTGGGGGGCG‑TGGTGGTGCTCTACGTCTGCTTCTCCGCGACCAAAGAGCGCTACACCTTTGAGGGGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGAc                                                             <  1:779636/234‑1 (MQ=255)
                                                                                                                                                                       gtttataagagAGGGGCG‑TGGTGCTGCGTTACGTCTGCTTCTGCCCGACCAAAGGGCGCTGCCCCTTTGGGGGGCAGCCGGGCTCGTCGGTGGCGAAAGCCCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGAc                                                             <  2:22544/228‑1 (MQ=255)
                                                                                                                                                                        gctggggctgggCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTGGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGAcc                                                            <  1:1202108/238‑1 (MQ=255)
                                                                                                                                                                             ggCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTa                                                                                                    >  1:455188/1‑193 (MQ=255)
                                                                                                                                                                                tGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATgg                                                                                <  1:972031/210‑1 (MQ=255)
                                                                                                                                                                                  aGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGAGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGTTGGATCACCCGGAGTTGGCGACCCAGTTtttact                                                 >  1:193019/1‑239 (MQ=255)
                                                                                                                                                                                    cggcg‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTtttacttt                                               >  1:756578/1‑239 (MQ=255)
                                                                                                                                                                                           ggtgCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGGGCGCGGCGGGGCGACGCCCTACTTCGTGAAATACGCGATGGATAACCCGGGGCTGGcccttataaccatctgacgccgac                                         >  2:972031/1‑216 (MQ=255)
                                                                                                                                                                                            caccgggttccccgtgctttcttcACGACCAAAGAGCGCTACCACTTGGGGGGGCAGGCGGGCTCGTCGGGGGCGAAAGACCTTAAGCTGCTGGTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGc                                          <  2:699446/219‑1 (MQ=255)
                                                                                                                                                                                             gcctgctCTACGTCGGCGTCCTCACGACCAAAGAGCGCTACACTTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCGGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACg                                                                                                                             <  2:1053750/151‑1 (MQ=255)
                                                                                                                                                                                                  ctctACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGGCCCAGTTTTTACTTTACGGCAGCCttg                                   >  1:1117026/1‑236 (MQ=255)
                                                                                                                                                                                                  ctctACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCCGTTtttacctgcctcttttacac                                   >  1:22544/1‑223 (MQ=255)
                                                                                                                                                                                                  ctctACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGGCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGAc                                                             <  1:553813/212‑1 (MQ=255)
                                                                                                                                                                                                    ctACGTCTGCTTCTTCACGACAAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCCCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGACCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGcca                                <  1:677239/239‑1 (MQ=255)
                                                                                                                                                                                                       cGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCCTCCTGTGCG‑CGTTCAAGATGATGGCGACCGGCTCCAACGTGGTGCGCGGCGGGGCGACGCCCTACTTCGTGAAATACGGGATGGATCACCCGGAGTTGCGGGCCCAGTTTTTACTTTACggagcccccgcccc                               >  2:1255272/1‑225 (MQ=255)
                                                                                                                                                                                                       cGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCgcgg                                                                                                                   >  1:1053750/1‑153 (MQ=255)
                                                                                                                                                                                                       cGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCCGTTTTTACTTTACGGCAGCCTCGccacca                             >  1:19072/1‑239 (MQ=255)
                                                                                                                                                                                                          ctgtttcctccgggccAAAGTGCGCTAGCCCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑GGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGt                          <  1:1237326/225‑1 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGCTGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGCAGCGCTACTTCGTGAAATACGGGATGGAAACCCCGGGGTTGGGGACCCAGTTtttacctttgccttaaaaacatctgccgctg                        >  2:553813/1‑216 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGCCGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCTCCCCCATGTTc                        >  1:878103/1‑239 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGGGACCCCGTTTTTACTTTACGGCAGCCTCGCCACCATGTTc                        >  1:591606/1‑239 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCCGTTTTTACTTTACGGCAGCCTCGccacctgtc                          >  1:699446/1‑233 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGcccccagtttc                        >  1:737055/1‑234 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCCCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGGGACCACGTTTTTACTTTACGGCAGCCTCGCCCCCATGTTc                        >  1:480051/1‑239 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCCCCAACGGGGTGCGCGGCGGGGCGACGCCCCACTTCTTGAAAACCGTGTTTGATCACCCGGAGTTGGCGACCCACTTTTTACTTTACGCCAccccccccccccagttc                        >  2:1112380/1‑225 (MQ=255)
                                                                                                                                                                                                            gCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATAGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGCCGCTCTACTTCGTGAAATAAGGGATGGATCACCCGGCGTTGGCGACCCAGCTTTTCCTTTACGGCAGCCCCGCCCCCCAGTTc                        >  2:308713/1‑239 (MQ=255)
                                                                                                                                                                                                             cttcttCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTg                                                                                                                                     <  2:342050/129‑1 (MQ=255)
                                                                                                                                                                                                             cttcttCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTg                                                                                                                                     >  1:342050/1‑129 (MQ=255)
                                                                                                                                                                                                             cttcttCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGGCGCTCTACTTCGTGAAATACGTGAAGGACCCACCGGCGTGGGCCACCCAGCTTTTACTTTACGGCAACCTCCCCACCAGGTTCg                       >  2:351662/1‑239 (MQ=255)
                                                                                                                                                                                                               tcttcACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCAtgtg                                                                                                                                                              <  2:783592/103‑1 (MQ=255)
                                                                                                                                                                                                               tcttcACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCAtgtg                                                                                                                                                              >  1:783592/1‑103 (MQ=255)
                                                                                                                                                                                                                    aCGCCCAAGGGGCGCTACCCTTTTGGGGTGCAGCTGGGGTGGGCGGGGGCGAAAGACCTTAAGCGGCTGCGGGGCAACAGCCAGTGGCGCATCAGGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCt                <  2:63869/239‑1 (MQ=255)
                                                                                                                                                                                                                        ccAAAGAGCGCTACACCTTTGAGGTGCAGGCGGGCTCGGCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTg             <  1:103828/238‑1 (MQ=255)
                                                                                                                                                                                                                          aaaGGGCGCTACACCTTTGGGGTGCAGCCGGGCTCGTCGGTGGCGAAAGGCCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGctc          <  1:939945/239‑1 (MQ=255)
                                                                                                                                                                                                                           aaGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGGGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGGGATGGATCACCCGGGGTTTGCGACCCAGTTTTTACTTTACGGCACCCCCGCCCCCAAGTTTGGCCCGCTTTGccc          >  2:1262527/1‑236 (MQ=255)
                                                                                                                                                                                                                                  cTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGGCCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGc   <  1:120261/238‑1 (MQ=255)
                                                                                                                                                                                                                                  cTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCGGTTTTTACTTTACGGCACCCTCGCCACCAAGTTCGGCCCCCTTTGCCCCCCACGcc  >  2:120261/1‑239 (MQ=255)
                                                                                                                                                                                                                                  cTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCAccc   >  1:956118/1‑236 (MQ=255)
                                                                                                                                                                                                                                  cTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCGTTCCTTCGTGAAATACGTGAGGGATAACCCGAGGTTGGGGACCCAGTTTTTACTTTACGCCAGCCCCCCCACCCAGTTCGCCCCGCttttcacccccccc   >  2:42700/1‑227 (MQ=255)
                                                                                                                                                                                                                                  cTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCCCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCCGTTTTTACTTTACGGCACCCTCGCCCCCATGTTGGGCCCGCTTTGCTCCTCCCg    >  2:1184470/1‑237 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
TGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCG‑CTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAT‑GTGCGTGCTGGGGCTGAGCGGCG‑TGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCG‑CGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCC  >  NC_000913/285786‑286247

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑