Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,008,373	0	G	T	37.5%	55.3 / 11.5	40	G33C (GGC→TGC)	yedE	putative selenium transporter YedE
Reads supporting (aligned to +/- strand):  ref base G (9/16);  new base T (15/0);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 3.64e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCATTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGC  >  NC_000913/2008124‑2008620                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          ggCATTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTGGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTTCATTCATTTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTgg                                                                                                                                                                                                                                                          <  1:41177/251‑1 (MQ=255)    aTTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAATTTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TTTTCTCTC‑TTCTTTAT‑ATTTTTGCa                                                                                                                                                                                                                                                        >  2:52826/1‑250 (MQ=255)     ttAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑ATTTTTGCAtt                                                                                                                                                                                                                                                      >  1:265635/1‑251 (MQ=255)                                          gagaAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTTTCTC‑TACTTACT‑TTTTTTGCATTTCTAGCACCTTTTTGGCTGTCTC‑GGGGGAATTTTccc                                                                                                                                                                                                                >  2:85865/1‑251 (MQ=255)                                                      cTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑CGTTTTCTCCTA‑TTACTCATTTTGCCAT‑ACTTGCACCCTTTGGGGTGTCACCGGGGCAATTT‑CCCGTTTGggcggtc                                                                                                                                                                                                    >  2:92052/1‑249 (MQ=255)                                                              tattTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCt                                                                                                                                                                                             <  2:265635/250‑1 (MQ=255)                                                                    tgtgAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACt                                                                                                                                                                                       <  1:383820/250‑1 (MQ=255)                                                                                    gAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑ATTTTGTCATTACTGGCACCTTTTGTGCTTTCAC‑GGGTGAATTTACCCGTTGTTGCGGCCACCCCCCGCAACTTTTTGGCGTCCATGCtt                                                                                                                                                                      >  2:455353/1‑250 (MQ=255)                                                                                         taaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑TTTTTGGCCTTACTGGCACCCTTTGGGGTGTCAC‑CGGTGAATTTATCCCTTGGGGCGGCCAGCCCCTGCTACTGGTCGGCggcccagcttaagcg                                                                                                                                                                 >  2:443087/1‑237 (MQ=255)                                                                                           aatTTGCATGGGGGATGATTCGCGTCTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACCCGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGATCATCGCGGCTGG‑TATTCTCTC‑TACTTGCT‑GGTTGGGCATGACTGGCACCTTGTGGGCTGGCGC‑GGATGAATTTGCCCGTTGGGGCGTCCAGATCCTGCAACTGTTCGGCGTCCAGCCTGAAGAGTg                                                                                                                                                               <  2:477350/251‑1 (MQ=255)                                                                                            atTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGAATTCT‑TC‑TTCTTACT‑ATTTTTGCATTTCTGGCCCCTTTTGGGGTGTCCC‑CGGGGAATTTACCCCTTTGGGCGGCCAGCTCCTGCAACTGTTTGGCGTCCATGCCGAagagggg                                                                                                                                                              >  2:159367/1‑248 (MQ=255)                                                                                            atTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTTTCTC‑TACTTACT‑ATTTTTTCCTTATTGGCCCCTTTTTGGCTGTCAC‑CGGGGGATTTACCCGTTTGGGCGGCCCGCTCCTGCAACCGTTCGGCGTCCATGCTGAagagggg                                                                                                                                                              >  2:289622/1‑248 (MQ=255)                                                                                                  aTGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGc                                                                                                                                                                        <  2:143165/235‑1 (MQ=255)                                                                                                  aTGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑ATTTTTGCATTACTGGCACCCTTTGGGCTGTCAC‑GGGGGAATTTACCCGTTGGGGCGGCCAGCCCCTGCAACTGTTCGGCGTCCATGc                                                                                                                                                                        >  1:143165/1‑235 (MQ=255)                                                                                                     gtgtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTTCT‑ATTTTTGCATTACTGGCACCTTTTGGGCTGTCCC‑GGGTGAATTTACCCCTTGGGGGGGCCCGCCCCTGCAACTGTTTGGGGTCCCTGCTGAAGGGGGGGG‑GTACttt                                                                                                                                                    >  2:551908/1‑251 (MQ=255)                                                                                                                  gACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTTTCCC‑TTATTACT‑TTTTTTGCATTTCTGGCACCCTTTGGGCTGGCAC‑GGGCGTATTTACCCGTTGGGGCGGGCACCTCCTGCAACCGTTCTGCGTCCATTCTTGAGAGGGGGG‑GTACTTTTAAAATATccctt                                                                                                                                       >  2:417940/1‑248 (MQ=255)                                                                                                                                   ttCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTTTCTC‑TAATTACT‑TTTTTTGCATTACTGGGACCTTTTTGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCCGCTCCTGCAAATGTTCGGCGTCCATGCTGAAGAGTGGGG‑GTACTTTTAAATTATCCATCTGGAAGGGTCGCCAtt                                                                                                                       >  1:437506/1‑250 (MQ=255)                                                                                                                                                ggCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCCACTGTTCGGCGTCCCTGCTGAAGAGTGGGG‑GTTCTTTTAAATTATCCATCTTGAAGGATCGCCATTAACCCGCCTCGACg                                                                                                         >  2:272555/1‑251 (MQ=255)                                                                                                                                                               atgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TTCTTCCT‑CTTTTTGCATTACTTGCACCTTTTTGGCTGTCAC‑GGGTGAATTTACCCCTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGGCCATGCTGAAGAGTGGGG‑TTTCTTTAAAATTATCCCTCTGGCAGGATCGCCATTAACCCCCCTCGACGGGATGATGATCCTcc                                                                                          >  1:185240/1‑250 (MQ=255)                                                                                                                                                               atgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACCTACT‑ATTTTGGCCTTACTGGCCCCTTTTTGGCTGTCAC‑GGGGGAATTTACCCGGTGGGGCGGCCCGCTCCTGCCACCGGTCGGCGGCCCTGCTGAAGAGGGGGG‑GTACTTTAAAATTATCCCTCTGGGAGGATCGCCATTAACCCGCCTCCACGGGAGGAAGATACTCg                                                                                          >  2:59626/1‑251 (MQ=255)                                                                                                                                                                 ggCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATc                                                                                                                                           <  1:479514/200‑1 (MQ=255)                                                                                                                                                                 ggCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TTCTTACT‑ATTTTTGGATTACTGGCACCTTTTTGGCTGTTAC‑GGGGGAATTTACCCGTTGGGGCGGGCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTTAAATTATc                                                                                                                                           >  2:479514/1‑200 (MQ=255)                                                                                                                                                                   cagcaATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATCTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTAt                                                                                      <  2:419444/251‑1 (MQ=255)                                                                                                                                                                            aaaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGGGGATTTTCCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGGCCCTGCTGAAGGGGGGGG‑GTACTTTAAAAATTTCCATCTGGAAGGATCGCCCTTTACCCGCATCGACGGGATGATGATCCTCGGTATGttttgggg                                                                             >  2:270751/1‑246 (MQ=255)                                                                                                                                                                              acacGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑ATTTTTGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGGGAATTTACCCGTTGGGGCGGCCCGCTCCTGCAACTGTTTGGCGTCCATGCTGAAGAGGGGGGGTTA‑TTTAAAATTTTTCCTCTGGAAGGATCGCCATTAACCCCCATCGACGGGATGATGATCCTCGGTATGTTTGGCGgct                                                                           >  2:548233/1‑251 (MQ=255)                                                                                                                                                                                   ccTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCttt                                                                      <  1:79452/251‑1 (MQ=255)                                                                                                                                                                                      ggTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGcc                                                                   <  2:185240/251‑1 (MQ=255)                                                                                                                                                                                      ggTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑ATTTTTGCATTACTGGCACCTTTTGGGCTGTCCC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGGGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAACCGCATCGACGGGGTGATGATCCTCGGTATGTTTGGCGGCTGCTTTGcc                                                                   >  1:311201/1‑251 (MQ=255)                                                                                                                                                                                       gTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTC‑TACTTACT‑ATTTTGGCCTTACTGGCACCTTTTGGGCTGTCAC‑GGGGGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGGAGAGGGGGG‑GTACTTTAAAATTATCCATCTGGAAGGATCGCCCTTAACCCGCATCGACGGGGTGATGATTCTCGGTATGTTTGGCGGCTGCTTTggcc                                                                  >  2:107272/1‑248 (MQ=255)                                                                                                                                                                                              aaaTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTTTCTC‑TTCTTACT‑ATTTTTGCATTTCTGGCACCCTTTGGGCTGTCCC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCCCCTGCCACTGTTCGGCGGCCATGCTGAAGAGGGGGG‑GTACTTTAAAAATATCCCTCTGGAAGGATCGCCCTTAACCCGCATCGACGGGATGATGATCCCCCGTATGTTTTGCGGGTGCGTTTCCGCAgcgcg                                                           >  2:545725/1‑250 (MQ=255)                                                                                                                                                                                                aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑GATTCTCTT‑TACTTACT‑CTTTTGGCATTACCGGCACCTTTTGGGCTGGCCC‑GGGGGAATTTAACCCTTGGGGCGGCCCGCCCCTGCCACTGTTCGGCGGCAATTCTGAAGAGGGGGG‑GTACTTTTAAATTATCCATCTGGGAGGATCGCCATTAACCCGCATCGGCGGGGAGATGATCCTCGGGATGTTTGGGGGCTGCTTTTCCGCCgcgccgt                                                         >  2:172336/1‑248 (MQ=255)                                                                                                                                                                                                            cATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTggc                                                                               <  1:544268/217‑1 (MQ=255)                                                                                                                                                                                                            cATCCCTGCGGGCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑CTTTTTCCCTTACTGGCACCTTTTGGGCTGTTAC‑GGGGGGATTTTCCCGTTGGGGCGGCCAGCGCCTGCAACTTTTCGGCGTCCCTGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAAGATCGCCATTAACCCGCATCGACGGGATGATGTACCTCGGTATGTTTTGCGGCTGCTTTTCCGCAGCCCTGGGGGCCAACAATg                                             >  2:154856/1‑251 (MQ=255)                                                                                                                                                                                                             aTCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGt                                            <  1:415388/251‑1 (MQ=255)                                                                                                                                                                                                                   gCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCaaa                                        <  1:59626/249‑1 (MQ=255)                                                                                                                                                                                                                                 ggg‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCa                        <  1:52826/251‑1 (MQ=255)                                                                                                                                                                                                                                        tctctc‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGt                   <  1:85865/250‑1 (MQ=255)                                                                                                                                                                                                                                           ctc‑tACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTa                  <  1:548233/248‑1 (MQ=255)                                                                                                                                                                                                                                             c‑tACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATc                <  1:270751/248‑1 (MQ=255)                                                                                                                                                                                                                                                 cttact‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCAtt          >  1:495548/1‑251 (MQ=255)                                                                                                                                                                                                                                                           ttGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGc  <  1:289622/250‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          GGCATTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGG‑TATTCTCTC‑TACTTACT‑ATTTTGGCATTACTGGCACCTTTTGGGCTGTCAC‑GGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGG‑TTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGC  >  NC_000913/2008124‑2008620

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