Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,024	Δ1 bp	coding (527/696 nt)	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,024	0	C	.	92.6%	210.1 / 9.1	54	coding (527/696 nt)	araD	L‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base . (30/20);  total (32/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 = 9.88e-01

GCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCA  >  NC_000913/65781‑66270                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        gctgtaattCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCaaat                                                                                                                                                                                                                                                 <  1:163546/243‑2 (MQ=255)                               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                                                                                                                                                                                                                                      <  1:176497/115‑4 (MQ=255)                                                                                                                                                    ggTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTgaga                                                                                                                                                                                                                                      >  2:176497/1‑112 (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                                                                                                  >  1:220875/1‑184 (MQ=255)                                                                                                                                                                                                                  ggCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGaa                                                                                                  <  2:220875/184‑1 (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)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        GCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCA  >  NC_000913/65781‑66270

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: ‑