Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,124,577	0	G	A	28.1%	45.4 / 12.6	32	S189L (TCG→TTG)	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base A (9/0);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 1.17e-02
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‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGTGAACATG  >  AE015451/1124306‑1124876                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                             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)                                                                                                                                                                                                                                                                                 cGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCCCATTGGTTAActta  >  2:156483/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                             GAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGTGAACATG  >  AE015451/1124306‑1124876

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: ‑