Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,124,5710TA41.7% 29.2 / 21.8 37Y191F (TAC→TTC) gcvT‑Iaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (15/6);  new base A (15/0);  total (30/7)
Fisher's exact test for biased strand distribution p-value = 3.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCAT  >  AE015451/1124311‑1124865
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      
gggCCTCGGCAGCAGGGAAGCCGGCAGCACGTGGCAAGTCGACACGGCGAACCTTGGAAAGTGCACACAGAAGGCTGGCGTCGATCGGCGTGGGCTCGCTGTTCATGTCGGGGCCATACAGGCACAGGCCAGCTTCCAGGGGCAGCGAGTCGGGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGGAGCCCGGGCGGCTGACGAAGCAGTCCACACCCAGCAGttt                                                                                                                                                                                                                                                                <  2:1047548/299‑1 (MQ=255)
       ctgccctacgcgattccgatatcGTGCCCCGTCGCAACGCCGCACCCTGAAAATTGCACACGGCAGGCGGGCTTCGATCGGGGTAGTCTCGCTGTTCATGTGGTGGCCATACAGGAACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGACATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCcgg                                                                                                                                                                                                                                                         <  2:1682225/279‑1 (MQ=255)
                 aaGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAGAGGCTGACGAAGCAGTCCACACCAAGCAGTTTTACCCGGCGGAACTGcc                                                                                                                                                                                                                                               >  1:563666/1‑300 (MQ=255)
                      ggCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTTTTCCGCTTGTAAACCAAACGGCTGACGAACCAGTCCCCACCCAGCACTTTTACCcggcagatcggaagagcc                                                                                                                                                                                                                                          >  1:74774/1‑287 (MQ=255)
                           caCGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCCGCGTGAAAAAAGAACGGCTGACGAAGCAGTCCAACCCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCAt                                                                                                                                                                                                                                     >  1:1710226/1‑301 (MQ=255)
                            aCGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCGGTCTTCTGGGTGTAACACAAGCGGCTGACCACGCAGTGCACACCCAGCAGTTGTACCCGGCGGagatgcagggaggtcctt                                                                                                                                                                                                                                    >  1:2232397/1‑284 (MQ=255)
                             cGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCCGTCTTCCCGGTGTAACAAGAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCcggagatcggaagagcacacgtctg                                                                                                                                                                                                                                   >  1:1682879/1‑279 (MQ=255)
                                                  cACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCGGTTTTCCCCTTGTATACCAAACGGCTGACGAAGCAGTCCACACCGAGCAGTTTTCCCCGGCGGAACTGAATGAAGGTCATGCCGGCCAcagatcggacgcgg                                                                                                                                                                                                               >  1:220432/1‑288 (MQ=255)
                                                      ttGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAACAAACGGCTGACGAAGCAGTCATCACCCAGCAGGTTTACCTGGCGTAACTGCATGACGGTCATGCCGTCCACTTCTGTAGCCAGGCGtt                                                                                                                                                                                                           >  1:479962/1‑301 (MQ=255)
                                                        ggAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGGCTGACAAAGCACTCCACACCCAGCAGATTTACCCGGCGGAACTGCATGAAGGTCATGACGGGCACTTCTGGAGCCAGGCCTTcc                                                                                                                                                                                                         >  1:1434754/1‑301 (MQ=255)
                                                          aaaTTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAGAAAAGCAGACGAAGCAGTCCACACCCAGCAGCTTTACCCGGCGGAACTGCATGAAGGGCATGCCGGCCAATTCTGGAGGCAGGCGTTCCaa                                                                                                                                                                                                       >  1:986310/1‑300 (MQ=255)
                                                                   acaGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTGCGCTTGAAACACAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGACGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGAc                                                                                                                                                                                              >  1:713133/1‑301 (MQ=255)
                                                                                ttCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACTAAGCAGTCCACACCCAGCAGTTTTACCCGGCGAAACTGCCTGAAGTTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCACCACCTTGAATCCATCCGGGccc                                                                                                                                                                                 >  2:55672/1‑300 (MQ=255)
                                                                                        gCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGACACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACGGCAGCCTGGCCttccacggg                                                                                                                                                                         >  1:2020489/1‑294 (MQ=255)
                                                                                                  gCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTCCTTGTAAAACAAAAGGGCGACGAAGCAGTCCACAACCACTCGTTCTACCCGGCGGCACCGCATGAACGTCATGCCTGCCACCTCTTGAGCCAGGCGTTCCAGCCCCTTTACCCCAGCCGCGCCTTGCCCGGCCAGCAggccc                                                                                                                                                               >  1:1929962/1‑298 (MQ=255)
                                                                                                              gTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCCCTTCTGGAGCCAGGCGTTGCAGCACCTTCACTGCAGCCGGGCCTTGCCGGGCCAGACGTGCACGCTCCTCCAAc                                                                                                                                                   >  1:108204/1‑301 (MQ=255)
                                                                                                                ggCCATACAGGCACAGGAGCGCTTCCATGCTCATCGAGGCGCGTGCTCACAAGACTATCGGCTGTACTTCCTGCTCGGGCATCAGGCTGAGGGCCAAGGCGTCGGCTGCGTTGACCGGTACGGAGCTCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAg                                                                                                                                                 <  2:108191/301‑1 (MQ=255)
                                                                                                                                                         cgTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCCGCTTGAAAAAAAAAAGGCTTACGAAGCATTCCACTCCCATCCTGTTTACCCGTCGGAACTGCACGCAGGTTATGCCGGCTACTTCGCAAGCCAGGCGTTCCAGCACCTTGCCTGCAGCCGTGCCTCGCAGGCCGACAACGTCACGCTCCTCGACCAGCGGCTGCACCTCGTAAccctcccccagttcggttttcactt                                                                                                        >  1:1944139/1‑277 (MQ=255)
                                                                                                                                                              gcCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAGCTGACGAAGCAGTCCCCACCCAGCCGTTTTACCCGGCGGAACTTCATGAAGGTCATGCCTGCCACTTCTTGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCTAGCAGGGCACGCTCCTCGAACAGCTGCTGCACCTCGCACCGGCTGCCGCTGTggtgttgccggttagcc                                                                                                   >  1:1709133/1‑286 (MQ=255)
                                                                                                                                                                                 cctcttccgctctcccacatgcggcgGGACCAGGCGTCGGGAGCGTTGACCGGTACCGAGATTTCGTAGCCGTCTGCTGCAGTGTAGCCCGAGCGTCTGACGAAGCAGTCCAACCCCATCAGTTTTGCCCGGAGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCATTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAgg                                                                                <  2:1522301/281‑1 (MQ=255)
                                                                                                                                                                                 aCTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGGCTGACGAAGCCGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGTCCACTTCTGTAGCCAGGCGTTCCAGCACCTTGACTGCATCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCCACGGCGGCCGATGTGGTTTTGCAGGTGAGCCCGGTCCCGCTCCTTGCAgg                                                                                >  1:1242418/1‑301 (MQ=255)
                                                                                                                                                                                    tCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTTGAggcgg                                                                             >  2:2264450/1‑301 (MQ=255)
                                                                                                                                                                                              ccATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCGTGAAAAAAGAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGGTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAggagatcggaagaga                                                                   >  1:1522890/1‑288 (MQ=255)
                                                                                                                                                                                               cATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGTAAAACAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTTCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTCCATGCGGCGTTgacgccc                                                                  >  1:2158272/1‑298 (MQ=255)
                                                                                                                                                                                                  cAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAACGACTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGTAGCCAGGCGTTCCAGCACCTTGACTGCATCCGGGCCTTGCAGGTCGAGCAGTGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGTGTTTGCAGGTGAGCCCGGTCCTGCTCCTTTCAGTCTGCGTTgacgtcccga                                                               >  1:1554314/1‑295 (MQ=255)
                                                                                                                                                                                                            gCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATAAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTTCAGCCGGGCCTTGCAGTGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTTAGCCAGGTCCTGCTCCTTGCAAGCGGCGTTGCCGACCAGGAACAgggggt                                                     >  1:2179135/1‑298 (MQ=255)
                                                                                                                                                                                                            gCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCGGCCGATGTGGGTTTGCAGGTTGGCAAGGTCCTTCTTCTTGCAGGCGGCGTTGACGACCAGGACAAGGgtgt                                                     >  2:660267/1‑301 (MQ=255)
                                                                                                                                                                                                                       ggCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCCGGCGGCGTTGACGACCAGGCACAGGGTGTCCTCGCCCAgg                                          >  1:2058092/1‑301 (MQ=255)
                                                                                                                                                                                                                        gCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGGTGGTTTGCATGTGAGCCAGGTCCTGCCCCTTGCAGGCTGAGTTGACGAACAggcaacgagtgtcccagcctaggg                                         >  2:2053/1‑279 (MQ=255)
                                                                                                                                                                                                                                        ccGAGATTTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTGTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTGCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtc                         <  1:660502/301‑1 (MQ=255)
                                                                                                                                                                                                                                           ataTCTCATAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAAATGCATGAAGGTCATGCCGGACACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCTAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGTGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtcgtc                      <  1:2265369/299‑1 (MQ=255)
                                                                                                                                                                                                                                             aTCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGGTGGCAACCATACGGtcgtcgg                    >  2:543769/1‑300 (MQ=255)
                                                                                                                                                                                                                                                  gTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCGGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTTGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTTCAGGCGGCGTTGACGACCAGGAAACGGGTGTCCTCGCCCAGGTTGGCAACCATCAgggcttcgtacatg               >  2:303477/1‑289 (MQ=255)
                                                                                                                                                                                                                                                   tAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACCGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATgg              >  2:1644418/1‑300 (MQ=255)
                                                                                                                                                                                                                                                         cGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTGGTCGGGaagggcgcc          >  1:1981657/1‑294 (MQ=255)
                                                                                                                                                                                                                                                           tCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCtt        <  1:543964/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ccGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCGCCTCGAACAGCGGCGGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCCCCTTGCAGGCCGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCAGCAGGTCGGCGAGAATGCCGCCTTGCTCAt  >  2:1784359/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      
TCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCAT  >  AE015451/1124311‑1124865

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