Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,428,569 T→C 11.5% noncoding (192/1542 nt) rrsD ← 16S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,428,5690TC11.5% 215.3 / 3.6 133noncoding (192/1542 nt)rrsD16S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (56/51);  new base C (16/2);  total (78/55)
Fisher's exact test for biased strand distribution p-value = 3.89e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

CCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACT  >  NC_000913/3428336‑3428788
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                               
cccGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGCGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCAAGGGTCCCCCCCTTTGGCC‑TTg                                                                                                                                                                                                                            >  2:661532/1‑238 (MQ=9)
 ccGAAGGCCTTCTTCATGCACGGGGCATGGCTGCATAAGGCTTGCGTCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑Cg                                                                                                                                                                                                                         <  2:205015/239‑1 (MQ=9)
 ccGAAGGCCTTCTTCATACCCGCGGCCTGGCTGCATCAGGCTTGCCCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑Cg                                                                                                                                                                                                                         <  2:455358/239‑1 (MQ=14)
     aGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGGCGCCTTGGTGGGCCGTTACCCCACCAACAAGCCAATCCCAACT‑GGGCACAACCCATGGCAAAAGGCCCGAGGGCCCCCCCCTTTGGCC‑TTG‑CGCCgg                                                                                                                                                                                                                     >  1:917091/1‑238 (MQ=11)
         cttcttCCTGCCCGCGGCATGGCTGCATCAGGCTTGCGCGCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATg                                                                                                                                                                                                                 <  2:474581/239‑1 (MQ=12)
         cttcttCATACCCGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATg                                                                                                                                                                                                                 <  2:188504/239‑1 (MQ=17)
         cttctgCATCCACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAAGATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGGGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATg                                                                                                                                                                                                                 <  2:117623/239‑1 (MQ=12)
          ttcttcAATCACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATg                                                                                                                                                                                                                 <  2:1021180/238‑1 (MQ=17)
          gcaggggatgtgtaaaagaggCCGCATCAGGCTTGCGCCCATTGTCCAATATTCCCCCCTGCTGCCTTCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGGGGCTGGTCATC‑CTCTTAGACCAGCTTGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATg                                                                                                                                                                                                                 <  2:1084947/219‑1 (MQ=11)
                       cgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCCCTTTGGTC‑TTG‑CGACGTTATGCGGGATTATCTAcc                                                                                                                                                                                                   >  2:215087/1‑239 (MQ=17)
                                    tCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCCCTTTGGCC‑TTG‑CGACGTTATGGGGGATTAGCTCCCGTTTCCAGTAGtt                                                                                                                                                                                      >  2:158980/1‑239 (MQ=12)
                                    tCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTAGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGAAGtt                                                                                                                                                                                      >  2:697961/1‑239 (MQ=17)
                                        gCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAAAGGCCCGCGGGCCCCCCCCTTTGGCT‑TTG‑CGACGTTTTGCGGGATTAGCTAACGCTTCCAAGAGTAATc                                                                                                                                                                                   >  1:1105082/1‑238 (MQ=11)
                                          ttgtctccctttttcAATTTTCCCCACTGCGTCCTCCCGTTGGCATCTGGGCCGGTGCTCAGGTCCAGCGTCGCTGGT‑ATCTCTCTCATCCCATTCAGGGATCTGCGCCTTGGTGAGGCGTTACCCCACCAACAAGCTAATCCCATCT‑GGTCTCATCCGATGGCAAGAGCCTCGAGGGTACCCCGCTTTGGTC‑TTG‑CGACGTTCTGTGGTATTAGCTACCGTTTCTACGATTTATcccc                                                                                                                                                                                <  1:206017/225‑1 (MQ=2)
                                          ttGCGCACATTGTGCAATATCCCCCACTGCCGCCTCCCGGAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc                                                                                                                                                                                <  2:176383/239‑1 (MQ=12)
                                          ttGCCCCCATTGTGCAATATGCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCAAACAACAAGCTAATCCAATTT‑GGGCACATACGATGGCAAGAGGCTCGAGGGTCCCCCTCTTTGGTC‑TTG‑AGAAGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc                                                                                                                                                                                <  2:568502/239‑1 (MQ=11)
                                            gcgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGTGGCCCGAGGGTCCCCCTCTTCGGTC‑TTC‑CGCCGTTATGCGGACTTAGCTACCGGTTCCagtagtaaatcccc                                                                                                                                                                               >  1:51513/1‑230 (MQ=11)
                                             cgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCt                                                                                                                                                                              >  2:403570/1‑238 (MQ=32)
                                             cgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGGGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGAAGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGCC‑TTG‑CGACGTATTGCGGGAGTACAAACCGGTTCCAGGTTTTATcccccc                                                                                                                                                                              >  1:931948/1‑237 (MQ=11)
                                                 ccttGTTCAATTTTCCCCCCTGCTGCCTCCCGTAGGAGGCCGGACCGTGTCTCAGTTCCAGTGGGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc                                                                                                                                                                         <  1:697961/237‑1 (MQ=255)
                                                  aTTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGCTGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑GTT‑CGCCGTTACTCGGCATTAGCATACGTTTCCAGTTGTTACACCcctccctcc                                                                                                                                                                        >  1:568502/1‑238 (MQ=11)
                                                  aTTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGTCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCa                                                                                                                                                                        >  2:975589/1‑239 (MQ=17)
                                                      tgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGGGTTGCGTGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGGGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑Ca                                                                                                                                                                   <  2:1224441/239‑1 (MQ=12)
                                                      tccagtctcccccGCTGCTCCGCCCTGTGGGAGTCTGCACCTTGTCTCCTTTCCAGTGTTGCTGGGCACC‑CTCTCAGACCAGCTGGGTATCGGCGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTGATCCCATCG‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTGTGGGT‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑Ca                                                                                                                                                                   <  1:983679/230‑1 (MQ=11)
                                                           tatTCCCTCCTGGTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc                                                                                                                                                                                <  2:89712/222‑1 (MQ=17)
                                                                cccccgtgtgctccccGTTGGGGTTGGGCCCGTGTGTCAGTTCTCGGGTGGCTGGTCACC‑CCCTCAGACCAGCTAGGGATCGTCGCCTTGGTCAGCCGTTACCCCTCCAACAAGGTAATCCCATCT‑GGGCAAATCCGATGGCAAGAGGCCCGAGGGTCCCCCGCTTTGGTC‑TTG‑CCACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGa                                                                                                                                                         <  1:1291763/226‑1 (MQ=11)
                                                                cccccgggtgcCCCCCGTAGGTGTCGTGCCTGGGTCCCGGTTCCAGTGTGGCTGGGCACC‑CTCTCAACCCAGCTTGGGAACGTCGCCTGGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCGCATCCGATGGCAAGAAGCCCGAGGGTCCCCCCCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCATTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGa                                                                                                                                                         <  1:1088046/231‑1 (MQ=11)
                                                                cccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTg                                                                                                                                                                                                                            <  2:686672/174‑1 (MQ=31)
                                                                cccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCCACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑AAGTTTCCCaa                                                                                                                                                          >  1:300828/1‑237 (MQ=18)
                                                                cccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCCAGGGTCCCCCTCTTTGGGC‑TTG‑CACCGATATGCGGAATTAACTACCGTTTCCCGACATTAACCCCCCCCCAAAGG‑CAGTTTCCCAGa                                                                                                                                                         >  1:639519/1‑239 (MQ=11)
                                                                cccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCCAGGGTCCCCCTCTTTGGCC‑TTC‑CGCAGTTATGCGGAATTAGCTACCGTTTCCCGAAGTAATCCCCCCCCCACAGG‑CATTTTCCCAg                                                                                                                                                          >  1:569698/1‑238 (MQ=11)
                                                                cccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGGGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTg                                                                                                                                                                                                                            >  1:686672/1‑174 (MQ=14)
                                                                     tgctgcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGCTGGCACGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTC‑CGACGTTATGCGGGATTGGGTACCGTTTCCCGTAGTTATCCCCCTcctataggcctgtctcttatacca                                                                                                                                                      >  1:89712/1‑216 (MQ=11)
                                                                       ctgcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGCCATTa                                                                                                                                                    >  2:804856/1‑237 (MQ=18)
                                                                        cgcCTCCCGTAGGGGTCCGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACt                                                                                                                                                  <  2:76223/237‑1 (MQ=18)
                                                                               cGTAGGAGTCTGGACCGGGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGGA‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCAAGTAGTTATCACCCCCCACCAGG‑CAGTTTGCCAGACATTAATCccccct                                                                                                                                           >  1:868929/1‑236 (MQ=11)
                                                                                        cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTc                                                                                                                                  >  2:295450/1‑238 (MQ=33)
                                                                                         tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCTTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCACACATTACTCACCCGTCTGCCACTCg                                                                                                                                 >  2:985243/1‑238 (MQ=18)
                                                                                         tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCATCAATTACTCACCCGTCCGCCACTCg                                                                                                                                 >  2:40540/1‑238 (MQ=18)
                                                                                         tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGCC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCCGTAGTTATccc                                                                                                                                                                                 >  2:1086388/1‑191 (MQ=18)
                                                                                               gtgtCTCAGTTCCAGTGGGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCACATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGCCCCCCTCTTTGGGC‑TTG‑CGACGTTATGCGGTATTAGACACCGTTTCCGGCAGTTATCCCCCCCCATCAGG‑ACGTTTCCCAGACATTACTCACCCCTCCCCCACCCgcccccaa                                                                                                                          >  1:1036193/1‑232 (MQ=11)
                                                                                                tgtCTCAGTTCCAGTGTTGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAACCCACCAACAAGCTAATCCCATCT‑GTGCACATCCGATGGCAAGAGGCCCGAGGGGCACCATCTTTGGTC‑TTG‑CGTCGCTATGCGGTATTGCCTACCGTTTCCCGTTGGTCTCCCCCTCCCACAGG‑CAGGTTCCCAGACATTACTCACCCCGACGCCACCCGT‑CAGCgca                                                                                                                        >  1:130036/1‑237 (MQ=11)
                                                                                                tgtCCCAGTCCCAGTGTATCTGGTCACC‑CTCTCAGCCCGGCTACGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAACTAATCCCATCT‑AGGCACATCCGATGGCAAGAGGCCCGAGGGTACCCCTCTACGGCC‑TTG‑CCACGTGCCGCGGTCTGAGCTACCGTTTCCCGTAGTTCTCCCCCTCC‑TCAGGTCAGTTCCCCAAACATTCCTCACCCACCCCCCAcccccaagggca                                                                                                                         >  1:1158522/1‑228 (MQ=2)
                                                                                                tgtCACAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCAGATGGCCAGAGGCCCGAGGTACCCCCTCTTTGGGC‑CTG‑CGCCGTTATGCGGGATTGGCTACCGTTTCCCGAAGACATACCCCTCCACCAGC‑ACGTTTCCCAGACATTACTCACCATTCCGCCACTCAT‑CTCCGaa                                                                                                                        >  1:966806/1‑239 (MQ=11)
                                                                                                             gtgtGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGCCGTTATGCGGGATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CCGTTTCCCAGCCATTACCCACCCGCCCGCCCCCCGT‑CACCGaaaaacaaaaccgtt                                                                                                           >  1:256433/1‑227 (MQ=11)
                                                                                                                  gCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CTGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CCCCGAACAAGCAAGCTGTTTcccg                                                                                                      >  2:658951/1‑237 (MQ=255)
                                                                                                                  gCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAACCTGTTTcccg                                                                                                      >  2:113514/1‑237 (MQ=37)
                                                                                                                  gCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CACCGAAACAGCAACCTGttt                                                                                                          >  2:1256806/1‑235 (MQ=37)
                                                                                                                  gCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGGTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGCCATTCCCCACCCGTCCCCCACCCGT‑CACCGAAACAGCAAACTGTTTCCTg                                                                                                      >  1:1261589/1‑239 (MQ=11)
                                                                                                                  gCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGAGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCCCATCCGATGGCAAGAGGCCCCAGGGTCCCCCTCTTTGGTC‑TTG‑AGCCGTTATGCGGTCTTACCTACCGGTACCTCTAGTTATCACCCTCCATCAGG‑CAGTTTACCAGCAATTACTCACCCGTCCGCCACTCGT‑AAGGGAAACAGCAAGCCCTTTCCta                                                                                                      >  1:1182673/1‑238 (MQ=11)
                                                                                                                  gCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTAGGTGAGACGTTACCCCACCTACTATCTAATCCCATCG‑AGGCACATCGGATGGCAAGAGGCCCGCAAGTCCCCCTCTATGGGC‑TGG‑CTAAGTTATGCGGGATTAGTGAACGTTGCCAGTAGTTATCACCCCCCGTCAGG‑CTGTTTCCCAGATCTTTCAAACCCGCACGCCACCCGA‑Cacaaaaaaatcagcagactcc                                                                                                         >  1:509602/1‑216 (MQ=0)
                                                                                                                        cATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCg                                                                                                                                        >  2:1078535/1‑200 (MQ=34)
                                                                                                                          tC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACg                                                                                                                                                                                                                      >  1:1215911/1‑121 (MQ=33)
                                                                                                                          tC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACg                                                                                                                                                                                                                      <  2:1215911/121‑1 (MQ=33)
                                                                                                                             ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGCCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGCTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCCCCCCT‑CAGCGAACCATCAAGCCG‑TTCCTGTTTcccagtg                                                                                            >  1:324384/1‑233 (MQ=11)
                                                                                                                                   acccgcttggggATTGTCGCTTGGGTTGGCAGTTCCCCCACCAATCAGCAGCTCCAATCG‑GGGCCTGGCGGAGGGCAAGAGGCCCGGGGTTCCCCCTCTTTGGTC‑TGG‑CGACGTTATGCGGTTTTAGCTCCCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCATGTTACCGTTCGACTTGc                                                                                      <  1:434587/230‑1 (MQ=11)
                                                                                                                                   accccgctgggcgtcgtctgtttGGTGGGCCGTTCCCCCCCCAACAAGCCAACCCCATCT‑GGCCACATCCGATGGCAAGAGGCCCGGGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTAGGCGGTATTAGCTACCGTTTCCAGTATTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTCCTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGc                                                                                      <  1:1256806/218‑1 (MQ=11)
                                                                                                                                   acccAGCTAGGGATCTTCGCCTTGGTGAGCCGTTCCCCAACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGc                                                                                      <  2:962573/237‑1 (MQ=255)
                                                                                                                                    aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCATTAGTTATCCCCCTCCATCAGG‑CAg                                                                                                                                                                  >  1:427260/1‑163 (MQ=18)
                                                                                                                                    aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAg                                                                                                                                                                  <  2:427260/163‑1 (MQ=33)
                                                                                                                                    aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACCCGT‑CACCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGCCTTCCa                                                                                     >  2:605237/1‑239 (MQ=255)
                                                                                                                                     ccAGCTAGGGATCGTCGCCTTGGTTGGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGc                                                                                      <  2:1158522/237‑1 (MQ=37)
                                                                                                                                     ccAGCTAGGGATCGTCGCCTTGGTGAGCCGTTCCCCCACCAACAAGCTAACCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGc                                                                                      <  2:300828/237‑1 (MQ=37)
                                                                                                                                         cTAGGGATCGTCGCCTTGGCGACCCGTTCCCCCACCAACAAGCTAACCCCATCT‑GGTCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTg                                                                                       <  2:1267971/232‑1 (MQ=255)
                                                                                                                                              gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGCCGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGC‑CAGCGAAAAAGCAAGCTGTTTCCTGTTACCGCTCCACTTGCATGTGTAAga                                                                            >  1:688916/1‑237 (MQ=255)
                                                                                                                                                  gtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTg                                                                        >  2:907066/1‑238 (MQ=37)
                                                                                                                                                   tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑AAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGCCCGgc                                                                       >  2:1180248/1‑238 (MQ=255)
                                                                                                                                                   tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCCCCCCTCGT‑CGTCGGAACAGCAAGCTGCTTCCTGTTACCGTTCGACTTCCATGTGTAcggccggc                                                                       >  1:831965/1‑238 (MQ=11)
                                                                                                                                                   tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGCC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCCCCACTCGTCCAGCAAAAC‑GCAAGCCGTTTCCTTTTACCGTTCGCCTTCCATTGGTTAAGccgcc                                                                       >  1:1300843/1‑235 (MQ=11)
                                                                                                                                                   tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGAGGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATACCCATCCCTCAGG‑CAGTTTCCCCACCATTACTAACCCGTCCGCCCCCCGC‑CACCGAAAACGCAAGACGTTTCCCGTTACCATTctccttcactgggttgggccaccc                                                                      >  1:670624/1‑216 (MQ=11)
                                                                                                                                                   tcgCCTTGGTGAGCCGTTACCCCACAAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGACCCGAGGGCCCCCCCCTTTGGTC‑TTG‑CGACGTTATGAGGTATTAGCTACAGTTTCCCGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCCGACATTACTCACACGTCGGCCACTCGT‑CCCCAGAACAGCAACGCGCTTCCTCTTACCACTCGagtagccagttaaaggcccgc                                                                       >  1:1188285/1‑218 (MQ=11)
                                                                                                                                                     ccgttggggggctggtcccccACCAACAAGCACATCCCATCT‑GGGCAAACCCGATGGCAAGGCGCCCGAGGGTCCCCCTCTATGGTC‑TTG‑CGACGTTATGCGGTATTAGTTGCCGTTTCCAGTAGTTATCCCCCTCCAGCAGG‑CAGTTTCCCAGAATTTACTCACCCTTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgccgc                                                                    <  1:1267892/222‑1 (MQ=11)
                                                                                                                                                             gAGCCGGTACCCCACCAAAAAGCTAACCCCATCT‑GGGCATATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTAGGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGtt                                                             <  2:487156/238‑1 (MQ=255)
                                                                                                                                                                ccGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CATCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGGGTTAGGCCTGCCGCCAGCGTTCaa                                                          >  2:468663/1‑238 (MQ=25)
                                                                                                                                                                ccGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCACCGTTCaa                                                          >  2:202670/1‑238 (MQ=25)
                                                                                                                                                                 cGTTCCCCCACCAACAAGCTAACCCCATTG‑GGGCACAGCGGATGGAAAGAGGCCCGAGGGTCCCCCTCTTTGGCC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGGAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAAt                                                         <  2:1182673/238‑1 (MQ=11)
                                                                                                                                                                  tttACCCCCCCACCAA‑CCCATCCCATCTGGGGCACAGCCGGTGGCAGGGGCCCCTGGGGTCCCCC‑CGTTGGTCTTTG‑CGACGTATGGCGGTATTAGCTCCCGTGTCCGGTAGTCACCCCCCTCCAACAGG‑CAGTTTCCCAGACATTCCTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCt                                                       <  1:975589/238‑1 (MQ=11)
                                                                                                                                                                      ccccccCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTCCTCCCCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGc                                                   <  2:565325/239‑1 (MQ=25)
                                                                                                                                                                      ccccACCAAAAAGCTAATCCCATTG‑GGGCACATCCGATGGCAAGAGGCCGGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTAGGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGc                                                   <  2:340663/239‑1 (MQ=255)
                                                                                                                                                                                aaGCTAACCCCATCT‑GGGCACTGCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCaaa                                         <  2:1317344/239‑1 (MQ=18)
                                                                                                                                                                                  gCTAATCCCATCT‑GGGCACACCCGATGGCAAGAGGCCCGAGCGTCCCCCGTTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTccc                                                                                                                                                            <  1:1060751/123‑1 (MQ=255)
                                                                                                                                                                                      aTCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCGGCCAGCGTTCAATCTGAGCCATGATCAAACCCTTc                                   >  2:891134/1‑239 (MQ=17)
                                                                                                                                                                                            tCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACt                                                                                                                                                  >  2:1060751/1‑123 (MQ=32)
                                                                                                                                                                                            tCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTa                             >  2:50010/1‑239 (MQ=32)
                                                                                                                                                                                            tCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CACCGAAACAGCAAGCTGTTTCCTGTTACCGTTCACCTGTCATGTGTAAGGCTTGCGCCCACCGTTCATACTAACCCAGAATCAACCCTTTCAAttaa                             >  1:535206/1‑237 (MQ=11)
                                                                                                                                                                                             cT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATTCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGCCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGGTACCGTTCGACTTGCATgtgt                                                                                >  2:843358/1‑187 (MQ=18)
                                                                                                                                                                                                 ggCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg                         >  2:449118/1‑239 (MQ=32)
                                                                                                                                                                                                    cctccggttGGCAAGAGGCCCGGGGTGCCCCCCTTTTGGTC‑TGG‑CGACGTTAGGCGGTTTTAGCTACCGTTTCCATTAGTTATCCCCCTCCATCAGG‑CAGTTCCCCGGACATTCCTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGttt                      <  1:1074279/231‑1 (MQ=11)
                                                                                                                                                                                                      ctCCGTTGGCACGTGCCCCTAGGTTCCCCTTCTTTGTTC‑TGC‑CTCACTTACGCGTTATTGCCTCCCTTTTCCAGCAGTTACCCCCCTCCAGCAGA‑CCGTCGCCCAGCCATGTCTCTCCCGTCCGCCCCTCGT‑CAGCGAAGCAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGa                    <  1:647029/238‑1 (MQ=2)
                                                                                                                                                                                                            ttggaaggcccccGAGCGTCCCCCT‑TTTGGTC‑TTGTCAGCGTTTAGCGGTATTCGCTTCCTTTTCCCCTTGTTACCCCCCCCCACCAGG‑AAGTTTCCCAGTCTTTACTCCCCCGTCCGCCACTCGT‑CGGCGAACCAGCAAGCTGTTTCTTGTTACCGTTCGACTTTCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCa              <  1:345406/229‑1 (MQ=11)
                                                                                                                                                                                                             gggaaGAGGCCCGGGGGCCCCCCTCTTTGGTC‑TTG‑CGACGTTAGGCGGTATTAGCTACCGTTTCCAGTATTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCa              <  2:716745/235‑1 (MQ=9)
                                                                                                                                                                                                             ggCAAGTGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCa              <  2:1134072/238‑1 (MQ=14)
                                                                                                                                                                                                             ggCAAGATGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCa              <  2:670624/238‑1 (MQ=14)
                                                                                                                                                                                                             ggCAAGAGCCCTGAGGGTCCCCCTCTTTGGTC‑GTG‑CGACGTTATGCGGTATTATCTTCCGTTTCCAGTAGTTACCCCCCTCCATCAGG‑CAGTTTCCCAGACATTCCTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCa              <  2:130036/238‑1 (MQ=255)
                                                                                                                                                                                                                aaGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTAACGTTCGACTTGCATGTGTTAGGCCTCCCCCAACCGTTCAATTGGAGCCATGATCAACCTCTTCACTTTAAAATTTTAATGCTCACAGa          >  1:573064/1‑239 (MQ=11)
                                                                                                                                                                                                                aaGAGGCCCGAGGGTCCCCCTCTTTGGCC‑TTG‑CGACGTTATGCGGAATTAGCTACCGCTTCCAGTCGTTATCCCCCTCAATCAGG‑CAGTTTCCCAGACATTCCTCAACCGTCCCCCACTCGT‑CAGCGAACCAGCAAGCTGTTTACTGTTCCCGTTCGACTTGCATGTGTTCACCCTGCCGCCAGCGTTCAATCGTGGCCATGTTAAAACTCTTCAATTTAAAAGTTTGATGCTAACAGa          >  1:797010/1‑239 (MQ=255)
                                                                                                                                                                                                                  gagGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCGGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGACAGTGTTAGGCCTGCCCCCCACGTTTAATCAGAGACATGATCAAACTCTTCACTTTAAAAATTTAGTTCTCAAAGAAt        >  1:1032113/1‑239 (MQ=11)
                                                                                                                                                                                                                      gccGAGGTTCCCCCTTTTTGTTC‑TTC‑CGACTTTATGCGTTATTACCTACCGTTTCCAATAGTTATCCCCCGCCAGTAGG‑CAATATCCCAGACATTACTCACCCGTCCCCCAAGCGT‑CAGCGAAACAGCAAGCTTTTTCCTGTTGCCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTaa     <  1:40540/237‑1 (MQ=11)
                                                                                                                                                                                                                        cGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACt  >  2:1009349/1‑239 (MQ=31)
                                                                                                                                                                                                                        cGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGGGTTAGCCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAAAATTAAACt  >  1:365888/1‑239 (MQ=14)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                               
CCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATC‑CTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT‑GGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTC‑TTG‑CGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGG‑CAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGT‑CAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACT  >  NC_000913/3428336‑3428788

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