Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,458,425 A→C 20.9% T156P (ACC→CCC)  pmbA → putative antibiotic peptide MccB17 maturation peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,458,4250AC20.9% 166.5 / 8.8 91T156P (ACC→CCC) pmbAputative antibiotic peptide MccB17 maturation peptidase
Reads supporting (aligned to +/- strand):  ref base A (45/24);  new base C (3/19);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 2.27e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TAGCGCATCATCCACCGATTTAAGCCCGCAGGCCATTGCCCGTACTGTACAGGCGGCGCTGGATATTGCCCGTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGT  >  NC_000913/4458190‑4458658
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
gggcctagcctcCCCCGATTTTAGCCCGCAGTCATTTGCCCGGACTGTACAGGCGGCGCTGGATATTGCCGTGTAGTCCCCCCCAGATCCCTGTGCCGGGGTGGCAGACAAAGAGCTGCTGGCCTTGGACGCACCAGATCT‑‑CGACTTGTTCCCCCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCAc                                                                                                                                                                                                                                               <  2:374567/227‑1 (MQ=255)
         atcCACCGATTTAAGCCCGCAGGCCATTGCCCGTACTGTACAGGCGGCGCTGGATATTGCCCGTTATACCTCGCCAAATCCCTGGGCCGGCGGGGCAGAACCAAAGCTGCTGGGCCTTGAAGCGCCCGACTC‑‑CGACTTGTTCCACCCCCCGCCAAGTTCCCCGGGTGAAGCCCTTGAACCGGCGGCCCGGCCAGAAAAGGCGGCATGTCAGGGCGGAAAAAACACCACCCAAACCGAAg                                                                                                                                                                                                                                  >  2:1305969/1‑239 (MQ=255)
         atcCACCGATTTAAGCCCGCAGGCCATTGCCCGTACTGTACAGGCGGCGCTGGATATTGCCCGTTATACCTCCCCAGATCCCTGTGCCGGCGTGGCCGACAAAGAGCTGCTGGCCTTTGCCGCACCAGCTCT‑‑CGGCTTGTTCCACCCTGCGGAAGTTTCCCCGGGTGAAGCCCTTTAACAGGCGGCCCGCGCAGAAAAGGCGGCACTGCGGGGGGCGCAACACAACACCAAAACCAAAg                                                                                                                                                                                                                                  >  2:100042/1‑239 (MQ=255)
            cACCGATTTAAGCCCGCAGGCCATTGCCCGTACTGTAAAGGCGGGGCTGGATATTGCCCGTTAGTCCTCGCCAGATCCCTGTGCCGGCGTGGCAGGAAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTg                                                                                                                                                                                                                              <  1:382425/239‑1 (MQ=255)
                          cGCAGGCCATGGCCCGTACTGTACGGGCGGCGCTGGATATTGCCCGTTATCCCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCCCCAATACCGAAGGTGGCAGCTTTAACAGc                                                                                                                                                                                                                <  1:968650/239‑1 (MQ=255)
                                ccATTGCCCGGACTGTACAGGCGGCGCTGGATATTGCCCGTTATCCCTCGCCAGATCCCTGGGCCGGCGTGGCAGACAAAGAGCTGCGGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTa                                                                                                                                                                                                          <  1:100042/238‑1 (MQ=255)
                                 cATTGCCCGTACTGTACAGGCGGCGCGGGATATTGCCCGTTATCCCTCGCCAGATCCCTGGGCCGGCGTGGCAGACAAAGGGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCCCCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACg                                                                                                                                                                                                        <  1:808761/239‑1 (MQ=255)
                                 agaagcgggtccTGTACGGGCGGCGCTGGATATTGCCCGTTATCCCTCGCCAGACCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACg                                                                                                                                                                                                        <  1:1065427/228‑1 (MQ=255)
                                           aCTGTACAGGCGGCGCTGGATATTGCCCGTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGGCAAACGCATTACCAAAACCGGAGGGGGAAGCTTTAACAGCCACCACCGTGCCAAAgg                                                                                                                                                                                              >  1:1183209/1‑238 (MQ=255)
                                               tACAGGCGGCGCTGGATATTGCCCGTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCCTTTCAGGCGGACAAAACCATCCCCAAAACCCAAAGAGGGAGCTTTAAAAACCACAACgcggtgaaaagtttt                                                                                                                                                                                          >  2:947897/1‑225 (MQ=255)
                                                    gcggcgCTGGATATTGCCCGTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCCCGCCGAACAAGCGGCATTTCAGGCGGACAAAAGCCACACCAATACCGAAAGAGGGAACCTTAAAACCCACCACAGCGGCaaaaa                                                                                                                                                                                              >  2:1065427/1‑228 (MQ=255)
                                                     ggcgttgggtttgtGCCTGTTATCCCCCCCCAGGCCCCTGGGCCGGGGGGGCAAACAAAGGGCTGTGGGCCTTGGCCCCCCCAGTACG‑‑CGACTTTCCCCCCCCTGGGGAGTTTTCCCCGGATGAGGCCATTGAACGGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAACCGCATCACCAATCCCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAAc                                                                                                                                                                                    <  2:219678/226‑1 (MQ=255)
                                                         gcTGGATATTGCCCGTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGGTTCCCCGGATGAAGCCATTTAACTGGCGGCCCGCGCCGAACACGCGGCATTGCAGGCGGGCAAACCCCTCCCCAACACCCGAAGTGGGCGCTTTTAAAACCCCCAAGGGTGCAAAATTTTTTGCCAACaca                                                                                                                                                                                 >  2:467301/1‑237 (MQ=255)
                                                            ggtgtttGCCCTTTATCCCCCACCAACCCCCGTGGCGGGGGGGCCAGAAAAAGCGTTGCGGGCTTTGGACGCACGAGACCC‑‑CGACTGTTGCCCCCCGGCGGACGTTTCCCCGGATGAAGCCATTGACCGGGCGGCCCGCGCAGAACGGGCGGCTTTGCCGGCGGAAAACCGAATCCCCAATACCGAAGGTGGCGGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACg                                                                                                                                                                             <  2:1338362/234‑1 (MQ=255)
                                                               tttttccctttttcccTCCCCAGCTCCCTGGGCGGGGGTGGCAAACAAAGGGCGTCGGGCCTTTGACGCCCCAGATCC‑‑CGACTTGTCCCACCCGGCGGAAGTTTCCCGGGAGGAAGCCATGGACCGGGCGGCCCGCGCAGAACAGGCGGCTTTGCAGGCGGACAAACGCATCACCAATACCGAAGGGGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGc                                                                                                                                                                           <  2:729744/224‑1 (MQ=255)
                                                                    cccGTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCAACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCCCTACGGTGTCAAAGTTTTTGGCAACCGCCACGGCATGtt                                                                                                                                                                      >  1:862187/1‑238 (MQ=255)
                                                                       gTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCCATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAg                                                                                                                                                                 >  1:1298183/1‑239 (MQ=255)
                                                                        ttATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACCAACGCATCACCAATACCGAAGGGGGCAGCTTTAAAAGCCACTACGGTGTCAAAGATTTTGGCACAAGCCACGGCATGTTTCAgg                                                                                                                                                                >  1:23384/1‑239 (MQ=255)
                                                                        ttATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAAACCCTCACCAATACCGAAGGGGGCAGCGTTAACAGCCACTACGGTGGTAAAGATTTTGGCAACAGGCACGGCATggtcaaag                                                                                                                                                                 >  1:22509/1‑232 (MQ=255)
                                                                         tataCCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTTGCCTTTGACGCACCAGATCT‑‑CGGCTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCCTTGAACTGCCGGCCCGCGCAGAACAAGCGGCATTGCAGGCGGACAAACGCCGCACCAATACCCAAAGTGGCGGGATTAACAGCCAATCAGGGGGTAAAAGTTTTGGCaaaaacaacgcacaggtagcagg                                                                                                                                                                >  2:450654/1‑218 (MQ=255)
                                                                          atccccagccggcccccGGTGCGGGGGTGGAAGACAAGAGGCTGCGGGCTTTTGACGACCCAGACCC‑‑GGATGTCTCCCCCCGGGAGGAAGTTCCCCGGAAGGAAGCCATTGAACGGGCGGCCCGCGCAAAACGGGCGTCATTGCAGGCGAACAAACGCATCCCCAATATCGAAGGGGGCAGCTTAAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGt                                                                                                                                                              <  2:554399/225‑1 (MQ=255)
                                                                                 gCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGCAAAACCCCCCACCAATACCGAAGGAGGCAGCTTTAAAAACCACCACGGTAGCAAAGTTTTTGGCAAAAGCCCCGGCATGTTCGCAGGTTACttc                                                                                                                                                        >  2:27960/1‑236 (MQ=255)
                                                                                 gCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGGGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTTTTGGCAACACCCACGGCAGGTTTCACGGTTACTGcc                                                                                                                                                       >  1:1277690/1‑238 (MQ=255)
                                                                                  ccagtccccTGTCCCGGCTTGCCAGAAAAAGGGCTGCGGCTCTTGGACGCACCAGATCC‑‑CGCTTTGTCCCCCCCGGCGGAAGTTCCCCGGGAGAAACCCATTGAACTGGCGGCCCGCGCAAAACGGGCGGATTTCCAGGGGGACAACCGCATCCCAAATACCGAAGTGGGCAGTTTAAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTc                                                                                                                                                      <  2:621164/233‑1 (MQ=255)
                                                                                     gATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTTTTGGCAACAGCCACGGCCAGTTGCAGGGGTACGGCCCAAc                                                                                                                                                   >  1:387858/1‑239 (MQ=255)
                                                                                             tgCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCCCTACGGGGTCAAAGTTTTTGGCCACAGCCACGGCATCTTGCAGGGTAACTCCt                                                                                                                                                       >  1:1225021/1‑227 (MQ=255)
                                                                                             tgCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCCAAGTTTTTGGCAACAGCCCCGGCATGTTGCAGGGGTACTGCTCAACCCGtatt                                                                                                                                            >  1:78163/1‑235 (MQ=255)
                                                                                             ggCCGGCGGGGCAGACAAAGAGCTGCTGGCCTTGGACCCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCAtt                                                                                                                                           <  1:366156/238‑1 (MQ=255)
                                                                                             gcccggtctGGCAAAAAGAGGGCTGCTGCCTTTGCGCGCCCGATCTTCCGACTTGTCCCACCCGGGGGAATTTTCCCGGGATAAAGCCATTAACCGGGCGGCCCGCGAAAAACGGGCGCCTTGGCAGGCGAAAAAACGCACCCCCAATACGGAAGGTGCAAGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCAtt                                                                                                                                           <  2:599486/231‑1 (MQ=255)
                                                                                              ggccggctttcaggcaaaaggtctgctgCCTTTTGCCGCATCAGGTCT‑‑GGATTTGTCCCCCCTGGGGGAATTTTCCCCGGATGAGGCAGTTAACCGGGGGGCCCGGGAAGAACAGGCGGACTTGCAGGCGGACAAACGCATCACCAATACCGAGGGGGGACGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCAtt                                                                                                                                           <  2:72385/217‑1 (MQ=255)
                                                                                              gCCGGCGTGGCAGACAAAGAGCTGCTGCCCTTTGACGCCCCAGATCC‑‑CGACTTGTCCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCAtt                                                                                                                                           <  1:947897/238‑1 (MQ=255)
                                                                                                  gCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCCTTGCAGGCGGACAAAAGCATCAACAAAACCGAAGGGGGCCGCTTTGAAACCCACTACGGGGTGAAAATTTTTGGGAACACGCAAGGGCTGTTGCAAGGTTACTGCCCAACGGGCCATTCGcc                                                                                                                                       >  2:1071287/1‑237 (MQ=255)
                                                                                                  gCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGct                                                                                                                                       >  1:326979/1‑238 (MQ=255)
                                                                                                  gCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCACCACCAATACCGAAAGGGGCAGCTTTAAAAGCCACCAAGGGGTCCAAGTTTTTGGCAACAACCCAGGCCAGGTGCCGGGGTTCAGCCCAAAGCGGCATTcccc                                                                                                                                       >  2:768470/1‑235 (MQ=255)
                                                                                                   cGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGGGTCATTCGccc                                                                                                                                      >  1:864944/1‑236 (MQ=255)
                                                                                                   cGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACCCCTCCCCAATACCCAAGGTGGCAGCTTTAACAGCCCCTACGGTGTCAAAGTTTTTTGCAACAACCACGGCCAGTTGCAGGGGTACTGCTCAACACGTCAATcctcc                                                                                                                                      >  2:628481/1‑234 (MQ=255)
                                                                                                                      gctgGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCCGCTGTGTAATTTCCg                                                                                                                   >  1:870985/1‑239 (MQ=255)
                                                                                                                       ctgGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGCCCAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCCATCCGGGGTCAAAGTTTTTGGCAACAGCCACGGCCTGTTGCCGGGTTACTGCTCAACGCGTCATTTGCTCTCCCGCTGTGTAATTGccca                                                                                                                  >  2:789421/1‑237 (MQ=255)
                                                                                                                        tgGCTTTTGCCCCACTAGATCT‑‑CGCCTTGCTCCCCCCGGGGGAGTTTCCCCGGGAGGAAGCCGTTGAGCGGGCGGCCCGCGAAGAACGGGCGGCATTGCAGGCGGACAAACGCACCACAAATCCCGAAGGTGGCAGTTTAAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCgaa                                                                                                                 <  2:387858/239‑1 (MQ=255)
                                                                                                                             tttGACGCACCAGATCC‑‑CGACTTGTCCCCCCCGGGGGAGTTTTCCCGGAGTGAAGCCATTGACCTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTAAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGaaaa                                                                                                            <  1:789421/239‑1 (MQ=255)
                                                                                                                               tGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGGGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTTTTGGCAACAGCCCCCGCAGTTTGCAGGGGTACTGCCAAACGCGTCATTTCCCCTCCACCTGTGGAATTGCCGGAGaaaaag                                                                                                          >  2:841631/1‑237 (MQ=255)
                                                                                                                                   gCACCAGATCT‑‑CGACTTGTTCCCCCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCg                                                                                                       <  1:1071287/238‑1 (MQ=255)
                                                                                                                                     aCCAGCTCC‑‑CGACTTGTTCCCCCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATCCCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGaaa                                                                                                    <  1:522466/239‑3 (MQ=255)
                                                                                                                                     aCCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCACACCGTCATTCGCGCTCCAGCTGTGTAATTGCCGAGGAAAATGGCGaca                                                                                                    >  1:1285730/1‑237 (MQ=255)
                                                                                                                                     aCCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGGGTCAAAGATTTTGGCAACAGCCAAGGCATGTTGCAGGGGTACCGCTCAACGCGTCATTTGCTCCCCAGCCGGGGAATTGCCGAAAAAAAAgggcaaa                                                                                                    >  2:645203/1‑234 (MQ=255)
                                                                                                                                     aCCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAAAAGCCACTAAGGGGGCAAAGTTTTTGGCAACAGCCCCGGGTTGTTGCAGGGGTACTGCCCACCCCGTCATTCCATCACCCCCCGTGTAATTGACGAAGaaaaaaggggta                                                                                                    >  2:1121591/1‑231 (MQ=255)
                                                                                                                                      ccAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTTCCGAAGAAAATGGCGatat                                                                                                   >  1:879439/1‑239 (MQ=255)
                                                                                                                                      ccAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGatat                                                                                                   >  1:754288/1‑239 (MQ=255)
                                                                                                                                      ccAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTTTAATTGCCGAAGAAAATGGCGatat                                                                                                   >  1:751043/1‑239 (MQ=255)
                                                                                                                                             tcccgctttttcccACCCGGCGGAAGTTCCCCCGGATGAACCCATTGAACTGGCGGCCCGGCGAGAACGGGCGGCATTGCAGGGGGACAAACGCATCCCCAATCCCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGc                                                                                              <  1:1019611/227‑1 (MQ=255)
                                                                                                                                              tcccgcTTGTCCCCCCCCGGGGGAATTTCCCCGGGAGGCACCCATTAAAGCGGGCGCCCGCGAAAACAGGGGGGAATTGAAGGCGAACAACCGCATACCCAATACGGAGGGGGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGc                                                                                              <  2:805964/234‑1 (MQ=255)
                                                                                                                                              cgcctcTTGTTCCCCCTGGGGGAGGTTTCCCGGGGTGAGGCAATTGAACTGGCGCCCCGCGCAAAAGGGGCGGATGGACAGCGGAAAAAACTCTACCCCAATCCCGAAGGTGGCAGCTTAAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCAGTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCg                                                                                             <  2:19438/234‑1 (MQ=255)
                                                                                                                                                  ccttgtccccccggggggAGGTTCCCCGGGAGGAACCCATGGAACGGCCGACCCGCGCGAACAGGGCGGAATGCCGGGCGAAAAAAGCATTCACAAATCCGAAAGGTGGCAGCTTAACCAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGAt                                                                                         <  2:751043/223‑1 (MQ=255)
                                                                                                                                                  atttGTCCCACCCTGCGGAAGTTTCCCCGGATGAGGCCATTGAACTGGCGGCCCGCGCAGAACAGGGGGCATTGCAGGAGGACAAACGCATCACCAATCCCGAAGGTGGCAGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGAt                                                                                         <  1:841631/237‑1 (MQ=255)
                                                                                                                                                   cTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCCGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCACAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGATCaa                                                                                                                                                    >  1:1210829/1‑179 (MQ=255)
                                                                                                                                                      gTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGGCAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTATTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGGTTCATTCGCTCTCCCGCGTGGTAATTGCCGAAAAAAAAGGCAAAATGGAGGGGGACTAc                                                                                      >  2:1186404/1‑238 (MQ=255)
                                                                                                                                                      gTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGCCAAACGCATCACCAATACCGAAGGGGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTACACGCGTCCTTCGCTCTCCAGCCGGGTACTTGCCCAACACAAAGGGGCTATGGAGCGTGATAAc                                                                                      >  2:101107/1‑238 (MQ=255)
                                                                                                                                                          cccccgggggagttctCCCGGGGTGAACCATTTAAACGGGCGGCCCGGGCAAAACGGGCGCCATTGCAGGCGGAAAACCGCACCCCAAATCCCGAAGGTGGAAGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTa                                                                                 <  2:864944/224‑1 (MQ=255)
                                                                                                                                                            cccTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCCCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTCCGCCTac                                                                                >  1:1279204/1‑238 (MQ=255)
                                                                                                                                                                gCGGAAGTTCCCCCGGATGAAGCCATTGACCTGGCGGCCCGCGGAGAACAGGCGGAATTGCAGGCGGACAAACGCATCCCAAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCAt                                                                           <  1:628481/239‑1 (MQ=255)
                                                                                                                                                                 cGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTa                                                                                                                                                            >  1:387577/1‑157 (MQ=255)
                                                                                                                                                                 cGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGCTATGGAGCGTGATTACGCCTAAACCAtt                                                                          >  1:863513/1‑239 (MQ=255)
                                                                                                                                                                       ttcccccGGTTGAGCCCATTGAGGGGGCGGCCGCCAGCGAAGAGGGGGTATTGCAGGCGGACAAACGCAACCCCAATGCCGAAGGTGGCAGTTTAACCAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCCAGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTTGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGt                                                                    <  2:544437/236‑1 (MQ=255)
                                                                                                                                                                        tcccccGGATGAAGCCATTGACCTGGCGCCCGGCGAAGAACAGGCGGCATTGCAGGCGGAAAACCGCAACCCCAATACCGAGGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGt                                                                    <  1:506005/236‑1 (MQ=255)
                                                                                                                                                                         tCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCACCAGCCCCGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCACCTGTGTAATTGCCGAAGAAAAAGGGGATATGGAGCGGGATTTCGCCTACAACAtttggcctgc                                                                  >  2:1096891/1‑231 (MQ=255)
                                                                                                                                                                         tCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTAATGCTCACAGCGTCATTCGCTCTCCAGCTTGGTAATTGCCGAAAAAAAAGGGGAAATGGGGGGGGATTAACCCCACACCATggggggggg                                                                  >  2:719325/1‑233 (MQ=255)
                                                                                                                                                                           cccGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACGGGCGGCATTGCAGGCGGACAACCGCATCACAAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGa                                                                <  1:645203/239‑1 (MQ=255)
                                                                                                                                                                                gggtGAAGCCATGAAACTGGCGGCCCGCGAAAAACAGGGGGCATTGCAGGGGGAAAAACGCATCCACAATTCCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGa                                                             <  1:569664/236‑1 (MQ=255)
                                                                                                                                                                                      gcccATTGAACTGGGGGCCCGCGCAGAACGGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGAt                                                        <  1:1032787/236‑1 (MQ=255)
                                                                                                                                                                                       gCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAAGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCCGGGTAATTGCCGAAGAAAATGGCGATATGGGGCGTGATTACCCCTACACCATTGGTCGTGCCAtagggga                                                         >  2:341168/1‑230 (MQ=255)
                                                                                                                                                                                          aTTGACCGGGCGGCCCCGGCAGACAGGGGGGCATTGCAGGCGGACAAACGCAGCACCAATCCCGAAGGGGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCa                                                   <  1:768470/239‑1 (MQ=255)
                                                                                                                                                                                          aTTGAACTGGCGGCCCGCGCAGAACGGGCGGCATGGCAGGCGGAAAAACCTTCCCCAAATCCCGAAGGTGGCAGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCa                                                   <  1:339155/239‑1 (MQ=255)
                                                                                                                                                                                            tGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGGTACTGCCCAACGCGTCATTCCCTCTCCAGCCGGGTAATTGCCGAAGAAAAAGGGGATATGGAGCGGGGTTACGCCTCAACCCTTGGtggggtgggaacggatgggcaaa                                                 >  2:1164824/1‑217 (MQ=255)
                                                                                                                                                                                            tGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGAGGCAGCTTTAACAGCCACTACGGTGTCCAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACCCGTCATTCGCTCTCCCGCTCTGTAATTGCCGAAGAAAATGGCGTTATGGAGCGTGATTTCGCCTACACCATTGGTCGGGCAAGAAGGGAGCTGCaaa                                                 >  2:204818/1‑239 (MQ=255)
                                                                                                                                                                                             gAACGGGGGCCCGGCGCAGAACAGGCGGCATTGCGGGCGGACAAACGCTTCACCAATGCCGAAGGTGGCCGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAAc                                                <  1:412942/239‑1 (MQ=255)
                                                                                                                                                                                                   gcggcCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGGTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCACCTGTGTAATTGCCGAAGAAAAGGGCGAGATGGAGCGTGATTACGCCTCACCCCTTGGGCGTGggaggacggagcgtccaacgccaaa                                          >  2:1193495/1‑220 (MQ=255)
                                                                                                                                                                                                     ggcCGGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCCCCAATTCCGAAGGTGGCAGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCagag                                         <  1:341168/238‑1 (MQ=255)
                                                                                                                                                                                                        ccGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATAAGGGATCTGCAAACGCCagatgggg                                     >  1:1015538/1‑234 (MQ=255)
                                                                                                                                                                                                        ccGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGGGTCAAAGTTTTTGGCAACAGCCACGGCCTGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAAAGGCGATATGGAGCGTGATTCAGCCTACACCATTGGCCGTGCGATGAGCGATCTCCAAACGCAAGAGgggg                                     >  2:33142/1‑239 (MQ=255)
                                                                                                                                                                                                        ccGCGCAAAACGGGTCGCATTGAAGGGGGAAAACGCTACACCCAATGCGGAGGGTGCCACTTTCAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGGTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTggg                                     <  2:326979/239‑1 (MQ=255)
                                                                                                                                                                                                                  acgggcggcATTGCAGGCGGACAACCCAATCACCAATACGGAGGGTGGCAGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCg                            <  1:1209340/236‑1 (MQ=255)
                                                                                                                                                                                                                        ggcATTGCGGCGGGAAAAACGCATCACAAATGCCGAGGGTGGCGCCTTTACCAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACtgt                       <  1:462736/238‑1 (MQ=255)
                                                                                                                                                                                                                           aTTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCCTTGGTCGTGCGATGAGCGATCTCCAAACGCCCGAGTGGGTTGGGGGCGACTGTGGTc                   >  1:165320/1‑239 (MQ=255)
                                                                                                                                                                                                                           aTTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGCCttggct                    >  1:411793/1‑233 (MQ=255)
                                                                                                                                                                                                                           aTTGCAGGCGAACAAACGCATCCCCAATACCGAAGGTGGCAGCTTTACCAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCt                    <  1:12459/238‑1 (MQ=255)
                                                                                                                                                                                                                             tGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGGGTAATTGCCGAAGAAAATGGCGATATGGAGCGGGATTACGCCTACACCATTGGTCGTGGGATGAGCGATATCGAAACGCCAGAGTTGGTTGGGGCCGACTGTGCTccc                 >  2:516108/1‑237 (MQ=255)
                                                                                                                                                                                                                               tggcgggcgcaaaaCGCACCCCCAATTGCGGAGGTGGGCAGCTTTACCAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCg                  <  2:1015538/228‑1 (MQ=255)
                                                                                                                                                                                                                               cggcgggcgcaaaacccaaccACCAATCCGGAGGGGGGCATCTTAACCAGCCACTACGGTTTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCg                  <  2:879439/219‑1 (MQ=255)
                                                                                                                                                                                                                                      ggACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGACGTGCGATGGGCGATCTGCAAACGCCAGAGCGGGTGGGGGCCGATTGGTCTCGCCGGACtt          >  2:967520/1‑239 (MQ=255)
                                                                                                                                                                                                                                          aaaCGCATCCCAAAGCCCGAAGGTGGCAGTTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATc      <  1:1096891/239‑1 (MQ=255)
                                                                                                                                                                                                                                              gtaTCCCCAAGACCGAAGGTGACAGCTTTAACAGCCCGTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGt  <  1:1104115/237‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
TAGCGCATCATCCACCGATTTAAGCCCGCAGGCCATTGCCCGTACTGTACAGGCGGCGCTGGATATTGCCCGTTATACCTCGCCAGATCCCTGTGCCGGCGTGGCAGACAAAGAGCTGCTGGCCTTTGACGCACCAGATCT‑‑CGACTTGTTCCACCCTGCGGAAGTTTCCCCGGATGAAGCCATTGAACTGGCGGCCCGCGCAGAACAGGCGGCATTGCAGGCGGACAAACGCATCACCAATACCGAAGGTGGCAGCTTTAACAGCCACTACGGTGTCAAAGTTTTTGGCAACAGCCACGGCATGTTGCAGGGTTACTGCTCAACGCGTCATTCGCTCTCCAGCTGTGTAATTGCCGAAGAAAATGGCGATATGGAGCGTGATTACGCCTACACCATTGGTCGTGCGATGAGCGATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGT  >  NC_000913/4458190‑4458658

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