Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,447,045	0	T	.	61.0%	3.9 / 25.1	41	coding (429/1086 nt)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base T (12/4);  new base . (1/24);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 2.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGTCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGA  >  NC_000913/2446745‑2447384                                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                                                                                     cACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGG‑ACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTC‑TTTTTCAGCGCACGCATC‑ACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑AAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAAt                                                                                                                                                                                                                                                                                                                                                           <  1:1283489/308‑1 (MQ=255)                                                               aGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑‑TTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTG‑CGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAAt                                                                                                                                                                                                                                                                                           >  1:1283507/1‑310 (MQ=255)                                                                                                              aCTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑AAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTgccgc                                                                                                                                                                                                                                                                                                                                               <  1:1283511/214‑1 (MQ=255)                                                                                                                         gATGGAGTCGCCCTC‑TTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑‑AAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑aa                                                                                                                                                                                                                         <  1:1283508/303‑1 (MQ=255)                                                                                                                          aTGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑‑AAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑‑AA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGAT‑ATCGCGCACGACGTA‑‑‑‑‑‑‑‑TTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGAc                                                                                                                                                                                                                                   <  1:1283510/295‑1 (MQ=255)                                                                                                                          aTGGAGTCGCCCTC‑TTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑AAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑‑AA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCG‑ATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCgg                                                                                                                                                                                                                                               <  1:1283509/284‑1 (MQ=255)                                                                                                                                                         tCAACTCGTCTAACGCGTCGAGTTTGTACGGGTCCGGGAACAAAAACGGA‑TTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTc                                                                                                                                                                                                                                                                                                                            <  1:430152/188‑1 (MQ=255)                                                                                                                                                           aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑AAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑‑TTCT‑CGGCGAGATA‑tt                                                                                                                                                                                                                                                                                                                                                                     >  1:1283521/1‑146 (MQ=255)                                                                                                                                                            aCTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑‑AAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑‑AA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑aa                                                                                                                                                                                                                         <  1:1283516/267‑1 (MQ=255)                                                                                                                                                            aCTCGTCTAACGCGTCGA‑TTTGTCGGGGTCCGGG‑ACAAAAACGGA‑TTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑‑TTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCgcggg                                                                                                                                                                                                                                                                                                                       <  1:430154/187‑1 (MQ=255)                                                                                                                                                              tCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGC‑‑AAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑aa                                                                                                                                                                                                                         <  1:1283517/266‑1 (MQ=255)                                                                                                                                                                   tAACGCGTCGATTTTGTCGGGGTCCGGGC‑AAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCAccgcc                                                                                                                                                                                                                                                                                                   <  1:1283520/202‑1 (MQ=255)                                                                                                                                                                   aaaCGCGTCGATTTTGTCGGGGTCCGGGC‑‑AAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAAT‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCG‑ATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTATTTTTTTTTTTTGTTTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGggtgg                                                                               <  1:430150/405‑1 (MQ=255)                                                                                                                                                                    aaCGCGTCGATTGT‑ACGGGGTCCGGG‑ACAAAAACGGA‑TTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAGCCGGAATGTCG‑CCCATCTGGGTCAAGGCAGCCACGGATTTCAAATACC‑AAAATTTTCT‑CGGCGAGATA‑‑‑TTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAGAAACGTCCACCGCCGCGa                                                                                                                                                                                                                                                                                               <  1:430151/208‑1 (MQ=255)                                                                                                                                                                            gATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTG‑CGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑AAAACGTCCtt                                                                                                                                                                                                                >  1:1283526/1‑263 (MQ=255)                                                                                                                                                                            gATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTG‑CGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGAT‑ATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑‑AAACGTCCTTAATCGCACTGTAATc                                                                                                                                                                                                 >  1:1283525/1‑276 (MQ=255)                                                                                                                                                                                                  aaaaaCGGATTTTGTCTCGACCTGCGACCAGTCTTTTGATATCCAG‑CGGAATGTCGCCCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATATTTTTTTGGCAAttagc                                                                                                                                                                                                                                                                                                                                                       >  1:1283529/1‑123 (MQ=255)                                                                                                                                                                                                    aaaCGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑‑TTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCaa                                                                                                                                                          <  1:1283524/292‑1 (MQ=255)                                                                                                                                                                                                                       cTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑‑AA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCA‑GGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCg                                                                                                                                                                                                           <  1:1283528/225‑1 (MQ=255)                                                                                                                                                                                                                        tGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCt                                                                                                                                                                                               <  1:1283530/238‑1 (MQ=255)                                                                                                                                                                                                                                    tttGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTgcgcca                                                                                                                                                                                                                                                                                                                                              >  1:1283533/1‑96 (MQ=255)                                                                                                                                                                                                                                      tGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAAT‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGt                                                                                                                                                                                                                                                                               >  1:1283535/1‑160 (MQ=255)                                                                                                                                                                                                                                      tGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGgag                                                                                                                                                                                                                                                                                                                  >  1:1283536/1‑124 (MQ=255)                                                                                                                                                                                                                                         tatCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCGCACTGTaa                                                                                                                                                                                                   <  1:1283531/219‑1 (MQ=255)                                                                                                                                                                                                                                                 cGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑aa                                                                                                                                                                                                                         <  1:1283532/188‑1 (MQ=255)                                                                                                                                                                                                                                                      tGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCg                                                                                                                                                                                                                                                                                >  1:1283542/1‑143 (MQ=255)                                                                                                                                                                                                                                                               aTCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCA‑GGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAA‑CGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTg                                                                                                                                                                                                                                      <  1:1283534/166‑1 (MQ=255)                                                                                                                                                                                                                                                                 cTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCTACGGCGAGATA‑‑TTTTTGGCAATAGCTCCTG‑CGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTATT‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCt                                                                                                                                                      <  1:430159/242‑1 (MQ=255)                                                                                                                                                                                                                                                                  tGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑‑TTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑‑AAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACg                                                                                                                                      >  1:430164/1‑253 (MQ=255)                                                                                                                                                                                                                                                                  tGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTG‑CGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑‑AAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGgag                                                                                                                    >  1:1283543/1‑271 (MQ=255)                                                                                                                                                                                                                                                                    ggTC‑AGGCAGCCACGGATTTC‑AATACC‑AAAT‑TTTCT‑CGGCGAGATA‑‑‑TTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑‑TTTTG‑TTCGt                                                                                                                                                                                                                                                           <  1:430160/141‑1 (MQ=255)                                                                                                                                                                                                                                                                         aGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑‑TTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTAT‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGAA‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTgag                                                                                                                                                                                            <  1:1283537/198‑1 (MQ=255)                                                                                                                                                                                                                                                                            cAGCCACGGATTTC‑AATACCAAAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCA‑GGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTg                                                                                    <  1:1283538/296‑1 (MQ=255)                                                                                                                                                                                                                                                                             aGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGT‑CACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑‑AAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGgagaagaat                                                                                                              >  1:1283544/1‑261 (MQ=255)                                                                                                                                                                                                                                                                             aGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCG‑A‑GGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAA‑CGCGCAGACCGTAT‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTc                                                                                             <  1:1283539/284‑1 (MQ=255)                                                                                                                                                                                                                                                                             aGCCACGGATTTC‑AATACC‑AA‑‑‑‑TTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTgccgccca                                                                                                                                                                                                                                                                                                                                            <  1:430161/61‑3 (MQ=255)                                                                                                                                                                                                                                                                              gCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCg                                                                                                                                                                                                                                                            <  1:1283540/132‑1 (MQ=255)                                                                                                                                                                                                                                                                                cACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑‑TTTTTGGCAATAGCTCCTGCCGCCACGCGCA‑GGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTAT‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGA‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGgtagt                                                                                                                                          <  1:1283541/238‑1 (MQ=255)                                                                                                                                                                                                                                                                                 aCGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGtttcg                                                                                                                                                                         >  1:1283545/1‑205 (MQ=255)                                                                                                                                                                                                                                                                                   ggATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑‑‑TTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGAAAAAAAAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTG‑TTTCGATC‑AAC‑AGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGTCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGa  <  1:430162/374‑1 (MQ=255)                                                                                                                                                                                                                                                                                    gATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATATTTTTTTGGCAAT‑GCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACC‑TA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑‑AAACG‑CCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATCAAACAAGCC‑ATGCTGGTGCCGgta                                                                                                                                            >  1:430165/1‑229 (MQ=255)                                                                                                                                                                                                                                                                                                  taccaattCT‑CGGCGAGATA‑‑‑TTTTGGCAATAGCTCCTGCCGCCACGCGCA‑GGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑‑‑‑TTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGGAAAAA‑AAAACGTCCTTAATCGCACTGTAATCCTg                                                                                                                                                                                              <  1:430163/166‑1 (MQ=255)                                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                                                                                     CACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTG‑CTCGACCTGCGACCAGTC‑TTTGATATCCAG‑CGGAATGTCG‑CCCATCTGGGTC‑AGGCAGCCACGGATTTC‑AATACC‑AAA‑‑TTTCT‑CGGCGAGATA‑TTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGG‑TTTCGCGGGCGGAAG‑AACGTCCACCGCCGCGATAATCGCGCAGACCGTA‑‑‑‑‑‑‑TTTTTG‑TTCGTAGGTGTAATCGGCATGGCCTGGACGG‑‑‑‑‑‑AAAACGTCCTTAATCGCACTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATC‑AACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGTCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGA  >  NC_000913/2446745‑2447384

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