Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 958,773 C→G 5.7% 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.7% 422.7 / 7.9 119intergenic (+32/‑39)serC/aroA3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base C (50/62);  new base G (3/4);  total (53/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 122 total aligned reads displayed.

CGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGC  >  NC_000913/958501‑959014
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCt                                                                                                                                                                                                                                                 >  4:156709/1‑275 (MQ=255)
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGACGAAATTTTGCTTAATCCCCACAGCCAGCCt                                                                                                                                                                                                                                                 >  1:75800/1‑275 (MQ=255)
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATGGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGACTGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACGGACTTCCTGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGGTTAATCCCACACgccagccg                                                                                                                                                                                                                                                 >  4:295660/1‑274 (MQ=255)
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGGTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGGCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCCCCGCCAGCCt                                                                                                                                                                                                                                                 >  2:152921/1‑275 (MQ=255)
  aTTTCTACCTCAATGCCGTGGCGAAAGCTAACCATTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCTATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                                                                                                               <  2:9157/275‑1 (MQ=255)
     tCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTggg                                                                                                                                                                                                                                            >  4:412520/1‑275 (MQ=255)
               tGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                  <  2:406617/275‑1 (MQ=255)
               tGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                  <  1:93576/275‑1 (MQ=255)
                gACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                  <  2:199186/274‑1 (MQ=255)
                  cggggCGAAAGCTAACCGTTCGCGGATGTACGTGCCGTGCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGGCTTTTGCTGCTGGCCGTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTGATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAGCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  4:493966/272‑1 (MQ=255)
                        gAAAGCTAACCGTTGGTGGATGAACGTGCCGTTCAAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTGGAAGAGTCGTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTa                                                                                                                                                                                                                          <  2:295661/274‑1 (MQ=255)
                              tAACCGTTCGCGGATTAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                   <  4:431386/275‑1 (MQ=255)
                                      cgcgGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                                                                                                                           >  2:144861/1‑275 (MQ=255)
                                            tGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAAt                                                                                                                                                                                                      >  4:83019/1‑274 (MQ=255)
                                                                 cGGAAAGTGGGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTTGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCg                                                                                                                                                                                 <  4:280221/274‑1 (MQ=255)
                                                                             ttGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGTCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATg                                                                                                                                                                     >  2:224708/1‑274 (MQ=255)
                                                                             ttGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                                                                                                                                    <  2:83020/275‑1 (MQ=255)
                                                                                   aaTTGTTCCTTGAAGCGTCTTTTGCTGCTGGCCGTCATGCAATGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACGTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTa                                                                                                                                                               <  2:417579/274‑1 (MQ=255)
                                                                                       gTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCACGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCTTGTCGATGGCACTATTCAt                                                                                                                                                          >  2:119573/1‑275 (MQ=255)
                                                                                           cttgtgggttgttttGTTGCTGGCCTTCATGCACTCAACGGTAACCGTGTGGTCGGCGGAATGCGCGTCTCATTTTATTACGCCAGGCGGCTGGAAGGCGTTAAAGCTCTGACAGGCTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCATCCTTTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGc                                                                                                                                                      <  4:145615/264‑1 (MQ=255)
                                                                                                              ctgGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCCCAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGATTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCCCTATTAATCTTCCCGTTTCCAagtgatgttc                                                                                                                                   >  4:1657/1‑267 (MQ=255)
                                                                                                              ctgGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCCCTATTAATCTGCCCGGTTCCAAGAGCGttta                                                                                                                                   >  4:372259/1‑274 (MQ=255)
                                                                                                                    ttCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcg                                                                                                                             <  4:144860/275‑1 (MQ=255)
                                                                                                                               aaaGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTg                                                                                                                   >  4:366412/1‑274 (MQ=255)
                                                                                                                                  ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCtgtg                                                                                                                                                                                                                                              <  6:86849/148‑1 (MQ=255)
                                                                                                                                      aCCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTCTTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggcggc                                                                                                           >  1:103640/1‑275 (MQ=255)
                                                                                                                                         gtgtGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt                                                                                                         >  2:212541/1‑274 (MQ=255)
                                                                                                                                              gTCGGGGGAATGAGCGCTTTTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGca                                                                                                    <  4:361377/274‑1 (MQ=255)
                                                                                                                                                  gcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGccaacg                                                                                                >  2:336479/1‑269 (MQ=255)
                                                                                                                                                      aaTGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCTCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCTATCGCTCGTGTCTATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCCTTAGCACACGGCaac                                                                                           >  2:302518/1‑274 (MQ=255)
                                                                                                                                                       aTGTGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCaaaa                                                                                          <  2:156710/275‑1 (MQ=255)
                                                                                                                                                        ttcgagcgtCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCGCAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCaaaa                                                                                          <  4:64086/266‑1 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                                                                                                              <  7:293441/122‑1 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                                                                                                              <  8:576615/122‑1 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                                                                                                              >  8:293441/1‑122 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                                                                                                              >  7:576615/1‑122 (MQ=255)
                                                                                                                                                             gcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                   <  6:78534/148‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGTCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCa                                                                              <  4:224707/274‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  5:664861/129‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               >  6:664861/1‑129 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCa                                                                              <  2:1657/274‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCa                                                                              <  4:63772/274‑1 (MQ=255)
                                                                                                                                                                     ttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                            <  7:934696/147‑1 (MQ=255)
                                                                                                                                                                         aaCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTGTTGTTGGCGGCATTAGCACACGGCAAAACAGCATTAACAAATatgt                                                                        >  1:112653/1‑274 (MQ=255)
                                                                                                                                                                         aaCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                                                                                                                >  4:480271/1‑235 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                                                                                                               >  8:37916/1‑107 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                                                                                                               <  7:37916/107‑1 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        <  6:976322/114‑1 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        <  6:320559/114‑1 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        >  5:976322/1‑114 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        >  5:320559/1‑114 (MQ=255)
                                                                                                                                                                           cGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGTCATTACCACCGTACAAAACAGACTTAACAActctt                                                                         >  4:40437/1‑271 (MQ=255)
                                                                                                                                                                             ccATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggcgg                                                                                                            <  2:480272/235‑1 (MQ=255)
                                                                                                                                                                                tGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAg                                                                 >  2:467254/1‑275 (MQ=255)
                                                                                                                                                                                         aaGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGt                                                        <  2:133046/275‑1 (MQ=255)
                                                                                                                                                                                          aGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCAACAGCCAGCCTGTGGGGTTTCTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATTACAGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGATTTATTGCTGGCGGCATTAGACCACGGCAAAACAGTATTACCCAATCTGATGGATAGCGATGAcagg                                                       >  4:348924/1‑272 (MQ=255)
                                                                                                                                                                                           ggCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgg                                                                                                                                                                                                                                             >  7:800532/1‑92 (MQ=255)
                                                                                                                                                                                               ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  5:166667/102‑1 (MQ=255)
                                                                                                                                                                                                  aaGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtg                                                                                                                                                                                                                                                >  6:517074/1‑82 (MQ=255)
                                                                                                                                                                                                       cTGAAAGAATTCATGGTTGAGTTCGAACGCCGGCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTg                                           <  4:212540/274‑1 (MQ=255)
                                                                                                                                                                                                        tGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTGTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACGATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGaa                                         <  2:372260/275‑1 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                            >  5:976642/1‑109 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                            <  6:976642/109‑1 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCg                                                                                                                                                                             >  5:284185/1‑140 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCg                                                                                                                                                                             <  6:284185/140‑1 (MQ=255)
                                                                                                                                                                                                                ttccttgTTGAGTTCGCACCCCCGTACGTTTAATGCCGAAATTTTGCCTAAGCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTaa                                 <  4:336478/269‑1 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                                                                                                                           >  8:916500/1‑72 (MQ=39)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                                                                                                                           <  7:916500/72‑1 (MQ=39)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTc                                                                                                                                                                                                                                 <  8:931988/82‑1 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTc                                                                                                                                                                                                                                 >  7:931988/1‑82 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               >  8:643273/1‑84 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  7:643273/84‑1 (MQ=255)
                                                                                                                                                                                                                      gTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGt                           <  1:120783/275‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGAACGCCGTCACGGTTAATGCCGAACTTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGATAGTTGAGTTCATGGAATCCCTGCCGTTACAACCCATCGCTCGTGTCGATGGCACTATNAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGATTTATTGCTGGCGGCATTATCACACGGCAAAACAGCAGTAACCAANCATATGGATAGCGATGATGATGGCCNCATGATGAATGCATTAACAGCggcagtgg                     >  3:98869/1‑268 (MQ=255)
                                                                                                                                                                                                                            ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGc                                                                                                                                                                   >  6:533194/1‑133 (MQ=255)
                                                                                                                                                                                                                            ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGc                                                                                                                                                                   <  5:533194/133‑1 (MQ=255)
                                                                                                                                                                                                                             tCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        <  6:1082439/63‑1 (MQ=255)
                                                                                                                                                                                                                             tCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTc                                                                                                                                                                                                                                 >  5:935459/1‑70 (MQ=255)
                                                                                                                                                                                                                             tCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGt                    <  2:366413/275‑1 (MQ=255)
                                                                                                                                                                                                                             accacctccttCACAGTTAGTGCCCGAATTTTTNTTACACCCAGCGGCCAGCCTTGTGTGTTTTTTTTTCTGTTGTANGGCGTTGAGTTCATTTAATCCCTGCCGTTACCACCCAGCTCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAACGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGt                    <  3:112653/265‑1 (MQ=255)
                                                                                                                                                                                                                                aaCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACCGCCATCCTTTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTTACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGt                           <  1:98869/265‑1 (MQ=255)
                                                                                                                                                                                                                                 tcGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGATGTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGc                <  2:537830/274‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                                                                                                                                     <  7:1072978/58‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGc                                                                                                                  >  3:34560/1‑173 (MQ=255)
                                                                                                                                                                                                                                         aCGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggc                                                                                                              <  1:34560/173‑1 (MQ=255)
                                                                                                                                                                                                                                          cGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAagag                                                                                                                                         >  7:661182/1‑145 (MQ=255)
                                                                                                                                                                                                                                          cGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGCAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAagag                                                                                                                                         <  8:661182/145‑1 (MQ=255)
                                                                                                                                                                                                                                               aaTGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGTGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGc  <  2:348925/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    cGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGAt                                                                                                                                                                      <  8:301493/106‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATct                                                                          >  4:439022/1‑191 (MQ=255)
                                                                                                                                                                                                                                                               ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCGGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATctgctg                                                                      <  2:439023/191‑1 (MQ=255)
                                                                                                                                                                                                                                                                tAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                                                                                                                >  6:512756/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 aaTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTgg                                                                                                               <  5:512756/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccccACAGGCTGGGTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                                                                                                                                 >  6:749851/1‑47 (MQ=255)
                                                                                                                                                                                                                                                                    ccccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                                                                                                                                 <  5:749851/47‑1 (MQ=25)
                                                                                                                                                                                                                                                                      ccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAt                                                                                                                     <  6:111207/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cacaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGtt                                                                                                                                     >  5:525970/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                       cacaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt                                                                                                         >  5:176103/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                       cacaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCATGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGtt                                                                                                                                     <  8:778794/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                               agccTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTa                                                                                                                                 <  7:657691/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
CGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGC  >  NC_000913/958501‑959014

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