Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,032	C→T	P173P (CCG→CCA)	araD ←	L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,032	0	C	T	96.4%	164.2 / ‑1.6	55	P173P (CCG→CCA)	araD	L‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (32/21);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGC  >  NC_000913/65811‑66276                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      ttataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCCAATGGGCCGTGCGAATGGACCAGGACGCCGGGc                                                                                                                                                                                                                         >  1:426025/2‑251 (MQ=255)                gTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGAGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACAGCCGGCAGATTATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCg                                                                                                                                                                                                          >  1:341356/1‑251 (MQ=255)                        ggCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTg                                                                                                                                                                                                  >  2:356135/1‑251 (MQ=255)                         gCTGTGGTTTTATACAGTCATTCCTGCACGTAATCTGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGt                                                                                                                                                                                                 >  2:64720/1‑251 (MQ=255)                         gCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATCGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGt                                                                                                                                                                                                 >  2:206653/1‑251 (MQ=255)                         gCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCCTCGATACCCTg                                                                                                                                                                                                  >  2:139186/1‑251 (MQ=255)                                    ataCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGt                                                                                                                                                                                      <  1:356135/251‑1 (MQ=255)                                        aGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCt                                                                                                                                                                                  <  1:12742/251‑1 (MQ=255)                                                  gCCCGTAATATGCCTTCGCGCCATGCTTACGCAGCTAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCCCGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGt                                                                                                                                                                        >  2:438626/1‑251 (MQ=255)                                                      gTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTAcc                                                                                                                                                                    >  1:310626/1‑251 (MQ=255)                                                         atatGCCTTCGCGCCATGCTTACGCAGATAGAGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATAGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGt                                                                                                                                                                 >  1:412876/1‑251 (MQ=255)                                                                     gcaaTGCTTACGCAGATAGTGTTGATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGAGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCata                                                                                                                                                     <  1:150033/248‑1 (MQ=255)                                                                     gcCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCata                                                                                                                                                     <  1:64720/251‑1 (MQ=255)                                                                         tGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTc                                                                                                                                                  >  1:199514/1‑250 (MQ=255)                                                                          gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTca                                                                                                                                                 <  1:5153/250‑2 (MQ=255)                                                                          gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGAATATCCGTAAACTGCGGCGCTAACTGACGGCAGAATACCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTGCGGCATTTTTGCCC‑ATGCAAATGGGCCGTTAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTca                                                                                                                                                 <  2:358278/250‑2 (MQ=255)                                                                           cTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTcacc                                                                                                                                               >  1:235404/1‑248 (MQ=255)                                                                            ttACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCg                                                                                                                                              >  2:38959/1‑251 (MQ=255)                                                                              aCGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCAGAACCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGAGACCTCTTCCAGCACGATGGCGTTATGCACCGAATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTAGAAGGTTTCTACGATGACGTTACCGGTTTCCCGCTCATATTcacagtt                                                                                                                                            >  2:41979/1‑245 (MQ=255)                                                                                  aGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAtt                                                                                                                                        <  1:428399/251‑1 (MQ=255)                                                                                       gtgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCCTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGc                                                                                                                                   >  2:124492/1‑251 (MQ=255)                                                                                        tgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGTTACCCTGTTTTTCGCAGGTTTCTACGATTGCGTTACCGGTTTCCCACTCCTATTCACCTTTGATTTCTGcc                                                                                                                                  >  2:412342/1‑250 (MQ=255)                                                                                        tgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCAGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCg                                                                                                                                  >  1:129173/1‑251 (MQ=255)                                                                                             attcagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGTCAGAATTTCACCATATAAGCGACCTCTTCCACCACGATGGCGTTATGCACCGCATCTTCGACCTTTTTGCCC‑ATGCAAATGGGCCGCGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGCCGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGt                                                                                                                             >  1:263057/4‑251 (MQ=255)                                                                                               tcagcagCGTTTGCTGAATATCCGGTAACTGCGGCGCTAACTGACGGCAGACTATCCACATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACAGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTc                                                                                                                            >  1:297911/2‑250 (MQ=255)                                                                                                           gctgcATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTg                                                                                                                <  1:38959/250‑1 (MQ=255)                                                                                                                   tCCGGTAACTGCGGCGCTAACTGACGGCAGCATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCG‑GGGaa                                                                                                        >  2:325874/1‑249 (MQ=255)                                                                                                                      ggTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTgaga                                                                                                                                                                                                                                            >  2:176497/1‑112 (MQ=255)                                                                                                                      ggTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTgaga                                                                                                                                                                                                                                            <  1:176497/115‑4 (MQ=255)                                                                                                                                  gcTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTCGAAAt                                                                                         >  1:251870/1‑249 (MQ=255)                                                                                                                                     aactgcgagccaacgatCCCCATATAAGCGACCTCTTCCAGCACGACGGAGTTATGCACCGCATCCTCGGCATTGTTGCCC‑ATGCAAATGGGCCGTGAGAAGGGTCCAGGACGCCTGGCAATTGCGCTGCCTCGATACCCTGTTATTCGAAGGTATCTACGATGACGCGATCGGGTTCCCGCGCATATTCACCGTTGATTTCTACGTCTGTCATTTTGCGG‑AGCAAGGGAATGGTACAGTAGAAATGGTCg                                                                                    >  2:106129/16‑251 (MQ=255)                                                                                                                                     aaCTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGt                                                                                      <  1:47674/249‑1 (MQ=255)                                                                                                                                     aaCTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTc                                                                                     <  1:445002/250‑1 (MQ=255)                                                                                                                                         gACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGc                                                                                  <  1:139186/250‑1 (MQ=255)                                                                                                                                             gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAACTAGTCGGCgtgggg                                                                            >  2:194296/1‑250 (MQ=255)                                                                                                                                               aGAAAATCACCATATAAGCGACCTCTGCCAGAACGATGGCGTGATGCACCGCAGCTTCTGCATTTTTGCCC‑ATGCAAATGTGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCTCTCATATTCACCGTTGATTGCTGCGTCGGTCATTTTGCGGGGGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCttgggtgt                                                                          <  2:93031/251‑2 (MQ=255)                                                                                                                                                 aaTATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGgtggtg                                                                        >  1:226186/1‑251 (MQ=255)                                                                                                                                                  atatCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGc                                                                       >  1:228097/1‑251 (MQ=255)                                                                                                                                                               aGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCGTCGGCATTTTTGCCC‑ATGCATATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTGTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGGCGGCGTGGGTGGTGCCGGTTGCTGGAAt                                                          <  2:328883/251‑1 (MQ=255)                                                                                                                                                                 cgtccTCTTCCAGCACGGTGGCGTTTTGCACCGCACCTTCGGCATTTTTGCCC‑AGGCAAATGGGCCGTGAGAATGGACCAGTACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGAGGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTTCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCCGCGTGGGTGGTGCCGTTTGCTGGAATCg                                                        <  2:341356/248‑1 (MQ=255)                                                                                                                                                                    ccTCTTCCAGCACGATGGCGTTATTCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                     <  2:415002/251‑1 (MQ=255)                                                                                                                                                                           cAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGc                                              >  1:170752/1‑251 (MQ=255)                                                                                                                                                                                   tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCa                                                                                                                                                       <  1:41307/139‑1 (MQ=255)                                                                                                                                                                                   tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCa                                                                                                                                                       >  2:41307/1‑139 (MQ=255)                                                                                                                                                                                   tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccc                                      >  1:157992/1‑251 (MQ=255)                                                                                                                                                                                    ggCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGaa                                                                                                        <  2:220875/184‑1 (MQ=255)                                                                                                                                                                                    ggCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGaa                                                                                                        >  1:220875/1‑184 (MQ=255)                                                                                                                                                                                            gCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCGAGAtagtgg                             >  2:413474/1‑251 (MQ=255)                                                                                                                                                                                               ccGCATCTTCGGCATTTTTGCCCC‑TGCAAATGGGCCGTGCGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGCAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTTCCGTAGGAATAGTCTGCGGGGGTGGAGGCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAgtgagtggcgg                          >  2:54415/1‑241 (MQ=255)                                                                                                                                                                                                     ccttcgtcattttGGCCCATGCAAATGGGCCGTGAGAATTGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTTCCGTTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGgcgc                     <  1:272575/242‑1 (MQ=255)                                                                                                                                                                                                       tCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGTCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGTAATCGACTGGCCCGCCTGCGCCCAGATAGtggcgtggcg                      <  1:364440/247‑1 (MQ=255)                                                                                                                                                                                                        cGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATCGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAg                  <  1:206653/250‑1 (MQ=255)                                                                                                                                                                                                            aTTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCCGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtggagggagggggt             >  1:443435/1‑240 (MQ=255)                                                                                                                                                                                                                   gcccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCTTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAAGGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTGGCTGGAATCGTCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGca        <  2:334298/249‑1 (MQ=255)                                                                                                                                                                                                                        tGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCAGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTgtgggtgcgcaatgc  >  2:155069/1‑237 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      GTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGC  >  NC_000913/65811‑66276

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