Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 958,773 C→G 5.5% intergenic (+32/‑39) 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,7730CG 5.5% 616.8 / 19.2 178intergenic (+32/‑39)serC/aroA3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base C (76/91);  new base G (4/7);  total (80/98)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-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)
                                                                                                                           aaaGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtg                                                                                                                                                                                                                                                                       >  8:488180/1‑149 (MQ=255)
                                                                                                                           aaaGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATTTCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTg                                                                                                                                          >  2:235270/1‑274 (MQ=255)
                                                                                                                              ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                                                                                                                                     <  6:180478/148‑1 (MQ=255)
                                                                                                                                 cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                                                                                                                                                  <  6:639701/148‑1 (MQ=255)
                                                                                                                                 cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                                                                                                                                                  <  5:216390/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)
                                                                                                                                  aCCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGt                                                                                                                                                                                                                                                                 <  5:309099/148‑1 (MQ=255)
                                                                                                                                     gtgtGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGatta                                                                                                                                                                                                                                                              <  6:619824/148‑2 (MQ=255)
                                                                                                                                     gtgtGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCg                                                                                                                                                                                               >  4:353221/1‑211 (MQ=255)
                                                                                                                                         ggTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                                                                                                                                                           <  2:353223/211‑1 (MQ=255)
                                                                                                                                          gTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                                         <  6:291988/148‑1 (MQ=255)
                                                                                                                                                 gAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                                                                                                                                      <  5:832785/128‑1 (MQ=255)
                                                                                                                                                 gAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                                                                                                                                      >  6:832785/1‑128 (MQ=255)
                                                                                                                                                    tGCGCGCCTCTATATATAACGCCAGGCCGCTGGGAGGCGTTAAAGCGATGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTGTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCaaaa                                                                                                                 <  4:330820/274‑1 (MQ=255)
                                                                                                                                                     gcgcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                                                                                                                                                             <  7:536170/149‑1 (MQ=255)
                                                                                                                                                       gcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgttg                                                                                                                                                                                                                                                   >  8:169071/1‑141 (MQ=255)
                                                                                                                                                       gcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgttg                                                                                                                                                                                                                                                   <  7:169071/141‑1 (MQ=255)
                                                                                                                                                       gcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagagag                                                                                                                                                                                                                                            <  5:384958/148‑1 (MQ=255)
                                                                                                                                                        cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                                          <  7:65349/149‑1 (MQ=255)
                                                                                                                                                        cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagt                                                                                                                                                                                                                                           >  6:589438/1‑148 (MQ=255)
                                                                                                                                                          cTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                                                                                                                                                        <  7:788606/149‑1 (MQ=255)
                                                                                                                                                            tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttgagt                                                                                                                                                                                                                                      <  8:471031/149‑1 (MQ=255)
                                                                                                                                                            tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttgag                                                                                                                                                                                                                                       <  6:629300/148‑1 (MQ=255)
                                                                                                                                                            tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagt                                                                                                                                                                                                                                           >  6:762028/1‑144 (MQ=255)
                                                                                                                                                            tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTATAGAgagttgag                                                                                                                                                                                                                                       >  5:85759/1‑148 (MQ=255)
                                                                                                                                                             tccagtataACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGt                                                                                                                                                                                                                                                                 <  6:200410/116‑1 (MQ=255)
                                                                                                                                                              tATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTc                                                                                                                                                                                                                                    >  7:825048/1‑149 (MQ=255)
                                                                                                                                                                tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggga                                                                                                                                                                                                                                                                 >  8:23155/1‑117 (MQ=255)
                                                                                                                                                                tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggga                                                                                                                                                                                                                                                                 <  7:23155/118‑2 (MQ=255)
                                                                                                                                                                tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                                                   <  6:849241/148‑1 (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)
                                                                                                                                                                       cGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                                                   >  6:408847/1‑141 (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)
                                                                                                                                                                                           ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCg                                                                                                                                                                                                        >  5:327538/1‑148 (MQ=255)
                                                                                                                                                                                           ttAAAGCGCCGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgt                                                                                                                                                                                                                                                     <  5:158504/103‑1 (MQ=255)
                                                                                                                                                                                                    tGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCg                                                                                                                                                                                               <  5:101022/148‑1 (MQ=255)
                                                                                                                                                                                                      agagagttgtgggTGGAGTGCGAACGCCGTTAGGGTTAATGCCGAAATTTTGCTGAATCCCCACAGCCAGGCTGTGGGGTGTTTAGTTCTGTTGGAGGGTGTTGAGTTCATGGCATCCCTGACGTTACAACCCATCGCGCGTGTCGAGTGCACTATTAATCGGCCAGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGATGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCGTATGCTGaa                                                                <  4:235269/262‑1 (MQ=255)
                                                                                                                                                                                                      aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGAt                                                                                                                                                                                             <  5:669877/148‑1 (MQ=255)
                                                                                                                                                                                                      aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATg                                                                                                                                                                                            <  8:825048/149‑1 (MQ=255)
                                                                                                                                                                                                            ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTa                                                                                                                                                                                      >  8:657016/1‑149 (MQ=255)
                                                                                                                                                                                                            ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTa                                                                                                                                                                                      >  7:412622/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)
                                                                                                                                                                                                                                   tCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAagag                                                                                                                                                                >  6:135983/1‑148 (MQ=255)
                                                                                                                                                                                                                                    cACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGc                                                                                                                                                               >  5:738333/1‑148 (MQ=255)
                                                                                                                                                                                                                                      cGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGt                                                                                                                                                             <  6:65067/148‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTa                                                                                                                                                        <  5:589438/148‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGAGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGAAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGCCg                       >  2:160511/1‑275 (MQ=255)
                                                                                                                                                                                                                                               ccGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGa                                                                                                                                                                                                                       <  3:189442/81‑1 (MQ=255)
                                                                                                                                                                                                                                               ccGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcgc                                                                                                                                                   >  7:244211/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAGTCTGCCCGGTTCCAAGAGCGTTTCTAACCgcgc                                                                                                                                                 >  8:794507/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttttGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAt                                                                                                                                                                                                                                                           <  2:208740/39‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCa                                                                                                                                 >  7:444059/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 cccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCa                                                                                                                                 >  7:190713/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ccgcagccagcCTGTGGGGTTTTTATGTTTGTTTGAGAGGTTTGAGTTCATGGAATCCCTTACTTTACAACCCATCGCTCGTGTCGCTGGCACTATTACTCTGCCCTGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGCCGATCGTACGCGTTGCGAAATTa  <  4:711996/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCa                                                                                                                                 >  6:618488/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt                                                                                                                                >  8:696206/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   cacaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt                                                                                                                                >  6:386722/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    acaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTa                                                                                                                              <  8:597094/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    acaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTa                                                                                                                              <  8:244211/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 < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: