Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913729,7100TG35.7% 150.6 / 20.5 98V43G (GTG→GGG) rhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (14/49);  new base G (34/1);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 5.77e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGCACCT  >  NC_000913/729461‑729959
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                             
atatatatAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCTCTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTg                                                                                                                                                                                                                                                              <  1:853587/251‑1 (MQ=255)
      ataAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAACCATGGTCAGGGTTCAGCCGTGGGAAGCCTCCGGGCCCCCACCCGCGCGGGCCGTTCGGGGTGCCCCGCCGGGGGGGCGTcc                                                                                                                                                                                                                                                       >  2:1159503/1‑251 (MQ=17)
      ataAATATTTTCATAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTcc                                                                                                                                                                                                                                                       <  2:683163/251‑1 (MQ=255)
          atatTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCCGGTTGGTGGGGCCCGGCGGGGGGACggcccgcc                                                                                                                                                                                                                                                    >  2:908934/1‑244 (MQ=18)
           tatTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGccc                                                                                                                                                                                                                                                  >  1:520216/1‑250 (MQ=255)
            atTTTCCTAAGTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCAt                                                                                                                                                                                                                                                <  2:348007/251‑1 (MQ=255)
              tttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCCGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACggcccgccctcc                                                                                                                                                                                                                                                >  1:363963/1‑240 (MQ=17)
              tttCCTAATTGTTCTTATCTGACAGATATCTCACGTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTGCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCTGCGTTCCCCGGCGGGGTAACGTCCg                                                                                                                                                                                                                                                      <  2:416321/244‑1 (MQ=255)
                tCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCg                                                                                                                                                                                                                                             <  1:944820/250‑1 (MQ=255)
                 ccTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCACCCGGGGTACGCATCGGGGCCCCCCCCCGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGTGACGTCCCGCCCTCCggg                                                                                                                                                                                                                                           >  1:849069/1‑250 (MQ=17)
                      ttGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCACTGTCCCGGGTTCAGCCGGGGGACGCCTCGGTGCCCCCCCCGGCGTGGCCTGTTTGGTGTGCCCCGCCGGGGTGACGTCCGGCCCACCGGTCAATc                                                                                                                                                                                                                                      >  1:928399/1‑251 (MQ=17)
                       tGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGGGCCCCCACCGGCGTGGCCTGCTCGGGGTGCCCCGGCGGGGGGGCGGCCGGCCCTCCGGTCGATcc                                                                                                                                                                                                                                     >  2:814617/1‑251 (MQ=17)
                        gTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATcc                                                                                                                                                                                                                                     <  1:588867/250‑1 (MQ=255)
                        gTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGACGGCCGGCCACCCGGTCacaccc                                                                                                                                                                                                                                    >  1:899663/1‑246 (MQ=255)
                           cTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCCGGCAATTCgctg                                                                                                                                                                                                                                 >  1:1096796/1‑251 (MQ=255)
                             tATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGTCCGGCCATCCGGTCAATCCgctgcc                                                                                                                                                                                                                               >  1:252232/1‑250 (MQ=255)
                                cTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACCTCCGGCCCTCCGGTCCAACCCCTGCTCgg                                                                                                                                                                                                                            >  1:820463/1‑251 (MQ=14)
                                  gACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGGACGCCTCGGTGCCCCCACCGGCGTGGCCTGTTCGGAGGGCCCCCGCGGGGGGGCGGCCCGCCCTCCCGTCAATCCgccgccccgtg                                                                                                                                                                                                                          >  2:608780/1‑242 (MQ=14)
                                       atatCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGTTGCAAAg                                                                                                                                                                                                                    <  2:850928/251‑1 (MQ=255)
                                        tatCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGGCCCGGCGGGGGGGCGTCCCGCCCCCCCGTTAATCCGCTGCCCGGGGCACAgg                                                                                                                                                                                                                   >  2:439851/1‑251 (MQ=14)
                                         atCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGGCCGGCCATCCGGGCAATCCGCTGCTCGGGGCAAAggg                                                                                                                                                                                                                  >  1:566536/1‑250 (MQ=255)
                                           ctcACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGGGCCCCGGCGGGGGGACGTCCGGCCCTCCGGTCAATCCGCCGCTCGGTGCAAAAGTcc                                                                                                                                                                                                                >  1:136714/1‑251 (MQ=14)
                                              aCTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTTAATCCGCTGCTCGGTGCAAAGGTCCTTc                                                                                                                                                                                                             <  1:1159503/251‑1 (MQ=255)
                                                ttAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCAACGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCCGCGGGGGGGCGGCCCGCCCTCCCGGCAATCCCCTGCTCGGTGCAAAAGtccctccc                                                                                                                                                                                                           >  1:547674/1‑251 (MQ=14)
                                                 tAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTcc                                                                                                                                                                                                            <  1:814617/249‑1 (MQ=255)
                                                    ggCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGGCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCCGCGGGGGCCGTTCGGGGTGCCCCGGCGGGGGGAAGGCCGGCCATCCGGGCAAACCCCCGCTCGGGGCAAAAGGCCTTCCCgggg                                                                                                                                                                                                       >  1:10638/1‑249 (MQ=14)
                                                        ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGGCGGGGGGACGTCCGGCCCTCCCGGCACTCCGCTGCCCGGGGCAAAGGTCCTTCCCGGTGGAAc                                                                                                                                                                                                   >  1:977037/1‑251 (MQ=14)
                                                          cTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGGCAATCCGCCGCTCGGGGCAAAGGGCCTTCCCCGTGAAAccc                                                                                                                                                                                                 >  1:142923/1‑250 (MQ=14)
                                                          cTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCAGCGGTGCCCCCCCCCGCGTGGCCTGTTCGGTGTGGCCCGGCGGGGTGATGGCCGGGCATCCCGTCCATCCCCTTCCCCGTGCCAAAGGCCCTTCCCGGGAAAccc                                                                                                                                                                                                 >  1:1068433/1‑250 (MQ=14)
                                                            tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCCTCCGGGCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGAc                                                                                                                                                                                               >  1:420212/1‑251 (MQ=255)
                                                            tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCg                                                                                                                                                                                                 <  1:232443/249‑1 (MQ=255)
                                                            tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTCCGGCCCTCCGGTCCATCCCCTGCTCGGTGCAAAGGGCCCTCCCCGTGAAACCGAc                                                                                                                                                                                               >  1:323960/1‑251 (MQ=14)
                                                             ataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGAc                                                                                                                                                                                               <  2:1096796/250‑1 (MQ=255)
                                                                aaTCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGGCCCGCCCACCCGGCAATCCCCTGCCCGGTGCAAAggggcctccccggggaaccgccaccc                                                                                                                                                                                           >  2:277517/1‑227 (MQ=14)
                                                                      tAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCCGCCGGGGTACGCCTCCGTGCCCCCCCCCGCGTGGCCTGTTCGGGGGGCCCCCGCGGGGGGACGTTCCGCCCCCCCGTCACTCCCCTGATCGGGGCAAAGGCCCGTCCCCGTGACAccccgaacgccccac                                                                                                                                                                                     >  2:25169/1‑238 (MQ=9)
                                                                          gTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccgg                                                                                                                                                                                 <  1:850505/251‑1 (MQ=255)
                                                                           tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCCGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGGCTGTTCGGTGTGCCCCGGCGGGGTGACGGCCGGCCCTCCGGTCAATCCGCCGCCCGGTGGCAAGGGCCTTCCCGGGGAAACCCACATCGCCCCgcccggc                                                                                                                                                                                >  1:1027303/1‑251 (MQ=14)
                                                                              cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggccc                                                                                                                                                                              <  1:1063391/250‑1 (MQ=255)
                                                                                   gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAATCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGTTCGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTg                                                                                                                                                                          <  2:994388/249‑1 (MQ=255)
                                                                                     cagcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGt                                                                                                                                                                      <  2:81567/251‑1 (MQ=255)
                                                                                       gcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc                                                                                                                                                                    <  1:317235/251‑1 (MQ=255)
                                                                                               aTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATtctctc                                                                                                                                                           <  1:908934/251‑1 (MQ=255)
                                                                                                 aaCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTccc                                                                                                                                                         <  1:439851/251‑1 (MQ=255)
                                                                                                      gATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCAcc                                                                                                                                                    <  2:1112539/251‑1 (MQ=255)
                                                                                                          aaaCAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCCGCCGGGGTACGCATCGGGGCCCCCCCCGGCGTGGCCTGTTCGGGGTGCCCCCGCGGGGGGACGGCCCGCCATCCCGGCAATCCGCGGCTCGGGGGGAAGGTCCCTCCCCGTGAAAACCACACCGCCCTGCACCGCCCGCCGCCCTTCATTTCTccccccccctacc                                                                                                                                               >  2:870343/1‑241 (MQ=9)
                                                                                                                ggATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAtt                                                                                                                                                                 <  2:1085514/229‑1 (MQ=255)
                                                                                                                ggATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGGCGGCCCGCCAACCCGGTAATCCCCCTCCCGGTGCCAAAGGCCCTCCCCGGGGAAACGCCATCCCCCTGCCCGGCCCCCCGCCGTTCCTTCTCCCCCGCACCTAcacccgtt                                                                                                                                          >  2:409711/1‑245 (MQ=9)
                                                                                                                ggATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAACATTGTCCAGGGTTCAGCCGGGGTACGCCTCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGGCCGGCCATCCGGGCAATCCGCTGCTCGGGGCAAAGGTCCTTCCCGGTGAAAACGACATCGCCCCGCCCCGGCCGCCGCCGtttctt                                                                                                                                                                 >  1:1085514/1‑225 (MQ=14)
                                                                                                                 gATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTTGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGGCAAACCGCCTCTCGGGGGCAAGAGCCTTCCCCGGGAAAACGAAATCGCCCCGCCCGGCCCCCCGCCCGTCATTCTCTTCCGCAcccaccaccgttt                                                                                                                                         >  1:209410/1‑240 (MQ=14)
                                                                                                                  aTGAGCAGATGGGCGGAAATCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTTTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGTTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCGCTCCCGCACCTACAGCAGTTAc                                                                                                                                        <  2:928399/251‑1 (MQ=255)
                                                                                                                  aTGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTAc                                                                                                                                        <  1:277517/251‑1 (MQ=255)
                                                                                                                      gCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGa                                                                                                                                    <  1:608780/251‑1 (MQ=255)
                                                                                                                      gCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGGGCAAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCCGCCCTTCATTCTCTCCCCCACCTACAGCAGGTACCggg                                                                                                                                    >  1:101398/1‑250 (MQ=14)
                                                                                                                        aGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGGCCCGCCCTCCGGGCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCCACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTTCCCCACCTACCGCAGTTACCGGAc                                                                                                                                   >  1:944981/1‑250 (MQ=14)
                                                                                                                        aGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTCCGGCCCTCCGGGCAATCCGCTGCTCGGGGCAAAGGTCCTTCCCGGTGAAACCCACACCGCCCCGCCCCGCCCGCTGCCGTTCATTCTCTCCCGCACCCACAGCAGTTAcccggc                                                                                                                                   >  1:21150/1‑246 (MQ=14)
                                                                                                                          aTGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaa                                                                                                                                <  2:566536/251‑1 (MQ=255)
                                                                                                                            gAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaaaa                                                                                                                              <  2:10638/251‑1 (MQ=255)
                                                                                                                                ggAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCTGtt                                                                                                                                          <  1:824550/235‑1 (MQ=255)
                                                                                                                                             cgcgGCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCAGGGTGCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGGCCGGCCATCCGGTCAAGCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTgggg                                                                                                             <  2:209410/251‑1 (MQ=255)
                                                                                                                                                 ccAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGcc                                                                                                         <  2:520216/250‑1 (MQ=40)
                                                                                                                                                   aGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGAAGGCCGGGCATCCGGGCAAAACGCTGCTCGGTGCCAAGGTCCTTCCCGGGGAAACCGACATCGCCCTTCCCGGCCCCCTGCCGGTCATTTCTTCCCGCACCCACAGCCGTTACCGGACAAAAACGCCCCCGCCGGGGGGGGGcccc                                                                                                       >  1:497326/1‑249 (MQ=14)
                                                                                                                                                    gggCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGACGTCCGGCCATCCCGTCAATCCCCTGCTCGGGGCAAAAGGCCTTCCCCGGGAAACCGACATCGCCCTGCCCCGCCCCCTGCCGTTTATTTTCTCCCCCACCCACAGCAGTTTCCGGGCAAAAAAGCCCGCGCCGGGGGGGGGccccc                                                                                                      >  1:552890/1‑248 (MQ=11)
                                                                                                                                                       cGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGGCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTTCCCGGCCCGCTGCCGTTCAGTCTCTCCCGCACCTACAGCAGTTACCGGATAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc                                                                                                   <  1:1034507/251‑1 (MQ=255)
                                                                                                                                                           aTGACGCAGTATGGCGGTAGCATTGTCCAGGGTGCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCTCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCgg                                                                                               <  2:820463/251‑1 (MQ=255)
                                                                                                                                                                     aTGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGGCCGGCCATCCGGGCAATTCCCTGCTCGGGGCAAAGGGCCTTTCCGGGGAAACCGACATCGCCCTGCCCGGCCCGCCGCCGTTCATTCTCTCCCGCACCCACAGCCGTTACCGGACCAAAACGCCCGCCCCGGTGGGGAGCCCCGGCCCCGGGCGGAAAATg                                                                                     >  1:391798/1‑251 (MQ=11)
                                                                                                                                                                      tGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAACTGc                                                                                    >  2:716730/1‑251 (MQ=255)
                                                                                                                                                                             aGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCgg                                                                              <  1:862546/250‑1 (MQ=40)
                                                                                                                                                                                   tGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATcc                                                                        >  1:833533/1‑251 (MQ=40)
                                                                                                                                                                                          ggTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACa                                                                 <  1:409711/251‑1 (MQ=38)
                                                                                                                                                                                              cAGCCGGGGAACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACCGCCGGCCCTCCGGTCAATCCGCTGCTCGGTGCAAAGGGCCTTCCCGGTTGAACCGACATCGCCCCGCCCGGCCCGCTGCCGTTCCTTCTCTCCCGCACCTAAAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGGGCCTCGGCCCCGACCGGAAAATGCCTGGGGAGATCCGCttcgcgctg                                                             >  2:1057430/1‑244 (MQ=14)
                                                                                                                                                                                                   ggggTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGCGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCCACCGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCAGCGGCAACCCGCTTACAGccccgcgg                                                        >  2:73365/1‑244 (MQ=14)
                                                                                                                                                                                                        ccGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAaca                                                   <  2:101398/250‑1 (MQ=37)
                                                                                                                                                                                                          gCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATg                                                                                     <  1:556463/215‑1 (MQ=38)
                                                                                                                                                                                                           cATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAg                                                                <  1:645387/235‑1 (MQ=37)
                                                                                                                                                                                                           cATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacacac                                                >  1:541919/1‑251 (MQ=37)
                                                                                                                                                                                                           cATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGGTCATTCTCTCCCGCACCTACAGCAGGTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGGTATCCGCTTAccg                                                                >  2:645387/1‑233 (MQ=255)
                                                                                                                                                                                                            aTCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTTCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTACCGCACCTACAGCAGTTACCGGGCAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGACTGCGGAAATCCGCTTACAGCTGCGCGATAACACACt                                               >  1:941540/1‑251 (MQ=255)
                                                                                                                                                                                                             ccgggggccccACCGGCGGGGCCTGTTCGGGGTGCCCCGGCGGGGTGGCGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTg                                              <  2:955148/244‑1 (MQ=255)
                                                                                                                                                                                                               ggTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGAt                                            <  1:716730/251‑1 (MQ=34)
                                                                                                                                                                                                                  gCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt                                         >  1:29499/1‑251 (MQ=34)
                                                                                                                                                                                                                  gCCCCCACCAGCGTGGACTGTTCGGTGTGCCCCGGCAGGGTGGCGTCCGGCCCTACGGTCAATCGGCTGCGCGGTGCAAAGGTCCTTCCCGGGGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGGTCATTCTCTACCGCATCTGCCGCAGCTGCCGGACAGAAAAGCACGCGGCGGTGGGGAGCCTCCGCCCCAGCGGGAAAATTCCTGCGGATATCCGCCTGCAGCGCCGCGCTAACCCACTGATaac                                         >  2:150525/1‑249 (MQ=2)
                                                                                                                                                                                                                        cccGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGa                                   <  2:21150/250‑1 (MQ=21)
                                                                                                                                                                                                                         cgggCGGGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGAt                                  <  2:363963/249‑1 (MQ=21)
                                                                                                                                                                                                                            gcgggggaggTTCGGTGGGCCCCGGCGGGGGGACGTCCGGCCAGCCGGTCAATCCGCTGCTCGGTGAAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATCCTCAGTGATAAc                               <  2:1068433/242‑1 (MQ=11)
                                                                                                                                                                                                                              gTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcgg                             <  1:1084240/251‑1 (MQ=34)
                                                                                                                                                                                                                                  ccTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGGTAACACACTGATACTCAGTGATCACggcggc                         >  2:445706/1‑251 (MQ=21)
                                                                                                                                                                                                                                    tGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCACAGGTCCTTCCCGGTGAAACCGACATCGCCCCGCCCGGCCCGCTGCCGTTCCTTCTCTCCCGCACCTACAGCAGTTACCGGACCAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAAAGCCTGCGGATATCCGCTTACAGCTGCGCGAGCACACACTGATGCTcccggataccagccgg                         >  2:31354/1‑234 (MQ=11)
                                                                                                                                                                                                                                       gcGGTTTGCCCCGGCGGGGTGACGGCCGGCCATCCGGTCAAGCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCATGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGTCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATTCTCAGTGATAACGGCTGCAGAAg                    <  2:399203/250‑1 (MQ=11)
                                                                                                                                                                                                                                        cGGTGTGCCCCGGCGGGGTGACGCCCGGCCATCCGGTCAAGCCGCTGCTCGGTGCAAAGTTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTTCCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCGGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGc                   <  1:1098229/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGcc                  <  1:870343/251‑1 (MQ=34)
                                                                                                                                                                                                                                              gCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGaa                     <  1:859899/243‑1 (MQ=35)
                                                                                                                                                                                                                                              gCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAACAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGGTAACACACTGATACTCAGTGATAACGGCGGGAGaaagatcggg             >  2:859899/1‑243 (MQ=11)
                                                                                                                                                                                                                                               ccccGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAtt            <  1:511159/251‑1 (MQ=34)
                                                                                                                                                                                                                                                cccGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCGCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACt                                               <  1:996928/215‑1 (MQ=25)
                                                                                                                                                                                                                                                 ccGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAtttt          <  1:425960/251‑1 (MQ=34)
                                                                                                                                                                                                                                                      ggggTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGc      >  1:1022203/1‑250 (MQ=34)
                                                                                                                                                                                                                                                       gggTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGCAc    <  2:1038546/251‑1 (MQ=33)
                                                                                                                                                                                                                                                         gTGACGGCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGCACCt  <  2:676560/251‑1 (MQ=18)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                             
ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGCACCT  >  NC_000913/729461‑729959

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