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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913958,7710AT 5.8% 413.8 / 10.3 119intergenic (+30/‑41)serC/aroA3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base A (50/62);  new base T (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 = 5.74e-01

Only 100 of 122 total aligned reads displayed.

CGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGC  >  NC_000913/958501‑959014
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCt                                                                                                                                                                                                                                                 >  4:156709/1‑275 (MQ=255)
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGACGAAATTTTGCTTAATCCCCACA‑GCCAGCCt                                                                                                                                                                                                                                                 >  1:75800/1‑275 (MQ=255)
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATGGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGACTGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACGGACTTCCTGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGGTTAAT‑CCCACACgccagccg                                                                                                                                                                                                                                                 >  4:295660/1‑274 (MQ=255)
cGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGGTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGGCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCCCC‑GCCAGCCt                                                                                                                                                                                                                                                 >  2:152921/1‑275 (MQ=255)
  aTTTCTACCTCAATGCCGTGGCGAAAGCTAACCATTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCTATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtgt                                                                                                                                                                                                                                               <  2:9157/275‑1 (MQ=255)
     tCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTggg                                                                                                                                                                                                                                            >  4:412520/1‑275 (MQ=255)
               tGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                  <  2:406617/275‑1 (MQ=255)
               tGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                  <  1:93576/275‑1 (MQ=255)
                gACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTAttt                                                                                                                                                                                                                                  <  2:199186/274‑1 (MQ=255)
                  cggggCGAAAGCTAACCGTTCGCGGATGTACGTGCCGTGCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGGCTTTTGCTGCTGGCCGTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTGATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAGCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  4:493966/272‑1 (MQ=255)
                        gAAAGCTAACCGTTGGTGGATGAACGTGCCGTTCAAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTGGAAGAGTCGTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTa                                                                                                                                                                                                                          <  2:295661/274‑1 (MQ=255)
                              tAACCGTTCGCGGATTAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                   <  4:431386/275‑1 (MQ=255)
                                      cgcgGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                                                                                                                           >  2:144861/1‑275 (MQ=255)
                                            tGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAAt                                                                                                                                                                                                      >  4:83019/1‑274 (MQ=255)
                                                                 cGGAAAGTGGGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTTGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCg                                                                                                                                                                                 <  4:280221/274‑1 (MQ=255)
                                                                             ttGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGTCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATg                                                                                                                                                                     >  2:224708/1‑274 (MQ=255)
                                                                             ttGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                                                                                                                                    <  2:83020/275‑1 (MQ=255)
                                                                                   aaTTGTTCCTTGAAGCGTCTTTTGCTGCTGGCCGTCATGCAATGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACGTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTa                                                                                                                                                               <  2:417579/274‑1 (MQ=255)
                                                                                       gTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCACGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCTTGTCGATGGCACTATTCAt                                                                                                                                                          >  2:119573/1‑275 (MQ=255)
                                                                                           cttgtgggttgttttGTTGCTGGCCTTCATGCACTCAACGGTAACCGTGTGGTCGGCGGAATGCGCGTCTCATTTTATTACGCCAGGCGGCTGGAAGGCGTTAAAGCTCTGACAGGCTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCATCCTTTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGc                                                                                                                                                      <  4:145615/264‑1 (MQ=255)
                                                                                                              ctgGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCCCA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGATTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCCCTATTAATCTTCCCGTTTCCAagtgatgttc                                                                                                                                   >  4:1657/1‑267 (MQ=255)
                                                                                                              ctgGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCCCTATTAATCTGCCCGGTTCCAAGAGCGttta                                                                                                                                   >  4:372259/1‑274 (MQ=255)
                                                                                                                    ttCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcg                                                                                                                             <  4:144860/275‑1 (MQ=255)
                                                                                                                               aaaGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTg                                                                                                                   >  4:366412/1‑274 (MQ=255)
                                                                                                                                  ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GGCTGGCtgtg                                                                                                                                                                                                                                              <  6:86849/148‑1 (MQ=255)
                                                                                                                                      aCCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTCTTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggcggc                                                                                                           >  1:103640/1‑275 (MQ=255)
                                                                                                                                         gtgtGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt                                                                                                         >  2:212541/1‑274 (MQ=255)
                                                                                                                                              gTCGGGGGAATGAGCGCTTTTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGca                                                                                                    <  4:361377/274‑1 (MQ=255)
                                                                                                                                                  gcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGccaacg                                                                                                >  2:336479/1‑269 (MQ=255)
                                                                                                                                                      aaTGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCTCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCTATCGCTCGTGTCTATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCCTTAGCACACGGCaac                                                                                           >  2:302518/1‑274 (MQ=255)
                                                                                                                                                       aTGTGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCaaaa                                                                                          <  2:156710/275‑1 (MQ=255)
                                                                                                                                                        ttcgagcgtCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCGCA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCaaaa                                                                                          <  4:64086/266‑1 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtgtg                                                                                                                                                                                                                                              <  7:293441/122‑1 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtgtg                                                                                                                                                                                                                                              <  8:576615/122‑1 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtgtg                                                                                                                                                                                                                                              >  8:293441/1‑122 (MQ=255)
                                                                                                                                                            cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtgtg                                                                                                                                                                                                                                              >  7:576615/1‑122 (MQ=255)
                                                                                                                                                             gcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                                                                                                                   <  6:78534/148‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGTCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCa                                                                              <  4:224707/274‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  5:664861/129‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               >  6:664861/1‑129 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCa                                                                              <  2:1657/274‑1 (MQ=255)
                                                                                                                                                                    tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCa                                                                              <  4:63772/274‑1 (MQ=255)
                                                                                                                                                                     ttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                            <  7:934696/147‑1 (MQ=255)
                                                                                                                                                                         aaCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTGTTGTTGGCGGCATTAGCACACGGCAAAACAGCATTAACAAATatgt                                                                        >  1:112653/1‑274 (MQ=255)
                                                                                                                                                                         aaCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                                                                                                                >  4:480271/1‑235 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtgt                                                                                                                                                                                                                                               >  8:37916/1‑107 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtgt                                                                                                                                                                                                                                               <  7:37916/107‑1 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        <  6:976322/114‑1 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        <  6:320559/114‑1 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        >  5:976322/1‑114 (MQ=255)
                                                                                                                                                                          aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        >  5:320559/1‑114 (MQ=255)
                                                                                                                                                                           cGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGTCATTACCACCGTACAAAACAGACTTAACAActctt                                                                         >  4:40437/1‑271 (MQ=255)
                                                                                                                                                                             ccATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggcgg                                                                                                            <  2:480272/235‑1 (MQ=255)
                                                                                                                                                                                tGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAg                                                                 >  2:467254/1‑275 (MQ=255)
                                                                                                                                                                                         aaGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGt                                                        <  2:133046/275‑1 (MQ=255)
                                                                                                                                                                                          aGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCAACA‑GCCAGCCTGTGGGGTTTCTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATTACAGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGATTTATTGCTGGCGGCATTAGACCACGGCAAAACAGTATTACCCAATCTGATGGATAGCGATGAcagg                                                       >  4:348924/1‑272 (MQ=255)
                                                                                                                                                                                           ggCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GGCTGGCTGTgg                                                                                                                                                                                                                                             >  7:800532/1‑92 (MQ=255)
                                                                                                                                                                                               ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  5:166667/102‑1 (MQ=255)
                                                                                                                                                                                                  aaGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCtg                                                                                                                                                                                                                                                >  6:517074/1‑82 (MQ=255)
                                                                                                                                                                                                       cTGAAAGAATTCATGGTTGAGTTCGAACGCCGGCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTg                                           <  4:212540/274‑1 (MQ=255)
                                                                                                                                                                                                        tGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTGTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACGATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGaa                                         <  2:372260/275‑1 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                            >  5:976642/1‑109 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCa                                                                                                                                                                                                            <  6:976642/109‑1 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCg                                                                                                                                                                             >  5:284185/1‑140 (MQ=255)
                                                                                                                                                                                                           cAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCg                                                                                                                                                                             <  6:284185/140‑1 (MQ=255)
                                                                                                                                                                                                                ttccttgTTGAGTTCGCACCCCCGTACGTTTAATGCCGAAATTTTGCCTAAGCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTaa                                 <  4:336478/269‑1 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GGCTGGCTGTgggg                                                                                                                                                                                                                                           >  8:916500/1‑72 (MQ=39)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GGCTGGCTGTgggg                                                                                                                                                                                                                                           <  7:916500/72‑1 (MQ=39)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTc                                                                                                                                                                                                                                 <  8:931988/82‑1 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTc                                                                                                                                                                                                                                 >  7:931988/1‑82 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               >  8:643273/1‑84 (MQ=255)
                                                                                                                                                                                                                 tCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                                                                                                               <  7:643273/84‑1 (MQ=255)
                                                                                                                                                                                                                      gTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGt                           <  1:120783/275‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGAACGCCGTCACGGTTAATGCCGAACTTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGATAGTTGAGTTCATGGAATCCCTGCCGTTACAACCCATCGCTCGTGTCGATGGCACTATNAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGATTTATTGCTGGCGGCATTATCACACGGCAAAACAGCAGTAACCAANCATATGGATAGCGATGATGATGGCCNCATGATGAATGCATTAACAGCggcagtgg                     >  3:98869/1‑268 (MQ=255)
                                                                                                                                                                                                                            ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGc                                                                                                                                                                   >  6:533194/1‑133 (MQ=255)
                                                                                                                                                                                                                            ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGc                                                                                                                                                                   <  5:533194/133‑1 (MQ=255)
                                                                                                                                                                                                                             tCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGttt                                                                                                                                                                                                                                        <  6:1082439/63‑1 (MQ=255)
                                                                                                                                                                                                                             tCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTc                                                                                                                                                                                                                                 >  5:935459/1‑70 (MQ=255)
                                                                                                                                                                                                                             tCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGt                    <  2:366413/275‑1 (MQ=255)
                                                                                                                                                                                                                             accacctccttCACAGTTAGTGCCCGAATTTTTNTTACACCCAGCG‑GCCAGCCTTGTGTGTTTTTTTTTCTGTTGTANGGCGTTGAGTTCATTTAATCCCTGCCGTTACCACCCAGCTCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAACGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGt                    <  3:112653/265‑1 (MQ=255)
                                                                                                                                                                                                                                aaCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACC‑GCCATCCTTTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTTACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGt                           <  1:98869/265‑1 (MQ=255)
                                                                                                                                                                                                                                 tcGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGATGTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGc                <  2:537830/274‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTa                                                                                                                                                                                                                                     <  7:1072978/58‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGc                                                                                                                  >  3:34560/1‑173 (MQ=255)
                                                                                                                                                                                                                                         aCGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggc                                                                                                              <  1:34560/173‑1 (MQ=255)
                                                                                                                                                                                                                                          cGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAagag                                                                                                                                         >  7:661182/1‑145 (MQ=255)
                                                                                                                                                                                                                                          cGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGCAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAagag                                                                                                                                         <  8:661182/145‑1 (MQ=255)
                                                                                                                                                                                                                                               aaTGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGTGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGc  <  2:348925/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    cGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGAt                                                                                                                                                                      <  8:301493/106‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATct                                                                          >  4:439022/1‑191 (MQ=255)
                                                                                                                                                                                                                                                               ttAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCGGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATctgctg                                                                      <  2:439023/191‑1 (MQ=255)
                                                                                                                                                                                                                                                                tAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                                                                                                                >  6:512756/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 aaTCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTgg                                                                                                               <  5:512756/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccccACA‑GGCTGGGTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                                                                                                                                 >  6:749851/1‑47 (MQ=255)
                                                                                                                                                                                                                                                                    ccccACA‑GGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                                                                                                                                 <  5:749851/47‑1 (MQ=25)
                                                                                                                                                                                                                                                                      ccACA‑GGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAt                                                                                                                     <  6:111207/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                       caca‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGtt                                                                                                                                     >  5:525970/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                       caca‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt                                                                                                         >  5:176103/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                       caca‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCATGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGtt                                                                                                                                     <  8:778794/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                               agccTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTa                                                                                                                                 <  7:657691/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
CGATTTCTACCGCAATGACGTGGCGAAAGCTAACCGTTCGCGGATGAACGTGCCGTTCCAGTTGGCGGACAGTGCGCTTGACAAATTGTTCCTTGAAGAGTCTTTTGCTGCTGGCCTTCATGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACA‑GCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACGGCAAAACAGTATTAACCAATCTGCTGGATAGCGATGACGTGCGCCATATGCTGAATGCATTAACAGCGTTAGGGGTAAGCTATACGCTTTCAGC  >  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: ‑