Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,658,163	0	C	T	29.6%	70.0 / 17.6	56	S32F (TCT→TTT)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base C (17/21);  new base T (16/0);  total (35/21)
Fisher's exact test for biased strand distribution p-value = 1.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGA  >  NC_000913/1657916‑1658405                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             aCATACAACATAATCCCACCTTATTACGCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCttt                                                                                                                                                                                                                                                           <  2:538733/251‑1 (MQ=255)     acaacaTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATTGGTAGACTCCCTCTATTGTTAGCGCGCGAAATATTCAATATATAAACTTTTATATAACGATAAAGCACAGGGAGTGAGTTATGTCCAAAAATGACCGAATGGTGGGCATCATCCGCAGAACACTCGTTAAATCTCCCGCGATAGGTTCTCTGGCGCTGGCTGCCGGCGGTTTTTCTTTG‑cg                                                                                                                                                                                                                                                       <  2:194740/250‑2 (MQ=255)       aacaTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCATCTATTGTTAGCACGCTAAATATTCAATATCTAAACTTTTATATAACGATAACGAACAGGGAGTGAGTTATGTCAAAAAATGAACGTATGGTAGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGTTAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTAttttctgtt                                                                                                                                                                                                                                                     >  2:679187/1‑245 (MQ=255)        acaTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGttt                                                                                                                                                                                                                                                   <  2:553859/251‑1 (MQ=255)            aaTCCCACCTTATTACTCATACCCTTCTATTGATATGGCTTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGCTAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGGCCGTTTAcc                                                                                                                                                                                                                                                >  2:607293/1‑251 (MQ=255)                      tattaCTCATACCCTTCTATTGATATGGCTTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGACATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGCTTTTCttttcctttttacccgccccatgc                                                                                                                                                                                                                                      >  2:182806/1‑230 (MQ=255)                       attaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATgc                                                                                                                                                                                                                                    <  1:171327/250‑1 (MQ=255)                             cATACCCTTCTATTGATATGGATTAATAATTCTTAACACAAAATGGGTAGACTCCCTCTATTGTTAGCGAGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCTGGCGGGTTTTATTTT‑CCGTTTACCCCTG‑GCAATgca‑acagcc                                                                                                                                                                                                                            >  2:693168/1‑250 (MQ=255)                                       tATTGATATGGATTATTAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTacaa                                                                                                                                                                                                                   <  1:54910/250‑1 (MQ=255)                                       tATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTACCC‑CGCGCCATGCA‑GCAGCAGCGGTgcaa                                                                                                                                                                                                                   >  2:54910/1‑250 (MQ=255)                                              atGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑TCGTTTAACC‑TGCGCAAAGCA‑GCAGCAGCGGGACAACAgggccc                                                                                                                                                                                                           >  1:190267/1‑247 (MQ=255)                                              atGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGccc                                                                                                                                                                                                            <  1:607293/250‑1 (MQ=255)                                                        aatTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATAGTCAATATATAAACTGTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACCCTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTGG‑CCGTGTAACC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTg                                                                                                                                                                                               <  2:151052/251‑1 (MQ=255)                                                                    aaaaTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGAGG‑GCTGGGGTGcc                                                                                                                                                                                  >  1:688494/1‑251 (MQ=255)                                                                    aaaaTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAAAAGGCCCGCG‑AA‑AAAgagggacgggggggc                                                                                                                                                                                   >  1:508466/1‑237 (MQ=255)                                                                       aTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTT‑CCGTTTTCCC‑TGCGCCATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAAGGG‑TCTGGGGGGCC‑TGt                                                                                                                                                                              >  1:1380/1‑251 (MQ=255)                                                                                    ctctATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑CGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCCCG‑GA‑AAAGAGG‑TCTGGGGGGGC‑TGTTTCGTccaccgtg                                                                                                                                                                 >  1:449644/1‑244 (MQ=255)                                                                                       tATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTACCC‑TGCGCCAAGCA‑GCAGCAGCGGGACAACAAGCCCGCG‑AA‑ACAGTGG‑TCTGGGGCGCC‑TGGTCCGTCaaacgtggt                                                                                                                                                               >  1:742120/1‑243 (MQ=255)                                                                                          tGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGc                                                                                                                                                           <  1:743883/250‑1 (MQ=255)                                                                                             tAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTTCCC‑TGCGCAAAGCA‑GCAGCACGGGGACAAAAGGCCCGCG‑AA‑AAAGTGGGACTGGGG‑GCC‑TGTTACGGCAACTGTGG‑GAGccccc                                                                                                                                                       >  1:74341/1‑248 (MQ=255)                                                                                                 gcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTACCC‑TGCGCCATGCA‑GCAGCAGCGGTACAACAGGGCCGCG‑AA‑AAAGGGG‑TCTTGGGGGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCtgtgg                                                                                                                                                   >  1:215708/1‑250 (MQ=255)                                                                                                     tAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTATCC‑TGCGCAATGCA‑GGAGCCGCGGTACAACATGCCCCCG‑AACAAAGGGG‑CCTGGGG‑GCC‑CGATTCGTCAACTGTGG‑TAGCCGCtggtcaatt                                                                                                                                               >  2:274173/1‑244 (MQ=255)                                                                                                         tatTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGAGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑ACGCTTTCCC‑TGCGCAATGGA‑GCAGAAGCGGGACAATACGCCCGCG‑AA‑AAAGTGG‑GCCGGGGTGGC‑TGTTCCGTCAAACGTGG‑GAGCCGCTGTG‑CACTTccgc                                                                                                                                          >  1:660441/1‑248 (MQ=255)                                                                                                             cAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGCTTACCC‑TGCGCAATGCA‑GCCGCAGCGGTACAACAGGCCCGCG‑GA‑AAAGGGG‑GCTGGGGGGCC‑TGTTCCGTCAACTGTGG‑GAGCCGCTGTG‑CACTTCGTCTAc                                                                                                                                       >  1:636657/1‑250 (MQ=255)                                                                                                                        cTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGAGAGGTTCTCTGGCGCTGGCTGCAGGCGGATTTTTTTTG‑CCGGTTACCC‑CGCGCAAAGCA‑AGAGCAGCGGCACAAACGAGCCGCG‑AA‑AAAGTGG‑GCTGGGGGGGC‑CGTTCCGGCCACTGTGG‑GAGCCGCTGTG‑cccagcggctaccatttaaag                                                                                                                              >  1:679925/1‑228 (MQ=255)                                                                                                                             atataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAAAAGGCCCGCG‑AA‑AAAATGG‑TCTGGGGGGCC‑TGTTTCGTCAACTGTGG‑TAGCCCCTGTG‑CACTTCGTCTACCTGTTAAAGATAATgg                                                                                                                       >  1:232656/1‑249 (MQ=255)                                                                                                                                 aaCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCC‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGAGG‑GCTGGGGGGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACCTCGTCTACATGTTAAAGATAATGAAGTg                                                                                                                   >  1:408732/1‑250 (MQ=255)                                                                                                                                 aaCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGAGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGGTTTTTTTTG‑ACGTTTAACC‑CGCGCAATGCA‑GGAGCAGCGGGACCACAAGCCCGCG‑TA‑AAAGTGG‑GCTGGGGGGCC‑CGTTCAGACAAATGTGT‑GAGACGCCGTG‑CACCTCGTGTTCATGtttacagtaatgaaatgt                                                                                                                  >  2:586006/1‑235 (MQ=255)                                                                                                                                   cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑cc                                                                                                               <  1:326632/251‑1 (MQ=255)                                                                                                                                    gATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTATGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑ACGTTTACCC‑CGCGCAATGCA‑GCCGACGCGGTACAACAAGCCCGCG‑AA‑AAAATGG‑GCGGGGGCGCC‑GGATCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTAAATGTTAAAGATAATGAAGTGA‑ccc                                                                                                              >  1:619852/1‑250 (MQ=255)                                                                                                                                     aTAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTg                                                                                                             <  1:70170/251‑1 (MQ=255)                                                                                                                                       aaaGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGATTTTTTTTG‑CCGTTTACCC‑CGCGCAATCCG‑GCAGCAGCGGTAAAACAGGCCCGCG‑CC‑AAAGGGG‑TCTGGGGGGCC‑CGTTCCCTCCACTGTTG‑TAGCCGCTGTG‑CACTTCGTCTACCTGTTAAAGATAATGAAGTGA‑CCAggg                                                                                                           >  1:30673/1‑251 (MQ=255)                                                                                                                                         aGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAAAAGGCCCGCG‑AA‑AAAGTGG‑GCTGGGGGGCC‑CGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑GACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTggggg                                                                                                         >  1:72347/1‑249 (MQ=255)                                                                                                                                          gAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACGGTGG‑TAGCCg                                                                                                                                                         <  2:664782/204‑1 (MQ=255)                                                                                                                                               gggAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTCCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACt                                                                                                   <  2:687498/251‑1 (MQ=255)                                                                                                                                                      agtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCCATGCA‑GCAGCCGCGGGACAACAGGCCCGCG‑CA‑AAAGGGG‑TCTGGGGTGCC‑TGTTTCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTCCCTGTTAAAGATAATGCAGTGA‑CCCGGGTGGGAACTGACAAt                                                                                             >  1:264116/1‑250 (MQ=255)                                                                                                                                                         tATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTg                                                                                                         <  1:680280/235‑1 (MQ=255)                                                                                                                                                                 aaaaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAAAAGGCCCGCG‑GA‑AAAGGGG‑TCTGGGGGGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTTAAGATAATGAAGTGA‑CCTGGGGGGAAACTGACAATACCGGCAGCg                                                                                   >  1:330052/1‑249 (MQ=255)                                                                                                                                                                    aTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTT‑CCGTTTACCC‑TGCGCAATGGA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGCGG‑TCTGGGGGGCC‑TGTTCCGTCAAAACTGG‑TAGCCGCTGTG‑CGCTTCGTCTACATGTTAAAGATAATGAAGTGA‑CATGGGTGGAAACTGCCAATACCCGCAGCgctg                                                                                >  1:504338/1‑246 (MQ=255)                                                                                                                                                                       aaCGAATGGTGGGCTTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTAc                                                                           <  2:190267/251‑1 (MQ=255)                                                                                                                                                                       aaCGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTACCC‑CGCGCAATGCA‑GCAACAGCGGGACAACAGGCCCCCGGAA‑AAAG‑GG‑GCTGGGGGGGC‑TGTTCCGGCAACTGTGG‑GAGCCGCCGTG‑CACTTTGGCTACATGTTAAAGATAATGAAGGGA‑CCTGGGGGGAAAATGACAATACCGGCAGCGagggggac                                                                           >  1:122226/1‑244 (MQ=255)                                                                                                                                                                          gAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑CCGTTTTCCC‑CGCGCCATGCC‑GCAGCAGCGGGAAAAAAGGACCGCG‑CA‑AAAGGGG‑TCTTGGGTGCC‑CGTTCCGTCCACTGTTG‑GCGCCGCTGTGGCACTCCG‑CTACATGTTAAAGATAATGAAGGGA‑CCTGGGTGGNAACCTACCATAACGGCAGCGaaggggaaggg                                                                        >  2:386680/1‑241 (MQ=255)                                                                                                                                                                                  ggCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGCTCTCTGGCGCTGGCTGCAGGCGGGTTTTCTTTG‑GCGTTTGCAA‑TGCGCAACGCAAGC‑GAAACGGGCCCACCAGCCCGCC‑AA‑AAAAGGG‑ACTGGGGGGCCATGT‑CCGGCAACGGGGG‑GAGCCGCTGTG‑GAATTCGTCCCCATGTTTAAGATAAAGGAGTGG‑ACTGGGTGGGAACTGACAATACAGGCacaggagagtacggacacccac                                                                 >  2:272976/1‑229 (MQ=255)                                                                                                                                                                                   gCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAg                                                               <  1:182806/251‑1 (MQ=255)                                                                                                                                                                                        aGCCGCCGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCCGCAACGGTGCAACAGGCCCGCG‑AA‑CAAGGGG‑TCTGGGGGGCC‑TGTTCCGTCAAC‑GTGGGTAGCCGCCGTG‑CACTTCCTCTACATGTTAAAGATAATGAAG‑GAGCCTGGGTGGAAACTGACAATACCGGCAGCGAGGAGTGCGGCACACATCAGGtg                                                            >  2:200539/1‑248 (MQ=255)                                                                                                                                                                                              aGAACACTCGTTAAATCTACCGCGATAGGTTCTCTAGCGCTGGCTGCAGGCGGTTTTTTTTTT‑CCGTTTACCC‑CTTGCAATGGA‑GCAGACGTGGTACAACAGGCCCCCG‑CA‑AAAACGG‑CCTGGGGTGCC‑TGCTCCGTCAACTGTAG‑TTGCCAGTGTG‑CACTTCGGCTCCATGTTAAAGATAATGAAATGA‑ACTGGGGGGGAACAGACAATACCGGCAGCGAAGAGTAGGGCCACCCTCAGGTACGCGCCTg                                                    >  1:204931/1‑251 (MQ=255)                                                                                                                                                                                                   ccTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTTCCC‑TGCGCAAGGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACGGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTTA‑CCTGGGTGGAAACTGACAATTCCGGCAGCGATGAGTACTGCAACCATCAGGTACTCGCCTGTTTgc                                               <  2:504338/250‑1 (MQ=255)                                                                                                                                                                                                                aCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTT‑CCGTTTACCC‑TCCGCAATGCA‑GCCGCAGCGGTAAAACAGGCCCGCG‑AA‑AAAGGGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CCCTTCGTCTACATGTTTAAGGTAATGAAGTGT‑CCTGGGGGGAAACTGGCAAAACCGGGAGCGGTGAGTACGGCAAACACCAGGGACGCCCCCGTTTTCGCGGTCGCTcccc                                  >  1:408815/1‑249 (MQ=255)                                                                                                                                                                                                                aCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTcc                                    <  1:274173/249‑1 (MQ=255)                                                                                                                                                                                                                     gagaGGTTCTCTGGCGCTGGCGGCGGGCGGTTTGTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCTCGCG‑AA‑AAATTGG‑TCTGGTGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTTTG‑CACTTCGTCTACATGTTAAAGATAATGAAGGGA‑CCTGGGTGGAAACTGCCAATACCGGCAGCGATGAGGACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgcc                             <  2:660441/248‑1 (MQ=255)                                                                                                                                                                                                                     gATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgcc                             <  2:122226/251‑1 (MQ=255)                                                                                                                                                                                                                     gATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGGTTACCC‑TGCGCCATGCA‑GCAGGAGCGGTACAAAAAGCCCGCG‑GA‑AAAATGG‑GCTGGGGGGCC‑CGGTCCGTCAACTGTTG‑TAACCGCTGGG‑CACCTCGTCTACATGGTTAAGATAATGAAGGGG‑ACTGGGGGGGAACCGACAATACCGGGAGGGACGAGGACAGCAACCATTCGGCACGCGCCTGTTTGGCCGGTCTCTTCATccgcc                             >  1:456395/1‑251 (MQ=255)                                                                                                                                                                                                                                    gctggctgCAGGCGGTTTTTCTTTG‑CCTTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAatcatc              <  1:679187/251‑1 (MQ=255)                                                                                                                                                                                                                                       ggctgCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGGTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATcc             <  2:330052/249‑1 (MQ=255)                                                                                                                                                                                                                                                 ggTTTTTCTTAG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCAGCG‑AA‑AAAGGGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATCCCGGCAGAGATGAGTACGGCAACCCACAGGTACGCGCCTGGTTGCGCGGTCGCACCATCCGCCGGCGTATTAATCATCCCGATTGCtgta  >  2:585676/1‑247 (MQ=255)                                                                                                                                                                                                                                                  gTTTTTCTTTG‑CCGTTTACCC‑TGAGCAATGCA‑GCAGCAGCGGTATAACAGGCCCGAT‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGGGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACACTCAGGTACGCGCCTGTTTGCGCGGGCGCTACATCCACCGGCGTATTAATCATCCCGATCGCTTGa  >  2:258436/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             ACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCC‑TGCGCAATGCA‑GCAGCAGCGGTACAACAGGCCCGCG‑AA‑AAAGTGG‑TCTGGGGTGCC‑TGTTCCGTCAACTGTGG‑TAGCCGCTGTG‑CACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGA  >  NC_000913/1657916‑1658405

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