Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,124,5750CA43.9% 24.5 / 32.5 41G190C (GGC→TGC) ‡gcvT‑Iaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (12/11);  new base A (18/0);  total (30/11)
Fisher's exact test for biased strand distribution p-value = 7.57e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGA  >  AE015451/1124286‑1124849
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                   
acacCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCG‑CCTCCCCGGTGTAGCCCAAGAGGCTgc                                                                                                                                                                                                                                                                          >  4:245901/1‑300 (MQ=255)
     ctctcccgaactggCGAATAGCGCCTCGGCACCCGGGAAGCCGGAAGGACGGGCACCGTCGGCACGGCGCACCTTGGA‑AAGTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGGGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAg                                                                                                                                                                                                                                                                     <  1:186035/289‑1 (MQ=255)
           cACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGGGTTGACCGGTACCGAGAGCTCGTA‑ACCGTCATCGCGGGTGTAGCCCGAGCgctcgaggaagaactcc                                                                                                                                                                                                                                                               >  2:504826/1‑285 (MQ=255)
              gTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTTACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTccccc                                                                                                                                                                                                                                                            >  1:178546/1‑298 (MQ=255)
                      gacgcgcagcagcccccaggaggCCGGCAGAACGTGAGCCGGCGCCACGGCG‑ACCGTGGAGAATTGCCCGGCGCAGGCTGGGTTCGATCGGCGAAGGCTCGCTGTTCATGTCGTGTCCAGACAGGCACAGGCCAGGTTGCCGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGGGGCAGTCCACACCcagcag                                                                                                                                                                                                                                                    <  2:245857/279‑1 (MQ=255)
                        ttCGCCTCGGCACCCGGGAATCCGGCAGAACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTTGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGtt                                                                                                                                                                                                                                                  <  3:178552/300‑1 (MQ=255)
                                               ggCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAg                                                                                                                                                                                                                           >  4:605166/1‑301 (MQ=255)
                                                              gtgggcagtagcaaCCTTGGA‑AATTGCCACCAGCGGGCTGGTTTCGATCGCCGTGTTCTAGCTGTTCATGCCGCGGCCATACCGGCACAGGCCGGCTTCCGGGGCCAGCGAGTCGCGTTCGCCAAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCGGGCGGCGGGCCAAGGCGTCGGCAGCGTTGCCCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACt                                                                                                                                                                                                            <  4:504952/288‑1 (MQ=255)
                                                                cGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTT‑GTGCTTGAAAAAAAAACGGCTGACGAAGCAGGGCAAACCCAGCAGTTTTACCCGGCGGCACTGACTGAAGGTCaagacagatacgag                                                                                                                                                                                                          >  3:335573/1‑288 (MQ=255)
                                                                 ggCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑GTGAAAACAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGGTTTACCAGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCt                                                                                                                                                                                                         >  3:292486/1‑301 (MQ=255)
                                                                                 gA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGATGACGAAGCAGTCCACACCCCGCAGTTTTACCCGTCGGAACTGCATGAAGGTAATGCCGGCCACTTCTGGAGCCAGGCGTTccc                                                                                                                                                                                         >  4:717857/1‑300 (MQ=255)
                                                                                          ccACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCCGTTTTACCCGGCGGAACTGCATGGAGGTCATGCCGGCACCTTCTGGACGCAGGCGTTCCAGCACC‑ttt                                                                                                                                                                                >  2:669467/1‑300 (MQ=255)
                                                                                             cagcagGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAAAGACTGAAGAAGCAGTCCACACCCAGCAGTTTTACCCGTCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCAAGGCGTTCCAGCACC‑TTGACt                                                                                                                                                                             >  3:278917/1‑301 (MQ=255)
                                                                                                      ggctTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGGCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTTACGAAGCAGTCCACACCCCGCAGTTTTACCCGGCGGCACTGCATGAAGGTCATGCCATGCACTTCTGGAGCCCGGCGTTCC‑GCACCATTGCCTGCAGcccgg                                                                                                                                                                    >  2:214148/1‑298 (MQ=255)
                                                                                                        ctTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAACAAAAAAGGCTGACGAAGCATTCCCCACCCAGCACCTATACACGGCATAACTGCCTGACTGTCATGCCGCCCAGTTCTGTATCCAGGCGTTCCACCACC‑TTGACTGCATCCGGGcc                                                                                                                                                                  >  4:179634/1‑301 (MQ=255)
                                                                                                         ttCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCGGACGAAGCAGTCCACACGCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACGTCTGGATGCAGGCGTTCCAGCACC‑TTGACTGAAGCCGGGCCt                                                                                                                                                                 >  3:447992/1‑301 (MQ=255)
                                                                                                             aTCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGCCTTCTGC‑TTGTAGAACAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGCCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGCCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGcg                                                                                                                                                             >  4:7207/1‑300 (MQ=255)
                                                                                                                          cacgcgatactgccgttctcaacggcaccaggcCGGCTTCCACGCGCAGCCAGTCGCGTGAGCGAAAGCCGATCTGCTGGACATCCGGCTGGGCGATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATGTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACCCCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGc                                                                                                                                                <  1:413801/273‑1 (MQ=255)
                                                                                                                              gTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAAGCAGAAGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGGATTAAGGGCATGCCAGGCACTGCTGGAGCCAGGCGTTCCAGCACC‑TTGAATGCAGCCAGGCCTTTCATGGCGaaccgtgcatgc                                                                                                                                            >  3:103249/1‑290 (MQ=255)
                                                                                                                                             aTACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTC‑CGCTTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATTCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCagatcggaagagcacacgt                                                                                                                             >  3:413811/1‑283 (MQ=255)
                                                                                                                                                   gcacaggcCGGCTTCTGGGCGCAGCGGGCTGCTTGCGCTCAAGCCGATGCGCTGTACGTCCGGCCCGTGCAGCACGCGTCTGGCCAAGGCGTAGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGCCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGACGGTCATGCCTGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCAcc                                                                                                                       <  1:292480/301‑1 (MQ=255)
                                                                                                                                                                         aGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCACCGGCTGCCGCTGTggg                                                                                                 >  4:615338/1‑301 (MQ=255)
                                                                                                                                                                            ggatcgcgcgcgccctacacGATTGGCTGTCCTTCCGTCTCGGCAATCAGTGGTCGGGCCAGTGCGTGGCCAGCGTTGACCGGTACCGGGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGCGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGTAACTGCATGAAGGTCATGGCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCGTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGttt                                                                                              <  2:7207/282‑1 (MQ=255)
                                                                                                                                                                                gctggtggccgaggCCGAACGGATGTACATCCGGATAGGCAAGCAGGCGTCGGGCGAAGGCGTGGGCAGAGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGAAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGAACCTCGCAACGGCTGCCGATGTGGGTTTGCAg                                                                                          <  3:125874/288‑1 (MQ=255)
                                                                                                                                                                                   gTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTCC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGCACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTg                                                                                       >  4:665843/1‑301 (MQ=255)
                                                                                                                                                                                   gTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCATGGCGAGCATGGCACGCTCCTCGAACAGCGGCCGCACCTCGCAACGGCTGccactcggtgtttcgagagc                                                                                       >  1:125870/1‑283 (MQ=255)
                                                                                                                                                                                       gcCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGTCCACTTCTGGAACCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGcc                                                                                   >  3:246798/1‑301 (MQ=255)
                                                                                                                                                                                                        gTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGTCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTTCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCTAACAGCGGCTGCACCTCGCCACGGCTGCCGATGTGGGTTTGCCGGTGAGCCAGGTCCTGCCCCtctcc                                                                  >  4:433433/1‑297 (MQ=255)
                                                                                                                                                                                                            ttCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAACAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAcgcg                                                              >  4:90017/1‑301 (MQ=255)
                                                                                                                                                                                                                     ggCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAACAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCTGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTgacg                                                     >  4:60415/1‑301 (MQ=255)
                                                                                                                                                                                                                     ggCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTgacg                                                     >  1:160603/1‑301 (MQ=255)
                                                                                                                                                                                                                           cAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGTCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGCTTTGCAGTTTGGCCAGGTCCTGATCCTTGAAGGCTGCGATGACGAcccgg                                               >  3:268100/1‑298 (MQ=255)
                                                                                                                                                                                                                            aGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTTGGTTTGCAGGTTATCCCGGTCCTGCTCCTTGCAGGCCGCGTTTACGACCAGGa                                              >  4:256901/1‑301 (MQ=255)
                                                                                                                                                                                                                                ggcggGCCAGGGCGGCGGCAGCGTTGACCGGTACCGAGATCGCGTA‑GCCGTCTTCGCCGGGGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAg                                          <  4:214186/301‑1 (MQ=255)
                                                                                                                                                                                                                                 gcggGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAgg                                         <  3:160609/301‑1 (MQ=255)
                                                                                                                                                                                                                                         aGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAGAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGAGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTc                                 >  4:128998/1‑301 (MQ=255)
                                                                                                                                                                                                                                             gTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTTAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGccc                             >  4:138957/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         gACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGTGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTCTGCGCGACCAGGTAGGCAAcc                 >  2:643319/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                tccgatctctcGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgt          <  1:443992/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ctcGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCCGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATccgtgcgtcagaa  >  2:302711/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                        ctcGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtagatcgga  >  3:444004/1‑293 (MQ=255)
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                   
ACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGA  >  AE015451/1124286‑1124849

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