Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,826,379	0	T	C	32.1%	94.5 / 10.3	54	E182G (GAG→GGG)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (26/10);  new base C (1/16);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 4.80e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTC  >  NC_000913/1826132‑1826618                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                gAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGGTCTCACGATACCGACACCCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGGTCC‑A‑GCGGACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑GAAAC‑GGACCACCAGGGT‑CCTGATGGAACCCC‑AGC‑GCCTccc                                                                                                                                                                                                                                              >  1:342418/1‑250 (MQ=255)  aaGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAg                                                                                                                                                                                                                                             <  1:274491/251‑1 (MQ=255)    ggACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGcc                                                                                                                                                                                                                                           <  1:439380/251‑1 (MQ=255)     gACGGAGAGTGGCGAGTGATTTGGGAAACCACCCGATAACGGGGCGGCGGTATTCTCACAATACCGACACTCTCACCAGATAGCAGCGCCGCAATGGCGCTGGCAGTTACTGCAAACTTGCG‑AAGATCGTTTTGCCCAAA‑GGGGGACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGTCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑CTACCATCAGGTTTCGTGA‑AGAACATC‑AGC‑GGCTCCAGccc                                                                                                                                                                                                                                          <  2:548254/251‑1 (MQ=255)      aCGGCGAGTGGCGAGTAATTTGCGAAACCACCAGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGcccc                                                                                                                                                                                                                                         <  2:400419/251‑1 (MQ=255)                gcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGGTCC‑A‑GCGGACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑GAAAC‑GTACCACCAGGGT‑CCTGATGGAACACC‑AAC‑GCCCCCAGCCCCCCCGCAcccc                                                                                                                                                                                                                               >  2:457948/1‑250 (MQ=255)                 cgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTTCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGGGC‑TGAAAAGGG‑TAAAC‑GTACCACCAGGGT‑CCTGATGGAAAACC‑AGC‑GCCTCCAGCCCCGCCCCAccccg                                                                                                                                                                                                                              >  1:31462/1‑249 (MQ=255)                                        aaaggggggcggcgGTGTTCTCACGATCCCGACCCTCTCCCCAGAGAGCAGCGCAACATTGGCGCGGGCATTTCCTGAAAACGGGCG‑AGGATCGTTT‑GCCCAAAAGGGCACCGCTTTGC‑AAGTTCCCG‑GCGTACGGGCCAGCGCGCCAAAAGG‑TCTGCCGG‑AAAAAGGGG‑AAACC‑GTACCCCAAGGTC‑CCGGGTGAAACACT‑CGC‑GCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTccc                                                                                                                                                                                                       <  2:171184/243‑1 (MQ=255)                                          accGGAGCGGCGGTGTTCCCACGATCCCGACACTCTCACCAGATAGCAGCGCAGCATTTGCGCTGGTATTTACTGCAACCGGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAATTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TAAAATGGG‑TAAAC‑GTACCACCAGTTT‑CCGGATGAAACACC‑AGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAg                                                                                                                                                                                                     <  1:188113/249‑1 (MQ=255)                                                      gttttCTCACGATCCCGCCACTCTCACCAGATAGCAGCGCAGCAATTGCGCGGGGAGTTCCTGCAAACGGGCG‑AAGATCTTTTTGCCAAAA‑GGGCACCGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCCCAAAAATG‑TCTGGCGC‑TGAAAGGGG‑TAACC‑GTCCCCCCAGTTT‑CCTGAGGAAACACC‑AGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGGCGCGt                                                                                                                                                                                         <  1:201770/248‑1 (MQ=255)                                                        gttcccccGATACCAACACTCACCCAAAATAGCAGCGCAGAAATTGCGCGGGCAGTTACGGAAAACGGGGG‑AAGATCGTTTTGCCAAAG‑GGGCACCGCTTT‑CCAAGGTCC‑AAGCGTACAGCCCGGCCCGCCAAAATG‑TCGGGCCC‑TAAAAGGGT‑TAAAC‑GCCCCCCCAGGTC‑CCTGATGAAAAACC‑AGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGtt                                                                                                                                                                                        <  2:310625/243‑1 (MQ=255)                                                        gtTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTg                                                                                                                                                                                       >  1:535521/1‑251 (MQ=255)                                                                    ccGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑ACGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑G‑GCGTACAGGCCAGCGCGCCACAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACTTCAGAAATTTCTCGTCCCAGGGACTGTCTCGTTGCGTTAATACAcc                                                                                                                                                                           >  2:507874/1‑251 (MQ=255)                                                                     cGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCg                                                                                                                                                                          >  1:267333/1‑251 (MQ=255)                                                                     cGACACCCTCACGAGATCGAAGCGCAGCAATTGCGGTGGATGTTACTG‑AAACGGGGGGAAGTTTGTTTTGCCCAAA‑GGGCAACGTTTTGCCAATTTCC‑A‑GCGTACAGGCCGGCGCGCAAAAATG‑TCGGGCCC‑TGAAATGGG‑TAACC‑GTACCCCCAGGTT‑CCTGAGGAAACCCC‑GGC‑GCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCg                                                                                                                                                                          <  2:699991/251‑1 (MQ=255)                                                                      gACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCg                                                                                                                                                                          >  1:498747/1‑250 (MQ=255)                                                                              cACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCgg                                                                                                                                                                  >  1:45070/1‑250 (MQ=255)                                                                               cccAGATTGCAGCCCACCAATTCCTCTGGCAGTTCCTGCAAATTGGCG‑AAGACTTTTTTGCCCAAA‑GGCCAACTTTTTGCAAAGTTCC‑A‑GCGTACAGGCCTGCGCTCCAAAATG‑TCGGGCGG‑TAAAAGGGG‑TAAAC‑GGTCCGCCGTGTA‑CCTATTGGACCACC‑GGC‑TCCTCCGCCCCCGCCGCACCCATCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACTCCGAACTGCGGAt                                                                                                                                                                <  1:262640/250‑1 (MQ=255)                                                                               cccAGATAGCAGCACAGCAATTGCGCTGGCATTTCCTGAAAACTGGCG‑AAGATCGTTTGGCCAAAG‑GGGCAACGCTTTGCAAAGTTCC‑A‑GCGTACGGGCCAGGGCGCCAAAATG‑TCTGGCGC‑TAAAATGGG‑TAAAC‑GCCCCCCCAGTTT‑CCGGAGGGAACACC‑GGC‑GCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGAt                                                                                                                                                                <  1:688261/250‑1 (MQ=255)                                                                                                caATTGCGCTGGAAGTACCTGCAAACGGGCG‑AAGATCGTTTTGCCCAAG‑GGGCACCGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCAAAAATG‑TCGGGCGC‑TGAAATGGG‑TAAAC‑GTACCCCCGGGTT‑CCTGATGGAACCCC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAAtt                                                                                                                                               <  1:29209/251‑1 (MQ=255)                                                                                                aaatggcgcTGGAAGTTACTGAAAATTGGCG‑AAGGTCGTTTTGCCAAAA‑GGGCACCGCTTTGCCAAGTCCC‑A‑GCGTACAGGCCAGCGCGCAAAAATG‑TCTGGCGC‑TGAATGGGG‑AAAAC‑GTCCCACCGGGTC‑CCTGAGGGAACCCC‑AGC‑CCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAAtt                                                                                                                                               <  1:301466/246‑1 (MQ=255)                                                                                                          gggatttacggaaaacggggg‑AAGATCTTTGTCCCCAAA‑GGGCAACGCTTGGCCAATTTCC‑G‑GGGAACAGGCCAGGGCGCAAAAATG‑TCGGGCGC‑GAAAAGGGG‑AAAAC‑TTCCCACCAGTTC‑CCTGAGGAAACACC‑AGC‑GCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTa                                                                                                                                     <  2:64018/231‑1 (MQ=255)                                                                                                          ggCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGt                                                                                                                                      >  1:99565/1‑250 (MQ=255)                                                                                                           gCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGt                                                                                                                                      >  1:425231/1‑249 (MQ=255)                                                                                                                              g‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGcac                                                                                                                  >  1:4140/1‑250 (MQ=255)                                                                                                                                  gttCGTTTTGCCAAAG‑GGGCACCGTTTGGCAAAGTTCC‑A‑GCGACCAGGCCGGCGCGCAAAAATG‑TCTGG‑GCAGGAAAGGGG‑AAACC‑GAACCGCCGGGTT‑CCTGGTGAACCGCC‑AGC‑CCCTCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGAt                                                                                                              <  1:353934/249‑1 (MQ=255)                                                                                                                                     cgtgttGCCCAGA‑GGTCAACGCTGTGCAAAGTACC‑A‑GCGAACAGGCCAGCGCGCAAAAATG‑TCTGGCGC‑T‑AAATTGG‑TAAACCACCCCCCCAGTCT‑CTTGATGAAACCCA‑AGC‑GCCGCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCt                                                                                                           <  1:110395/247‑1 (MQ=255)                                                                                                                                       ttttGCCCAAA‑GGCCACCGCTTTGCCAAGTTCC‑A‑GCGAACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAACC‑GTCCCACCAGGTT‑CCTGATGAAACACC‑GGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCtt                                                                                                          <  1:224577/250‑1 (MQ=255)                                                                                                                                            cccAAG‑GGGCAACGCTTTGCCAATTTCC‑A‑GCGTACAGGCCGGCGCGCAAAAATG‑TCGGGCGC‑TGAAATGGG‑TAACC‑GTCCCCCCGGTTT‑CCGGATGGAACACC‑AGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCt                                                                                                    <  1:584460/251‑1 (MQ=255)                                                                                                                                            cccAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCt                                                                                                    >  1:539287/1‑251 (MQ=255)                                                                                                                                              ccaaagggcgctatttcgCCAATTTCC‑A‑GCGTACGGGCCGGCGCGCCAAAATG‑TCTGGCCC‑TAAAATGGG‑AAAAC‑GACCCCCCAGCTT‑CCTGATGAAAC‑CC‑AGCCGCCTCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTg                                                                                                   <  1:686812/234‑1 (MQ=255)                                                                                                                                              cAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑ACCGGCGC‑CGAAGTGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑ACC‑GCCTCAAGCACCGCCGCACCCAGCCCCGTCAGAAATTTCTCGCCCCAGGGAATGTCGCGGTGCGGTAATACACCTAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGCTGCCAGCGCCCAGATTCTTTGCCCTgg                                                                                                  >  2:363756/1‑251 (MQ=255)                                                                                                                                                aa‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTc                                                                                                >  1:706132/1‑251 (MQ=255)                                                                                                                                                aa‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTc                                                                                                >  1:714285/1‑251 (MQ=255)                                                                                                                                                               gCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTcc                                                                                   >  1:182826/1‑250 (MQ=255)                                                                                                                                                                        c‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGtt                                                                         >  2:588823/1‑251 (MQ=255)                                                                                                                                                                             gaaaagggcaaGCGCGC‑AAAATGGTCTGGCGC‑TAAAAGGGG‑TAACC‑GTCCCACCGGGTC‑CATGATGGACCCCC‑AGC‑GCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCa                                                                      <  2:31462/241‑1 (MQ=255)                                                                                                                                                                              gTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGa                                                                                                                                                                                                  >  1:420789/1‑126 (MQ=255)                                                                                                                                                                                aCAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGtt                                                                   >  2:617847/1‑251 (MQ=255)                                                                                                                                                                                aCAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGtt                                                                   >  1:565225/1‑251 (MQ=255)                                                                                                                                                                                 cAGGCCAGCGCGCCAAAAGG‑TTGGGCGC‑TGAAATGGG‑TAAAC‑GTCCCCCCAGGTC‑CCTGATGGAACCCC‑GGC‑GCTCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTc                                                                  <  1:651513/251‑1 (MQ=255)                                                                                                                                                                                  aGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTcg                                                                 >  1:589992/1‑251 (MQ=255)                                                                                                                                                                                    gCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGGGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTcgcg                                                               >  2:123600/1‑251 (MQ=255)                                                                                                                                                                                         cgcgcCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGccc                                                          >  1:309880/1‑251 (MQ=255)                                                                                                                                                                                            ccaaaaaTG‑TCTGGCCC‑TAAAAGGGG‑AAAAC‑TCCCCCCCAGGTT‑CCTGATGGAACCCC‑AGC‑CCCCCCGCCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGAc                                                       <  1:123600/248‑1 (MQ=255)                                                                                                                                                                                                  aTG‑TCGGGCGC‑AGAAATGGG‑TAACC‑GCCCCCCCGGGTT‑CCTGA‑GGAACCCCCAGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGttt                                                 <  1:333266/251‑1 (MQ=255)                                                                                                                                                                                                       gtGGCGC‑TAAAATGGG‑TAACC‑GTACCACCGGGTT‑CCTGATGGACCACC‑AGC‑GCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTcc                                             <  1:457948/250‑1 (MQ=255)                                                                                                                                                                                                         ggCGC‑T‑AAATGGGGTAACC‑TTACCACCAGGTT‑CCTGATGAAACCCC‑CGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTcc                                             <  1:588823/249‑1 (MQ=255)                                                                                                                                                                                                            gc‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCttt                                        >  2:163225/1‑251 (MQ=255)                                                                                                                                                                                                            gc‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCttt                                        >  2:711007/1‑251 (MQ=255)                                                                                                                                                                                                                     ggg‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAAt                                >  2:671517/1‑251 (MQ=255)                                                                                                                                                                                                                             c‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGcc                         >  1:614481/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccaggtccctgagaaaacccc‑AGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGAcc                  <  2:611101/232‑1 (MQ=255)                                                                                                                                                                                                                                                      ccccc‑AGC‑CCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTc  <  1:445473/248‑1 (MQ=255)                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                GAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCG‑AAGATCGTTTTGCCCAAA‑GGGCAACGCTTTGCCAAGTTCC‑A‑GCGTACAGGCCAGCGCGCCAAAATG‑TCTGGCGC‑TGAAATGGG‑TAAAC‑GTACCACCAGGTT‑CCTGATGGAACACC‑AGC‑GCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTC  >  NC_000913/1826132‑1826618

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