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

GAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGTGAACAT  >  AE015451/1124306‑1124875
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                             
gAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGAACGGTACCGAGAGTTCGTC‑GCCGTCTTCGCCGGTCTAGCCCGACGGGCTGAccaaaccgtccaaccccagc                                                                                                                                                                                                                                                                                >  2:742956/1‑282 (MQ=255)
             ggacgcGGGGAGCCGGCACAACGAGCACCGACGGCACGGCGCGCCTTGGAGATTGACGACACGAGGATGCCATCGACCGGCGTAGTCTCTATGACCAGCTCGTGGCCATGAAGGAACAGGCCAGCTTCCAGGCGCAGCGATGT‑GCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCATTCCACACCCAGCAGTTTTACCcgga                                                                                                                                                                                                                                                                   <  2:297050/296‑2 (MQ=255)
                     gAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGGCGGCCTTCTGC‑TTGTAGAACAAGAGGCTGACGAAGCAAGCCACACCCAGCCGTTTTACCCGGCGGAACTGc                                                                                                                                                                                                                                                           >  4:566352/1‑301 (MQ=255)
                                     gaaccgtagacacggctcaacgtGGAAAGTGCCGACAGCAGGCTGTCTGCGATCGGCGTAGTCTCGCTGTTCATGGCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGGGGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTGTTACCCGGCGGAACTGCATGAAGGTCATGCCgg                                                                                                                                                                                                                                           <  4:743216/279‑1 (MQ=255)
                                           tCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGCTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCCGATGACGAAGCAGTCGCCACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATggggagagggg                                                                                                                                                                                                                                     >  2:77140/1‑291 (MQ=255)
                                                                         ccgctcttccgcacGATCGGCGTAGTCTCGCTGTTAATGTCGGGTCCATACAGGCACAGGCCAGCTTCCAGGCGCATCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTg                                                                                                                                                                                                       <  1:614527/288‑1 (MQ=255)
                                                                                        aTCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCTTCTTCTGC‑TTGTAAAAAGAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTgagatcggaagagcac                                                                                                                                                                                        >  3:614529/1‑286 (MQ=255)
                                                                                                 gTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGACAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGAATGAAGGTCATGCCGGCCACTGCTGGAGCCAGGCGTGACAGCACCTTGACTGCAGCCGGGcgtatcagggcggtcc                                                                                                                                                                               >  4:526895/1‑286 (MQ=255)
                                                                                                                      ccATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑GTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTCTCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGTGCAGGGCACGCTCCTCTAACatcg                                                                                                                                                          >  4:688360/1‑298 (MQ=255)
                                                                                                                                        gCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTC‑CGGTTGAAGAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGc                                                                                                                                        >  4:565009/1‑301 (MQ=255)
                                                                                                                                                      cGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCGTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCCGCTCCACCTCgatcggaagagctccctgtacg                                                                                                                          >  2:26213/1‑280 (MQ=255)
                                                                                                                                                        ctcagcgacgcttttcagatatacgggtGTACTGCCGGCTGGGACATCAGGCGGCGGGCCAAGGGTTCGGCATCGGTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGttt                                                                                                                         <  2:348263/274‑1 (MQ=255)
                                                                                                                                                                           ccggtcctggTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTg                                                                                                     <  4:726847/292‑1 (MQ=255)
                                                                                                                                                                            gATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGCCTTCTGC‑TTGAAAAAAAAAAGGCAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGTAACTGCATGAAGGTCATGCCGGCCACTTCTGTAGCCAGGCGTTCCAGCACCTTTACTGCAGCCGGGCCTTGCATGGCGGTCAGGGCACCTTCCTCGACCAGCGGCTTCACCTCCCAACGCCGTACGATGTGGGTTTGCA‑TGTGAGCCAGGTCCTGc                                                                                                    >  4:623167/1‑301 (MQ=255)
                                                                                                                                                                              tCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGtttagatcggaagagcacacgtctg                                                                                                   >  4:348347/1‑279 (MQ=255)
                                                                                                                                                                                  cTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGAA‑GGTGAGCCAGGTCCTGCTCCttt                                                                                              >  3:66567/1‑300 (MQ=255)
                                                                                                                                                                                   tGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTGGCACGGT‑AGCCAGGTCCTgagagcggg                                                                                             >  2:726587/1‑293 (MQ=255)
                                                                                                                                                                                              cTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAACAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGAGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGGTCCAGCACCTTAACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGAGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCATGCACCTTGCAGGCGGCGTTg                                                                                  >  4:779416/1‑301 (MQ=255)
                                                                                                                                                                                                     aTCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTTACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGttaaatcggg                                                                           >  3:621671/1‑293 (MQ=255)
                                                                                                                                                                                                                          cagacgtgtgctcgtccgatctagaTCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGACCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGATCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCa                                                      <  4:264418/279‑1 (MQ=255)
                                                                                                                                                                                                                               ggttcagacgtgtgctcttccgatctA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTg                                                 <  3:88560/276‑1 (MQ=255)
                                                                                                                                                                                                                                         ttcagacgtgtgctctgccgatctttCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATc                                        <  4:368020/277‑1 (MQ=255)
                                                                                                                                                                                                                                           acgtgtgctcttccgatctgTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAg                                      <  3:300714/282‑1 (MQ=255)
                                                                                                                                                                                                                                              cccctacacgacgctcttccgatctggTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtc                                   <  2:128837/276‑1 (MQ=255)
                                                                                                                                                                                                                                                agaTCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTACTTGCAGGCGTCGTTGACGACCAGGAACAGGGTGTCCTCGccaaagatcggaagagcggcgtggag                                >  2:264352/1‑277 (MQ=255)
                                                                                                                                                                                                                                                         a‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTgagatcggaagagcgtcgtgtggggaa                       >  1:88557/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              gTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAgagatcggaagagcgtcgtg                   >  1:300715/1‑282 (MQ=255)
                                                                                                                                                                                                                                                                       ttttgaGCCCGAGCGGCTGAAGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTGCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTg          <  2:666499/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ggTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtcagatcggaagagcacacgtctgaac          >  4:128866/1‑276 (MQ=255)
                                                                                                                                                                                                                                                                         tgtAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTgag        >  4:666721/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                               ccGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGTGAACAt  <  4:156521/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                             
GAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGTGAACAT  >  AE015451/1124306‑1124875

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