Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,492,537	C→T	intergenic (‑12/+50)	idnT ← / ← idnO	L‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,492,537	0	C	T	97.7%	141.4 / ‑3.0	43	intergenic (‑12/+50)	idnT/idnO	L‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/21);  minor base . (1/0);  total (22/21)
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.82e-01

CGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGA  >  NC_000913/4492292‑4492780                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATACGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTtta‑tt                                                                                                                                                                                                                                                >  2:278824/1‑251 (MQ=255)         aaaTCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCTGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAg                                                                                                                                                                                                                                        <  2:64671/251‑1 (MQ=255)          aaTCAGCTTGCCAAACATGGCCCCAAAACCGAGGAATATGGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGa                                                                                                                                                                                                                                       <  2:533986/251‑1 (MQ=255)              aGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGaaa                                                                                                                                                                                                                                     >  2:285485/1‑249 (MQ=255)                 ttGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCCGCTACCAGACCGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGAGCGACGCCTGCCGCAATAATGATTAATGGCATAGACTTTATCTTTTA‑TTAGCGCAAGaaacgc                                                                                                                                                                                                                                  >  2:300778/1‑246 (MQ=255)                  agccgatcATGGCCCCGAACCCGAGGTTATTTCCCACCCCGACGGGCGTGCTGCCGATACCGTTTTGTATAGATTGGAGGAGGGCCTGGGCATCCATCCCTTCGGCAAATCCGTCGACGGCAGCTACCAGAACGAGGACAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCTCGGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAagcag                                                                                                                                                                                                                                <  2:287147/243‑1 (MQ=255)                               cccGAAACCGAGAATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGc                                                                                                                                                                                                                  <  2:218730/251‑1 (MQ=255)                                     cccGAGGGTCATGGCCTCCCCCCCGGGCGTTCTCCCGATACCATTTTGTATAGAGTGCAGGCCGGCCTGTACATCCATCCCTTCTGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCTCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCATCGCTTTCTCGCTGCACGTAg                                                                                                                                                                                                            <  2:250979/250‑1 (MQ=255)                                        gAGGTTCATTACCAGCCCGCCGAGCGTGCTGCCGATACCACTTTGTATAGAGAGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATACGATGACGGCAGCTACCAGAACGCGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCACTAATGATTAATGGCATAACCTTTATCT‑TTACTTAATGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgagg                                                                                                                                                                                                         >  2:582284/1‑248 (MQ=255)                                         aGGATCATTGCCGGCGCGCCGAGCGTGCGGCCTATACCATTTTGTATAGAGTTCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAggggg                                                                                                                                                                                                        <  2:423323/251‑1 (MQ=255)                                                    cAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAAc                                                                                                                                                                                              >  2:58730/1‑250 (MQ=255)                                                             gAGCGGGCTGCCGATACAATTTTGTATAGAGGGCAGGACGGCCTGTGCATCCACCCCTTCGGAAAATCCGACGGCGGCGGCTCCCAGAACGAGGGAAATAAAGCCGTTAATTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGGTAAACAGCCACTAAc                                                                                                                                                                                    <  2:383817/251‑1 (MQ=255)                                                                        cGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATc                                                                                                                                                                         >  1:134840/1‑251 (MQ=255)                                                                               aTTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaac                                                                                                                                                                  >  2:513750/1‑251 (MQ=255)                                                                                ttttGTATAGAGTGCAGGACGGCCTGGGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTGAAACAGCCACTAACATGCCGCCATCaacaaac                                                                                                                                                                  >  1:232204/1‑250 (MQ=255)                                                                                   tGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAaacaac                                                                                                                                                               >  2:282409/1‑250 (MQ=255)                                                                                     tataGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTtta‑tta                                                                                                                                                                                                                                               <  2:233904/168‑1 (MQ=255)                                                                                     tataGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTtta‑tta                                                                                                                                                                                                                                               >  1:233904/1‑168 (MQ=255)                                                                                      ataGCGTGCAGGTCGGCCTTTGCATCCTGCCTTTCGGCAAACCCGACGACGGCATCTACCAGAACGAGGGCAAGAAAGCCGTCAACTTTAAAGCCGATCATCAGAATAAGATGCAGCGCGACGCCTGCCGAAATAATGATTAATGGCTTAGTCTTTTTTTTTTA‑TGAGTGCCAGAAAGGCAGCGTTTTCGCGCTGTACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGt                                                                                                                                                           <  1:311598/251‑1 (MQ=255)                                                                                           gTGCAGGACGGCCTGTGCATCCACCCCTTGGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTAAACTTAAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAGCAACAGGTGGCCg                                                                                                                                                      <  2:359920/251‑1 (MQ=255)                                                                                                ggACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCAACGCCTGCCGCAATAATGATTAAGGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTAc                                                                                                                                                 <  1:30984/251‑1 (MQ=255)                                                                                                    ggCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTCCCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGaa                                                                                                                                              >  2:33078/1‑250 (MQ=255)                                                                                                    ggCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGaa                                                                                                                                              >  1:20015/1‑250 (MQ=255)                                                                                                         gtgCATCCATCCCTTCGGCAAATCCGTCGACGGCAGCTACCAGACCGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTTGCCGTTTACGAAATCag                                                                                                                                         <  1:256228/250‑1 (MQ=255)                                                                                                               ccatccCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCtt                                                                                                                                   <  1:495738/250‑1 (MQ=255)                                                                                                                  acccTTCGACAAATCTGACGAGGGAAGCTACCAGAAGGAGGGAAATAAAGCTGTTAACTTTAAAGACGATCATCAGAATAAGAAGTAGCGCGACGCTTGCCGCAATAATGATTAATGGCATAGTTTTTATTTTTTA‑TTAGTGCCAGAAAAGCAGTGTTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTga                                                                                                                               <  1:278824/250‑1 (MQ=255)                                                                                                                                   ccgacgGCAGCTACCAGACCGAGGGCAATAAAGCCGTTAACTTTAAAGCCAATCATCAGGATAGTAAGCAGCGCAACGCCTGCCGCAATAATGATTAATGGCATATTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGTACACCgcag                                                                                                              <  2:296438/250‑1 (MQ=255)                                                                                                                                          cAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAAt                                                                                                       >  2:246130/1‑251 (MQ=255)                                                                                                                                                                     cGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGCCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTa                                                                                                                                                  <  2:62263/181‑1 (MQ=255)                                                                                                                                                                     cGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTa                                                                                                                                                  >  1:62263/1‑181 (MQ=255)                                                                                                                                                                      gTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTgcg                                                                           >  1:154768/1‑251 (MQ=255)                                                                                                                                                                                  gacGATCAACCGGATAAGAAGCAGTGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGAGCCAGAAAAGCAGAGCATTCGCGCTGCAAGCAGCGGGGCTGATAAACAGCCCCTCACATGCCGAAATCATCAAACAACAGGTGGCCGTTTACGAAATAAGAGGCTTTTGAAGAAAGGAACCCCGCAGCACCAATCAGTTCCTGCGGATCTCCACAGCGTGCGGCGGCGGtgcgg                                                               >  2:295334/3‑247 (MQ=255)                                                                                                                                                                                     gATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCaga                                                                                                                                        <  1:143904/175‑1 (MQ=255)                                                                                                                                                                                     gATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCaga                                                                                                                                        >  2:143904/1‑175 (MQ=255)                                                                                                                                                                                           caGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAAc                                                      >  2:122784/1‑251 (MQ=255)                                                                                                                                                                                              gATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGGCCGTTTGCACAAcc                                                     >  2:151197/1‑249 (MQ=255)                                                                                                                                                                                               aTAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAgg                                                  <  1:151197/251‑1 (MQ=255)                                                                                                                                                                                                              aCGCCTGCAGCAATAATGATTAATGGCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTTAAGGGTTCGt                                   >  1:49263/1‑251 (MQ=255)                                                                                                                                                                                                                              ggATTAATGGCATAGGCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCttt                   <  1:94161/250‑1 (MQ=255)                                                                                                                                                                                                                                       gCATAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTc          <  1:636042/251‑1 (MQ=255)                                                                                                                                                                                                                                        cATAGTCTTTATCTTTTA‑TTAGGGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCt         <  1:33078/251‑1 (MQ=255)                                                                                                                                                                                                                                         aTAGTCTTTATCTTTTA‑TTAGTGCCAGAAAAGCAGCGATTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGGCATTTCTg        >  1:224743/1‑251 (MQ=255)                                                                                                                                                                                                                                                ttATCTTTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGATTCGGCCTCAAACCGGGGTTTAGTCATTTATGTGTTGa  >  1:631763/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTA‑TTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGA  >  NC_000913/4492292‑4492780

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