Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,644	0	C	T	24.5%	102.7 / 14.8	49	intergenic (‑16/‑184)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base C (8/29);  new base T (12/0);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.37e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGT  >  NC_000913/4223397‑4223893                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTg                                                                                                                                                                                                                                                               <  1:14011/250‑1 (MQ=255)     aaCCTGACGCTCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCATTTCCGTGATTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCg                                                                                                                                                                                                                                                          <  2:270833/251‑1 (MQ=255)              cACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                                                 <  1:149813/251‑1 (MQ=255)               aCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTTGCGGG‑GATGGGGGCGCCCATGAC‑CATCTCTTTTTTTTTTATCGTGGGAAttat                                                                                                                                                                                                                                                >  1:175256/1‑248 (MQ=255)               aCGAAACCCTGCTGTTGCATCGAGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCCGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCACCGGCGGGTTTTTTTGCGCGGTTCGCGGG‑GGTGGGGGGGGCCCTGAC‑AGTCTCTTTTTTTTGTTTCGTggtacacat                                                                                                                                                                                                                                                >  1:355811/1‑244 (MQ=255)                    aCCCTGCTGTTGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACTGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCGCTTTTTTCGGTATCGTGGAATtcattttca                                                                                                                                                                                                                                           <  2:356916/251‑1 (MQ=255)                         gctgttgtagtgTGGTTAGCAGGCGGTGGGTCGTGGAAATGGGTAACCCGGCTTGTAGCGCCAGTGACGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGTCGCGTCTCGTGGG‑AATGGGTGAGACCAGGAC‑AGTCTCCTTTTTCTGTAACGTGGAAAtcatttacattt                                                                                                                                                                                                                                        <  2:118602/238‑1 (MQ=255)                         gctgTTGCTTCGTGGTTAGCAGGCGGTGGGTCGTGTAATTGGGTACCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt                                                                                                                                                                                                                                      <  2:297489/251‑1 (MQ=255)                              tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCCCGTTTCGCGGG‑GATGGGGGGGGCCCAGAC‑AGTCTCTTTTTTTTGGTTTGTGGGAATCatttttatttttttttt                                                                                                                                                                                                                                 >  1:175860/1‑246 (MQ=255)                               gCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTGCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGACACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTTGCAACGGCGGTTTGTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGtt                                                                                                                                                                                                                                <  2:304564/251‑1 (MQ=255)                                cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTGCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTa                                                                                                                                                                                                                               <  1:394611/251‑1 (MQ=255)                                         tAGCAGGCGGTGGGTCGGGGAATTGGGTAGCCCGTCGGGTTGCGCCAGTTACGTGAGTGCCACACTGCCAGTGGATGCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTGTTTGTGCTGTATCGTGGAAATCATTTTCATGTTTATTGTTAGCTAATG‑c                                                                                                                                                                                                                      <  2:260527/250‑1 (MQ=255)                                          aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑Ca                                                                                                                                                                                                                     <  2:147938/250‑1 (MQ=255)                                                 ggTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTac                                                                                                                                                                                                             <  2:243657/251‑1 (MQ=255)                                                   ggggTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTactg                                                                                                                                                                                                           <  2:353247/250‑1 (MQ=255)                                                      gTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGTTTTGCGGG‑AATGGGGGGGGCCATGAC‑AGGCTCCTTTTTCTTTTTTGTGGGAATCATTTTTCTTTTTTTTTTTAGCTAATG‑CAAAAGTTactgaac                                                                                                                                                                                                        >  1:329643/1‑251 (MQ=255)                                                                           gCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTg                                                                                                                                                                                <  1:71415/251‑1 (MQ=255)                                                                                          tCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGGCGCGGTTCGCGGG‑GATGGGTGCGGCCATGGC‑AGTCTCTTTTTTTTGTTTCGTGGGAATCATTTTTATTTTTTTTTTTTGCTTATG‑CAATAGTTAATGAAA‑TGAT‑TCGCTGAG‑TTAATGTTTAAAAAAACTCAT‑Gtt                                                                                                                                                                >  1:160354/1‑251 (MQ=255)                                                                                          tCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTCGCGGG‑GATGGGGGCGGCCATGAC‑AGACTCTTTTTTTTGTTTCGTGGGAATCCTATTTATTTTTTTTTTTATCTTATA‑CCATAGTTTCTGGAC‑TGAT‑ACGATGAG‑tttattttttacaaatatcttttt                                                                                                                                                                 >  2:12171/1‑229 (MQ=255)                                                                                                     cacaCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCTCGGG‑AAAGGGGGGGACCACTAC‑CGGCTCTTTTTTTTGTTTTGTGGGAATTAATTTCATTTTTTTTGTTTGCTTATG‑CTATAGATTCTGAAC‑TGAT‑ACTCTGAG‑GTAATTTTGAACAAATCTCAT‑TTTgtgggggggtc                                                                                                                                                     >  2:92379/1‑241 (MQ=255)                                                                                                        acTGCCAATGGATTCGGCAATTCACTCCAGTTATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCTCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAT‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCgat                                                                                                                                                  <  2:326759/251‑3 (MQ=255)                                                                                                            ccATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTTCGCGTTTCGCGGG‑GGTGGGGGGGCCCATGAC‑AGTCTCTTTTTTTTGTATAGTGGGAATCCTTTTTATTTTTATTTTTAGTGAATT‑CACTAGTTACTGAAC‑TGAT‑CCGAAGAA‑TTAATATTGTACAAATCCCAT‑TTTGCTTGGTGGGCGCctttt                                                                                                                                              >  2:434265/1‑250 (MQ=255)                                                                                                              aTTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTTCTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTAcc                                                                                                                                            <  2:415259/251‑1 (MQ=255)                                                                                                                tGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGG‑AATGGGGGCGACCCTGAC‑AGTATCTTTTTTTTGTATCGTGGAAATCATTTTCATTTTTTTTGTTAGCTACTG‑CAATAGTTACTGAAC‑TGAT‑TCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGGGGGGGTCGCTTTTAcccc                                                                                                                                          >  1:270371/1‑249 (MQ=255)                                                                                                                 ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGAATGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGCTCCCCAGG‑GATGGGGGCGGACATAAC‑AGTCTCTTTTTTTTGTATCGTGGAAATCACTTTCCATTTTATTTTTAGCTAATGTCAA‑AGTTATTGAAC‑TTAT‑CCAGTCAG‑GTAATGTTGACCAAATATCAA‑ATTGCGTGGTGGTCGCTTTTTCcacc                                                                                                                                         >  1:132445/1‑250 (MQ=255)                                                                                                                        gCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑ACAGGGTGGGACCATTAC‑AGTCTCTTTTTTCTTTTTCGTGGAAAATATTTTCATTTTTTTTTTTTGCTTATG‑CAATAGTTTCTTAAC‑TGAT‑TCGATGGG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGGGGTCGCTTTTTCcacaaaagcgg                                                                                                                                  >  1:8446/1‑244 (MQ=255)                                                                                                                            tcaaCTCCAGTAATTTCAGGCCACGAGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGATTTTTGTCGCGGTTAGCGGG‑AATTGGTGGGGCCCAGAC‑AGTCTTTTTTTTCTTTTTCGTGGGAATAATTTTCCTTTTTTTTTTTAAGTAATA‑CACTAATTTATGGAA‑TGATTCCAGTGAG‑TTA‑TGTTGAAAAAATCTCAT‑TTTGCGTGGGGGTCGCtttttcccacaacgcgtttat                                                                                                                              >  2:436663/4‑234 (MQ=255)                                                                                                                               aCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGGCCCTGAC‑AATCTCTTTTTTTTGTATCGTGGGAATCCTTTTTTTTTTTATTTTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGGTGAG‑TTAATGTTTAACCAATCTCAT‑GTTGTGTGGGGGTCGCTTTTACCACAGAAGCG‑GTTATGcc                                                                                                                          >  1:203834/1‑251 (MQ=255)                                                                                                                                  ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGAt                                                                                                                                      <  1:167381/237‑1 (MQ=255)                                                                                                                                  ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTAGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGCTTCGCGGGGAATGG‑TGCGACCCAGACCAG‑CTCTTTTTTTTGTATGTTGGGGATCATTTTTTTTTTTTTTTTTTTATAATG‑CAATAGTTACTGAACATG‑T‑TCGAATAG‑TTAATGTTGTAAAAATCTCATTGTTCTGTG‑TGGGCGCTTTTACcacatatagtacggaatagcg                                                                                                                        >  2:167381/1‑234 (MQ=255)                                                                                                                                      tAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCTGGTTTTTTGCCGCGGTTCGCGGG‑AATGGGTGCGGCCCTGGC‑CGTCTCTTTTTTTTGTATCGTGGAAATCAATTTTATTTGTATTTTTTGCGATTG‑CCAACGTTTATGAAA‑TGTT‑CCGATGAGGTCATTGTTGA‑CAAAACTCAT‑TTTGCGTGGTGGTCGGTTTTACCACACATGGGGTTTATGC‑AGTATgg                                                                                                                   >  2:120038/1‑251 (MQ=255)                                                                                                                                            cAGGCCACGCGGTAAAGACTGGACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCGAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCGTTTACCACAGATGCG‑TTTATGCCAGTATGGTttgt                                                                                                              <  2:175256/250‑1 (MQ=255)                                                                                                                                             aGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTttgttg                                                                                                            <  1:12171/251‑1 (MQ=255)                                                                                                                                                        tAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTttatt                                                                                                  <  2:160354/250‑1 (MQ=255)                                                                                                                                                                gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTggg                                                                                         <  1:389564/251‑1 (MQ=255)                                                                                                                                                                gAACCTGTCCAGGCGCTGGTGCGTGGGCAACGGCGGGTTTTCTGCCGCGTGTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTg                                                                                           <  2:12160/249‑1 (MQ=255)                                                                                                                                                                  aCCTGTCCAGGCGTTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATGTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATTTGGGtt                                                                                       <  2:433701/251‑1 (MQ=255)                                                                                                                                                                                                 cggGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaaa                                                        <  1:404128/251‑1 (MQ=255)                                                                                                                                                                                                    gttttctttctcgTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCGTTTTTCTGTATCGTGGAAATCATTTTCATGTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAAGCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGCGCAGCAAAGTg                                                     <  2:248515/240‑1 (MQ=255)                                                                                                                                                                                                        tCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATGTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAG‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaa                                                  <  2:232083/250‑1 (MQ=255)                                                                                                                                                                                                                  ttCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                        <  1:333044/250‑1 (MQ=255)                                                                                                                                                                                                                     gcgGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTgcgc                                     >  1:77310/1‑250 (MQ=255)                                                                                                                                                                                                                            aTGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCCTTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa                               >  1:432936/1‑250 (MQ=255)                                                                                                                                                                                                                                gTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTGATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                          <  1:436663/251‑1 (MQ=255)                                                                                                                                                                                                                                       cATGAC‑AGTCTCTTTTTTCTGTATCGCGGAAATCATTTTCCTTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTc                   >  2:401398/1‑251 (MQ=255)                                                                                                                                                                                                                                           aC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt               <  2:8446/251‑1 (MQ=255)                                                                                                                                                                                                                                            acagtatctggtgtCTGTATCGTGGAAATCATTTTCATTTTTAGTGTTAGCTAATG‑CAAGAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGGTGAACAAATCTCAT‑GTTGCGTGGTGGGCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCAGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTTCGTGCGCAGTTAAATGACCGTATTCTGGt               <  2:355811/238‑1 (MQ=255)                                                                                                                                                                                                                                                  tctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTgg          >  1:213914/1‑250 (MQ=255)                                                                                                                                                                                                                                                         tCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGAGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt  >  1:413439/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAAC‑TGAT‑CCGATGAG‑TTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACAGATGCG‑TTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGT  >  NC_000913/4223397‑4223893

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