Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 958,771 A→T 5.5% intergenic (+30/‑41) serC → / → aroA 3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913958,7710AT 5.5% 615.0 / 17.3 178intergenic (+30/‑41)serC/aroA3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base A (76/90);  new base T (4/6);  total (80/98)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01

Only 100 of 182 total aligned reads displayed.

TTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGCCGATCGTACGCGTTGCGAAATTA  >  NC_000913/958505‑959037
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                          
ttCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACagccagccg                                                                                                                                                                                                                                                                        >  2:190764/1‑270 (MQ=255)
           tGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                                         <  2:303865/275‑1 (MQ=255)
                tGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgtt                                                                                                                                                                                                                                                    <  2:466861/275‑1 (MQ=255)
                            aCCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAAGTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                                          >  2:536125/1‑273 (MQ=255)
                                                                      cgcTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCACTGACGTTAAAACCAAACGTTCGTGTCg                                                                                                                                                                                               >  4:61033/1‑274 (MQ=255)
                                                                         attaaaaaTTTTTCCTTGAAGGGTTCTTTGCGCCTGGCCTTTCTGCACTGCAAGGTCACCGTGTGGTCGGCGGCATGCGCGCTTCTATATATAACGCCAGGCCGCTGGAAGGCGTCAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                                                                                                                                                           <  4:100369/270‑1 (MQ=255)
                                                                          tGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                                                                                                                                                           <  2:61034/274‑1 (MQ=255)
                                                                                      ccTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACCGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgt                                                                                                                                                                                                                                                     >  4:26147/1‑204 (MQ=255)
                                                                                          gAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTa                                                                                                                                                                                                                                                 <  2:26147/204‑1 (MQ=255)
                                                                                               gTCTTTTGCTGCGGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTcc                                                                                                                                                                     <  2:374967/275‑1 (MQ=255)
                                                                                                tCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                                          >  2:633253/1‑205 (MQ=255)
                                                                                                    ttGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttgagt                                                                                                                                                                                                                                      <  4:633251/205‑1 (MQ=255)
                                                                                                          ctgGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgt                                                                                                                                                                                                                                                     >  1:116020/1‑184 (MQ=255)
                                                                                                          ctgGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATc                                                                                                                                                                                                                            >  2:166111/1‑209 (MQ=255)
                                                                                                              ccTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTa                                                                                                                                                                                                                                                 <  3:116019/184‑1 (MQ=255)
                                                                                                              ccTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTg                                                                                                                                                                                                                        <  4:166110/209‑1 (MQ=255)
                                                                                                                ttCATGCCGTGAAAGGTCATCGTGTTGGCCGCGCAATGCGCGGTTCTGTTTATAACGGCAGGCCGTTGTAAGTCGGTAAAACGCTGACAGACTTCCTGGTTTAGTTCGTACGCCGTTAATGGTAATGCCGAAATTTTGCTTAATCCCCACAGCCGGCCGGTGGGGGTTTGATTTTTGTAGTAGAGAAGTGGGTTCAAGGAATCCCTGCCCTTACAACCCGTCGCTCGTGTCGATGGTACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcg                                                                                                                                                    <  4:536123/275‑1 (MQ=255)
                                                                                                                ttCATGCCCTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcg                                                                                                                                                    <  2:647168/275‑1 (MQ=255)
                                                                                                                      cACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTaga                                                                                                                                                                                                                                               <  4:816366/178‑1 (MQ=255)
                                                                                                                          gAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACagccagcc                                                                                                                                                                                                                                                                         >  6:352883/1‑148 (MQ=255)
                                                                                                                           aaaGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtg                                                                                                                                                                                                                                                                       >  8:488180/1‑149 (MQ=255)
                                                                                                                           aaaGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATTTCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTg                                                                                                                                          >  2:235270/1‑274 (MQ=255)
                                                                                                                              ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                                                                                                                                     <  6:180478/148‑1 (MQ=255)
                                                                                                                              ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                                                                                                                                     <  5:912214/148‑1 (MQ=255)
                                                                                                                                 cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                                                                                                                                                  <  5:216390/148‑1 (MQ=255)
                                                                                                                                 cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                                                                                                                                                  <  6:39683/148‑1 (MQ=255)
                                                                                                                                 cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgggg                                                                                                                                                                                                                                                                  >  5:625017/1‑148 (MQ=255)
                                                                                                                                 cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgggg                                                                                                                                                                                                                                                                  >  6:754129/1‑148 (MQ=255)
                                                                                                                                  atcGTGTGGTCGGCGGAATGCGCGCTTCTTTTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                                                                                                                                                >  7:775287/3‑149 (MQ=255)
                                                                                                                                  aCCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                                                                                                                                                <  7:467135/149‑1 (MQ=255)
                                                                                                                                  aCCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                                                                                                                                                <  7:268020/149‑1 (MQ=255)
                                                                                                                                  aCCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGt                                                                                                                                                                                                                                                                 <  5:309099/148‑1 (MQ=255)
                                                                                                                                     gtgtGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCg                                                                                                                                                                                               >  4:353221/1‑211 (MQ=255)
                                                                                                                                       gtgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                                                                                                                                                            <  6:817913/148‑1 (MQ=255)
                                                                                                                                         ggTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                                                                                                                                                           <  2:353223/211‑1 (MQ=255)
                                                                                                                                            cggcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                                                      <  7:441074/149‑1 (MQ=255)
                                                                                                                                            cggcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCt                                                                                                                                                                                                                                                       >  6:630081/1‑148 (MQ=255)
                                                                                                                                                 gAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                                                                                                                                      <  5:832785/128‑1 (MQ=255)
                                                                                                                                                    tGCGCGCCTCTATATATAACGCCAGGCCGCTGGGAGGCGTTAAAGCGATGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTGTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCaaaa                                                                                                                 <  4:330820/274‑1 (MQ=255)
                                                                                                                                                     gcgcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                                                                                                                                                             <  7:536170/149‑1 (MQ=255)
                                                                                                                                                       gcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgttg                                                                                                                                                                                                                                                   >  8:169071/1‑141 (MQ=255)
                                                                                                                                                       gcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagagag                                                                                                                                                                                                                                            <  5:384958/148‑1 (MQ=255)
                                                                                                                                                        cgcTTCTATTTATAACGCCATGCCTCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagt                                                                                                                                                                                                                                           >  6:50542/1‑148 (MQ=255)
                                                                                                                                                        cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                                          <  7:65349/149‑1 (MQ=255)
                                                                                                                                                          cTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                                                                                                                                                        <  7:788606/149‑1 (MQ=255)
                                                                                                                                                            tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttgagt                                                                                                                                                                                                                                      <  8:471031/149‑1 (MQ=255)
                                                                                                                                                            tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttgag                                                                                                                                                                                                                                       >  6:641360/1‑148 (MQ=255)
                                                                                                                                                            tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTATAGAgagttgag                                                                                                                                                                                                                                       >  5:85759/1‑148 (MQ=255)
                                                                                                                                                              tATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTc                                                                                                                                                                                                                                    >  7:825048/1‑149 (MQ=255)
                                                                                                                                                                tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggga                                                                                                                                                                                                                                                                 >  8:23155/1‑117 (MQ=255)
                                                                                                                                                                tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggga                                                                                                                                                                                                                                                                 <  7:23155/118‑2 (MQ=255)
                                                                                                                                                                tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                                                                                                                                                  <  7:488180/149‑1 (MQ=255)
                                                                                                                                                                tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                                                                                                                                                  <  7:803136/149‑1 (MQ=255)
                                                                                                                                                                 ggaccggcccACGCCGCGGGAATGCGTTAAATCGCTGAAAGACTCAATGGTGGGGTTGGAAGGCGTTAACGGTTAATGCCGAAACTTGGCTTAATCCCCACAACCATGCGTTGTGGTTTTGATTTCTGTTGTGGGGATTTGACTTCGTCGCATCCCGGACGTTACAACCCATCGCTCATGTCGATGGCACTATTACTCTGACCGGTTCCAAAGGCGGTTCTAACCGCGCTTTATTGCTGACGGGATTATCACACGGCAAAACAGTATTAACCAAt                                                                                                   <  4:217508/267‑1 (MQ=255)
                                                                                                                                                                     aaCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACCACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCCGCCTTACCCAACGGCAAAACAGTATTAAACAAActg                                                                                                >  4:712147/1‑274 (MQ=255)
                                                                                                                                                                       cGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTTTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGc                                                                                                                                         >  1:196149/1‑231 (MQ=255)
                                                                                                                                                                         ccATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGa                                                                                                                                                                                                                               <  2:807069/143‑1 (MQ=255)
                                                                                                                                                                           aTGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCt                                                                                                                                                                                                                         <  2:999/147‑1 (MQ=255)
                                                                                                                                                                              ccGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGt                                                                                                                                                                                                                    >  8:203832/1‑149 (MQ=255)
                                                                                                                                                                                 cgggAAGGCGGGAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCAAGCCTGTGGGGttt                                                                                                                                                                                                                                                               >  4:390843/3‑103 (MQ=255)
                                                                                                                                                                                 cTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCTGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGGATTAACCAATCTGCTGGATAGCGAt                                                                                    >  2:42223/1‑274 (MQ=255)
                                                                                                                                                                                  tGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACTCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGctgtc                                                                                   >  2:217509/1‑270 (MQ=255)
                                                                                                                                                                                   ggAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTATAGAGAGTTGAGTTCATGGAATCCCTGACGTTACa                                                                                                                                                                                                                <  6:85759/148‑1 (MQ=255)
                                                                                                                                                                                     aaGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAAc                                                                                                                                                                                                              >  6:669877/1‑148 (MQ=255)
                                                                                                                                                                                           ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgt                                                                                                                                                                                                                                                     >  6:158504/1‑103 (MQ=255)
                                                                                                                                                                                           ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCg                                                                                                                                                                                                        >  5:327538/1‑148 (MQ=255)
                                                                                                                                                                                           ttAAAGCGCCGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgt                                                                                                                                                                                                                                                     <  5:158504/103‑1 (MQ=255)
                                                                                                                                                                                                gcgcTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                                                                                                                                                            <  7:172337/91‑1 (MQ=255)
                                                                                                                                                                                                    tGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCg                                                                                                                                                                                               <  5:101022/148‑1 (MQ=255)
                                                                                                                                                                                                      agagagttgtgggTGGAGTGCGAACGCCGTTAGGGTTAATGCCGAAATTTTGCTGAATCCCCACAGCCAGGCTGTGGGGTGTTTAGTTCTGTTGGAGGGTGTTGAGTTCATGGCATCCCTGACGTTACAACCCATCGCGCGTGTCGAGTGCACTATTAATCGGCCAGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGATGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCGTATGCTGaa                                                                <  4:235269/262‑1 (MQ=255)
                                                                                                                                                                                                      aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGAt                                                                                                                                                                                             <  5:669877/148‑1 (MQ=255)
                                                                                                                                                                                                      aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGAt                                                                                                                                                                                             <  6:625017/148‑1 (MQ=255)
                                                                                                                                                                                                      aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATg                                                                                                                                                                                            <  8:825048/149‑1 (MQ=255)
                                                                                                                                                                                                         gACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGc                                                                                                                                                                                          <  6:255278/148‑1 (MQ=255)
                                                                                                                                                                                                            ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACt                                                                                                                                                                                       <  6:612896/148‑1 (MQ=255)
                                                                                                                                                                                                            ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTa                                                                                                                                                                                      >  8:657016/1‑149 (MQ=255)
                                                                                                                                                                                                                tGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGt                                                                               <  4:212256/248‑1 (MQ=255)
                                                                                                                                                                                                                  gTTGAGTTCGAAGGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAACCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAAGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAGCCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGt                                                  <  1:1220/275‑1 (MQ=255)
                                                                                                                                                                                                                   tttaGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGtt                                                 <  2:55601/272‑1 (MQ=255)
                                                                                                                                                                                                                      aGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGcc                                                                                                                                                                            >  7:677738/1‑149 (MQ=255)
                                                                                                                                                                                                                      aGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCAATATGCTGAATGCATTAACAGCGTTAgg                                              >  1:153327/1‑275 (MQ=255)
                                                                                                                                                                                                                          cGAAAGCCGGAACGGTTAATGCCGAAATTTTGATTAATCCCCGCAGCCGGCCAGGTGGGTTTTTATTTCTGTTGTAGAAGGTTGGGTTAAGTGAATCACTGCCGTTACAACCCAGCGCGCGTGGCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAGCCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCGGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGt                                           <  4:42222/274‑1 (MQ=255)
                                                                                                                                                                                                                             aCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTcc                                                                                                                                                                     >  8:177362/1‑149 (MQ=255)
                                                                                                                                                                                                                               gCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCaa                                                                                                                                                                   <  8:313541/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAaga                                                                                                                                                                 <  8:726371/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAaga                                                                                                                                                                 <  8:412622/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAag                                                                                                                                                                  <  5:754129/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGc                                                                                                                                                               >  5:738333/1‑148 (MQ=255)
                                                                                                                                                                                                                                           aaTGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTa                                                                                                                                                        <  5:589438/148‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGAGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGAAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGCCg                       >  2:160511/1‑275 (MQ=255)
                                                                                                                                                                                                                                               ccGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGa                                                                                                                                                                                                                       <  3:189442/81‑1 (MQ=255)
                                                                                                                                                                                                                                                cGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcgc                                                                                                                                                   >  6:845140/1‑148 (MQ=255)
                                                                                                                                                                                                                                                 gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAGTCTGCCCGGTTCCAAGAGCGTTTCTAACCgcgc                                                                                                                                                 >  8:794507/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttttGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAt                                                                                                                                                                                                                                                           <  2:208740/39‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaccccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                                                                                                                                                           >  7:192652/3‑94 (MQ=255)
                                                                                                                                                                                                                                                                  ccgcagccagcCTGTGGGGTTTTTATGTTTGTTTGAGAGGTTTGAGTTCATGGAATCCCTTACTTTACAACCCATCGCTCGTGTCGCTGGCACTATTACTCTGCCCTGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGCCGATCGTACGCGTTGCGAAATTa  <  4:711996/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCa                                                                                                                                 >  6:618488/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt                                                                                                                                >  8:696206/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    acaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTa                                                                                                                              <  8:244211/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    acaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTa                                                                                                                              <  8:597094/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                          
TTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGCCGATCGTACGCGTTGCGAAATTA  >  NC_000913/958505‑959037

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