Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,016,051	0	T	C	33.3%	109.8 / 15.4	66	D40G (GAT→GGT)	ysgA	putative dienelactone hydrolase
Reads supporting (aligned to +/- strand):  ref base T (29/13);  new base C (0/21);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 4.61e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGAT  >  NC_000913/4015802‑4016295                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        acGGCCAGCACCCGCGAGTCAGGCACTTTTTCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAAGTT‑CC‑TGCACGCCAATGACC‑ACTGGCAGGGGGCC‑‑A‑TCGCTTTGCTTGG‑GTC‑T‑GGCAGGGTAAGCAGGCATGTTATc                                                                                                                                                                                                                                                        <  2:388199/250‑1 (MQ=255)    gCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATcccc                                                                                                                                                                                                                                                     <  2:334797/251‑1 (MQ=255)      cAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCtt                                                                                                                                                                                                                                                   <  1:517608/251‑1 (MQ=255)       aGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTtg                                                                                                                                                                                                                                                  <  1:470782/251‑1 (MQ=255)           ccTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGaa                                                                                                                                                                                                                                              <  1:458774/251‑1 (MQ=255)                 aGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GAt                                                                                                                                                                                                                                       >  1:134204/1‑250 (MQ=255)                           ggcagattcaggACCGCTAAGCAACGTGTGGATATCGGCAAAATCATTCGGAGCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTGG‑GGC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTTTGAA‑TG‑GATAGATGTG‑AAGc                                                                                                                                                                                                                           <  2:427526/240‑1 (MQ=255)                            tttGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑GGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGc                                                                                                                                                                                                                           <  1:329148/250‑1 (MQ=255)                            tttGCTACCAGACAGCTAAGCATCGTGGGGAGATCGGCAAAATCATTCGGACCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCAGGCACGCCAAAAATGT‑CC‑TGCACTACAATGACC‑ACTTGCAGTGGGCC‑‑A‑GCGCGTTGCTTTG‑GGC‑T‑GGCATTGTAAGCAGTCATGTTAACCCCTTGTGAA‑TG‑GATAGAGGTG‑AAg                                                                                                                                                                                                                            <  2:270222/249‑1 (MQ=255)                                cTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCCGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAC‑GG‑GCTAGAGGTG‑AAGCCCGcc                                                                                                                                                                                                                      >  1:628737/1‑251 (MQ=255)                                 tACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGcc                                                                                                                                                                                                                      >  1:445849/1‑250 (MQ=255)                                  aCCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCcac                                                                                                                                                                                                                    >  1:184348/1‑251 (MQ=255)                                   ccAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCcac                                                                                                                                                                                                                    >  1:545954/1‑250 (MQ=255)                                             aagaaaCGGGGGGATACCGGAAAAATCATTCGGATCGCCTTCGCGAAAGAAAAGTTCAGGTGCGATACCAAGATCCCCCCCCAGCGCCAGACGGCGACAAATGCCGGGGATATGTTCAGGCACGCCAAAAATTT‑CC‑TGCACTACAAGGCCC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCTTGTTACCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGcgtc                                                                                                                                                                                                           <  1:242191/245‑1 (MQ=255)                                                cAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCg                                                                                                                                                                                                      >  1:120423/1‑251 (MQ=255)                                                 aaCGTGGGGATATCGGCCAAAGCATTCGGCTCGCCTTCGAGGAGGTAAAGTTCAGGTGCGCTAGCCAGCTACCCCTCCACCGCCCGACGGCGACAAATGTAGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCGCTACAATGCCC‑ACTGGCAGTGGGAC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTCTCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGGCCGCCACAATCGCGTCGTCCgg                                                                                                                                                                                                     >  2:409870/1‑251 (MQ=255)                                                       gggATATCGGCAAAATAATACGGCTGGCCTTCGCGCACGCAAAGTTCAGGTGCGATAGCTAGATCCCCCTCCAGCGCACGCCGGCGACCAATGTCGCGGCTATGTTCATGCACGCCACAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGATATACCCTTGTGAA‑GG‑GAGAGAGGTG‑AAGCCCGCAACCATCGCGGCGGCCggggacat                                                                                                                                                                                               >  2:338727/1‑247 (MQ=255)                                                       ggaatctcgaaaaaaTCATGCGGACCCCCTCCGCGAAAGTAAATTTAAGGGGCGAAAGCCAAACCCCCCCCCAGCCCCAGCCGGCGAAAAAGGTCGCGGATTGGTTCATGCCCGC‑AAAAATTT‑CCCTGCACAACATGGCCC‑ACGGGCGGGGGGCC‑‑A‑TCGTTTTG‑TTTG‑GGC‑TTGGAAGGGTAAGCAGGCTTGTTCCCCCCTGGGGAA‑GG‑GAAAGAGGGA‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTg                                                                                                                                                                                               <  2:665621/240‑1 (MQ=255)                                                                gCAAAATCATTCGGATAGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGt                                                                                                                                                                                      >  1:205163/1‑251 (MQ=255)                                                                  aaaaTCATTCGGACCCCCTTCGCGAAAGAAAAGTTCAGGTGCGAGAGCAAGATCCCCCCCCAGCGCCAGACGGCGAAAAATGTCGCGAAATTGTTCATGCACGCAAAAAATTC‑CC‑TGCACTACAA‑GACCCACTGGCGGTGGGCC‑‑A‑GCGTTTTGCTTGG‑GTC‑T‑GGTATGGAAAGCAGTCAGGTTACCCCCTGGGGAG‑GG‑GAAAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCg                                                                                                                                                                                    <  1:104967/251‑1 (MQ=255)                                                                     agcatccGGATCGCCTCCGCGGAAGAAAAGTTCAGGGGCGATACCCAGATCCCCCCCCAGCGCCAGACGGCGACAAATGCCGCGGATATGTTCATGCCCGCAAAAATTTT‑CC‑TGCACTAAAATGCCC‑CCGGGCAGTGTGCC‑‑A‑CCGTTTTGCTTGG‑GTG‑T‑GGCTTGGTAAGCAGGCTTGTTACCCCTTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAg                                                                                                                                                                                 <  1:354076/245‑1 (MQ=255)                                                                              aTCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAC‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCg                                                                                                                                                                        >  2:103295/1‑251 (MQ=255)                                                                                 gCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATg                                                                                                                                                                     >  1:143313/1‑251 (MQ=255)                                                                                        cgaaaataaagtgagggggCGATCGCCAGACACCCCTCCAGCGACAGCCGGCAAAAAATTCCGCTGTTATGTTCACGCACGACAAAAATTC‑CC‑TGCTCAAAAATGCCC‑AGTGGGGGGGGGCC‑‑C‑TCGCTTTGCTTTG‑GTCGT‑GG‑TTGATAGGCGGGCATTTGACCCCTTTGAGAT‑GG‑GATAGAGGGG‑AGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGGCGAAGCGAGAGGCGATGCAGCAgg                                                                                                                                                              <  2:637850/233‑1 (MQ=255)                                                                                           aagaaaggTTCAGGTGCGATAGCCAGATACCCCCCCGGCGCCAGAGGGCGAAAAATGTCGCGGATATGTTAATGCACGCAAAAATTTT‑CC‑TGCACTACAAGCCCC‑CCGGGCGGGGGGCC‑‑A‑CCGCTTTTTTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTACCCCCTTGTGAA‑GG‑GATAGAGGGG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGc                                                                                                                                                           <  1:648536/244‑1 (MQ=255)                                                                                           aaGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTtgtg                                                                                                                                                                                                                                                >  1:241937/1‑169 (MQ=255)                                                                                           aaGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAAGGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTtgtg                                                                                                                                                                                                                                                <  2:241937/169‑1 (MQ=255)                                                                                             gTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCa                                                                                                                                                          >  2:118505/1‑250 (MQ=255)                                                                                                  gTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                                                                                                                                     >  2:553739/1‑250 (MQ=255)                                                                                                         ggCGATAGCCAGACACCCCCCCAGCGCAGGAGGGCGACAAAGGTCGAGAATATGTCCATGCCCGC‑AAAAATTT‑CCCTGCATTACAATGGCC‑A‑‑GGGACTTGGCCCTG‑TCTTTTTGCTTTG‑GTC‑G‑GGCTGGGAAAGCAGTCATGTTGGCCGCTTGT‑AA‑GGTGTAAGGGGTT‑GAGCCCGCCGCAATCGCTTCGTCCGGGGTCTGTGCTTTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgtt                                                                                                                                             <  1:297609/250‑1 (MQ=255)                                                                                                             aTAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgt                                                                                                                                         >  1:368331/1‑251 (MQ=255)                                                                                                                    gATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                                                                                                                  >  2:10693/1‑251 (MQ=255)                                                                                                                    gATACCCCTCCAGCGCCAGACGGCGACAAATGCCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGACCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                                                                                                                  >  2:36583/1‑251 (MQ=255)                                                                                                                       ccccccccacccccAGCCGGCGACAAAGTTCGCGGATATGTCCATCCCCGCAAAAAATTC‑CC‑TCCACTCCAATAACC‑CCTGGCAGGGGCCC‑‑A‑TCTTTTTGTTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGGAAA‑GG‑AATAG‑GGTGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                                                                                                                               <  1:463588/239‑1 (MQ=255)                                                                                                                       aCCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                                                                                                                               >  1:330895/1‑251 (MQ=255)                                                                                                                        ccccTCCAGCGCCAGACGGCGAAAAATGTCGCGTATATGTTCATGCACGCAAAAAATTT‑CC‑TGCACTACAATGCCC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGAAAGCAGGCATGTTACCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATgg                                                                                                                              <  1:432481/251‑1 (MQ=255)                                                                                                                                ccgcCAGACGGCGAAAAATGTCGCGGATATGTTCATGCCCGC‑AAAAATTT‑CCCTCCACTACATGGCCC‑TCTGGCAGTGGCCC‑‑T‑TCTTTTTGTTTTG‑GTC‑T‑GGCATGGAAAGCAGGCTTGTTACCCCCTGGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTTTTGTTGCCATGGTATTCTcc                                                                                                                      <  1:259463/250‑1 (MQ=255)                                                                                                                                   cccggcggagACAAATGGCGCGAATGTGTTAAGGCACGC‑AAAAATTT‑CCCTGCACAACAATGCCC‑CCGGGC‑GTGGGCC‑‑CCTCGTTTTGTTTTG‑GCC‑G‑GGAATGGAAGGCGGACTTGTTCCCCCTTTGTGAA‑GA‑AAGAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACTATGGTCGAAGCGTGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTGGCCATGGTATTCTCCGTa                                                                                                                   <  2:422433/242‑1 (MQ=255)                                                                                                                                            gACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata                                                                                                           >  2:71602/1‑250 (MQ=255)                                                                                                                                                   gTCGCGGATTTGTTCATGCACCAAAAAAATTT‑CC‑TCCACTACAAGGCCC‑ACTGGCATGGGGCG‑‑A‑GCGTTTTGCTTTG‑GTC‑T‑GGCGTGGTAAGCAGGCATTTTACCCCCTTGTGAG‑GG‑GATAGTGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTAGAAAAATGtt                                                                                                   <  1:209497/251‑1 (MQ=255)                                                                                                                                                        ggATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgcgc                                                                                               >  2:177464/1‑250 (MQ=255)                                                                                                                                                              ttccATGCGCGCCAAAATTTG‑CC‑TGCACGACAATGCCC‑GCGGGCAGTGGGCC‑‑T‑GCGCTTTGCTTGG‑GTT‑T‑GGCTTGAGAAGCGGGCATGTTACCCCCTTGTGAA‑GG‑GAGAGGGGGG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGtt                                                                                         <  1:668976/247‑1 (MQ=255)                                                                                                                                                                gcaggcccgaaaaaaaTTC‑CC‑TGCCCAAAAATGCCC‑ACGGGCGGGGGCCC‑‑A‑TCGCTTTGTTGTG‑GTC‑G‑GGCGGGGAAAGCAGGCATGTTACCCCCTTGGGAA‑GA‑GAAGGGGGGG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa                                                                                       <  2:239969/239‑1 (MQ=255)                                                                                                                                                                 caggcccgcaaaaaaTTC‑CC‑TGCACAAAAATGCCC‑CCTGGCGGTGGGCC‑‑G‑TCGCTTTGTTTTG‑GTC‑T‑GGCAGGGAAAGCAGGCATTTTACCCCCTTGTGAA‑GG‑TATAGAGGGG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAAc                                                                                      <  1:428074/240‑1 (MQ=255)                                                                                                                                                                    ccACGCAAAAATTTT‑CC‑TGCACTAAAATGACC‑ACGGGCGGTGGCCC‑‑A‑T‑GGTTTGCTTTGTGTC‑T‑GGCATGGTAAGCAGGATTTTTTTCCTCTTGTAAA‑GG‑GAGAGAGTTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtat                                                                                   <  1:488031/249‑1 (MQ=255)                                                                                                                                                                     gaacccAAAAATTTTCC‑TGCCATACAAGGACA‑CCGGGCAGGGGGCC‑‑A‑TCGTTTGGTTTTG‑GTC‑G‑GGAAGGGAAAGAAGCCATTTTACCCCCTTGTAAA‑GG‑GATAGAGGTG‑AACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtata                                                                                  <  1:256611/247‑1 (MQ=255)                                                                                                                                                                             aatttt‑CC‑GGCACTACAAGGCCC‑ACGGGCGGGGGGCC‑‑T‑TCGTTTTGTTTTG‑GTG‑T‑GGCATGGTAAGCAGGCATGTTACCCCTTTGTGAG‑GG‑GATAGGGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATg                                                                         <  1:188474/248‑1 (MQ=255)                                                                                                                                                                              aaTTT‑CCCTGCACTACAAGCACC‑ACGGGCAGTGGCCC‑‑A‑CCGTTTTGCTTGG‑GTG‑T‑GGCAGGGAAAGCGGGCAGGTTACCCCCTTGGGAG‑GG‑GAGAGAGGGA‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCGCCGTACCCGTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATg                                                                         <  1:71602/251‑1 (MQ=255)                                                                                                                                                                              aaTTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGc                                                                        >  2:492924/1‑251 (MQ=255)                                                                                                                                                                                    cc‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaaca                                                                    >  2:179102/1‑250 (MQ=255)                                                                                                                                                                                                  gccc‑ACGGGCAGGGGCCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGAAAGCGGGCATTTTA‑CCCCTTGGGAAAGG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTg                                                      <  1:118505/249‑1 (MQ=255)                                                                                                                                                                                                       ccTGGCAGTGGGCC‑‑A‑TCGTTTTTCTTGG‑GTC‑G‑GGCAGGGTAAGCAGGCATGTTCCCCCGTGGTGAA‑GG‑GAAGGAGGTA‑AGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc                                                  <  1:382133/250‑1 (MQ=255)                                                                                                                                                                                                             aGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                                                                       >  1:318656/1‑192 (MQ=255)                                                                                                                                                                                                              atgcgac‑‑A‑TCGCTGTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCAt                                           >  2:533524/7‑251 (MQ=255)                                                                                                                                                                                                                     ccctccccttggtttgg‑gTC‑T‑GGCATGGAAGGCAGCCTTTTTACCCCCTGGTGAG‑GG‑GATAGGGGTG‑AGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATGGCCTGAATCGGCTCATCttt                                       <  1:148636/235‑1 (MQ=255)                                                                                                                                                                                                                         tCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTCTCCCCTTGGGAA‑GG‑AAAAGGGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTAt                                    <  2:368331/250‑1 (MQ=255)                                                                                                                                                                                                                                 tttgg‑gTC‑T‑GGCAGGGGAAGCAGGCATGTTATCCCTTTGGGAG‑GG‑GAGAGGGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCt                            <  1:139570/246‑1 (MQ=255)                                                                                                                                                                                                                                  tttg‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGc                          >  1:44006/1‑251 (MQ=255)                                                                                                                                                                                                                                  tttg‑GTC‑G‑GGCATGGAAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑AATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTg                           <  1:71097/250‑1 (MQ=255)                                                                                                                                                                                                                                        tC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGGACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTACTGCAGAAt                     >  2:561122/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ggtaaaGCAGGCATGTTACCCCCTTGGGAA‑GG‑AATGGGGGGG‑AGCCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGt              <  1:103295/247‑1 (MQ=255)                                                                                                                                                                                                                                                 gggaaagcagtgttGTTACCCCTTTGTGAA‑GG‑TATAGGGGGG‑AATCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGt              <  1:530347/238‑1 (MQ=255)                                                                                                                                                                                                                                                       gaaGGCATGTTCCCCCCTTGTGAG‑GG‑AATAGAGGTA‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTTATTTa        <  1:492924/249‑1 (MQ=255)                                                                                                                                                                                                                                                         aGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTAtg      >  1:178031/1‑251 (MQ=255)                                                                                                                                                                                                                                                             aTGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGAt  >  1:573443/1‑251 (MQ=255)                                                                                                                                                                                                                                                             aTGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTACGGGGTCTGAACCATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAAACCCGATTGTTGTGATGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCCTGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTCTGCAGTCCTGCAGAATGAAGGGTGCTTTATGTGAt  >  1:62748/1‑251 (MQ=255)                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        TCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTT‑CC‑TGCACTACAATGACC‑ACTGGCAGTGGGCC‑‑A‑TCGCTTTGCTTTG‑GTC‑T‑GGCATGGTAAGCAGGCATGTTATCCCCTTGTGAA‑GG‑GATAGAGGTG‑AAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGAT  >  NC_000913/4015802‑4016295

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