Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,415	C→T	A31A (GCC→GCT)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,415	0	C	T	97.4%	108.4 / ‑6.5	38	A31A (GCC→GCT)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/8);  total (29/9)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGA‑TCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTG  >  NC_000913/4617189‑4617652                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              cGCATTTCAGAGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCaaatccccg                                                                                                                                                                                                                       <  2:586572/251‑7 (MQ=255)       aCAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGc                                                                                                                                                                                                                 >  1:548546/1‑251 (MQ=255)                       gTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCt                                                                                                                                                                                                 >  2:36444/1‑251 (MQ=255)                                   ccTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGCCGAGAAAGTAA‑TTGCTCTGGGTCATCAGGCCAAAACCCCGGTCGCCAATACCGCAGCTATATGTATCTATCCTc                                                                                                                                                                                     >  2:604532/1‑251 (MQ=255)                                       aaTTGTTATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGAAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCGGCGTGCACTGAAATTGATGGACCTGACCCCCCTGAATGACGACGACACCGACGAGAAATTAA‑TTGCCCTGTGTCATCAGGCCAAAACCCCGGTCGTCAATACCGCCGCTATCTGTATCTATCCTCGCt                                                                                                                                                                                  <  1:330293/250‑1 (MQ=255)                                           gtgATGTAAATCGAAGTGTGTGGCGGAGTAGATGTTAGAATACTAACAACCTCGAAAGGTAAATTTTATTGGCGACAAGCCAGGAGAATGAAAGGACTGTTCTAAAAGCAAGCAGCCGGCTTGCACTGAAATGGAGGGCCCTGACCCCCCGGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCGGGCAAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATc                                                                                                                                                                             <  1:621642/251‑1 (MQ=255)                                            tgATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGCAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATcc                                                                                                                                                                            >  1:53973/1‑251 (MQ=255)                                            tgATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATcc                                                                                                                                                                            >  1:504358/1‑251 (MQ=255)                                                  atATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTg                                                                                                                                                                      >  1:367427/2‑251 (MQ=255)                                                  atATCAAAGTGTGGTGGGGATGAGATGTTAGAATCCTAACAAACGCGCAAGGGTATTTTTATTGGCGACAAGCCAGGAGAATGAAATGATTACTCTGAAACCAAGCACCCTGCGTGCACTGAAATTGTTGGACCTGACCACCCTGAATGACAACGACCCCGACAAGAAAGTAT‑TTGCTCGGTTTCATCAGGCCAAACCCCCGGTCGGCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAtt                                                                                                                                                                       <  1:378102/249‑1 (MQ=255)                                                    aTCGAAATGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAATCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGCAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCt                                                                                                                                                                    >  2:211407/1‑251 (MQ=255)                                                       gAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGc                                                                                                                                                                 >  1:146221/1‑251 (MQ=255)                                                            gtgtTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGCCGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTCACTATcc                                                                                                                                                                                       >  2:215847/1‑224 (MQ=255)                                                            gtgtTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAAc                                                                                                                                                            >  1:284801/1‑251 (MQ=255)                                                            gtgtTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATGGGCGAAAAGCCGGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGAGGGACCTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTAT‑TTGCTCTGTGTCATCAGGCCAAACCCCCGGGCGGCAATACCGCCGCTATCTGTATCTATcc                                                                                                                                                                                       <  1:215847/224‑1 (MQ=255)                                                            gtgtTGCGGAGAAGATGTTAGAATACTAACAAACTCGCAAGGGGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACAGAAATTGATGGACCTGAACACCCTGAATGACGACGACACCGACGAGGAAGTAA‑TTGCTCTGTGTCATCCGGCCAACACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTACCCCGATTGCTCGCaaaa                                                                                                                                                             >  1:333255/1‑250 (MQ=255)                                                                tggggagaggATGTTAAAATACTAAAAACCCCGAAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAAGGACTGATCTAAAAGCAAGAAGCCTGCGTGCACTAAAATTGATGGACCTGACCCCCCTAAATGACGACGACACCGAGGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGCCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAAcact                                                                                                                                                         <  1:602927/241‑4 (MQ=255)                                                                                 aatcaaaacaaccgctaaaGGTAAATTTTATTGGCGAAACGCAAGGAAAATGAAAGGACTGATTTGAAAGTAAGCAGCCTGCTTACATTGAAATTGAGGGCCGTGACCACCCTGAATGACACCGACACCGACGAAAAATTAA‑TTGCTTTGTGTCACCATCCCAAACCCCCGGTCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAcccc                                                                                                                                       <  1:386647/234‑1 (MQ=255)                                                                                          acaCTCGCACGGTGACTTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGCTCTGAAAGAAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTAGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCg                                                                                                                              >  2:635825/3‑251 (MQ=255)                                                                                                                           cAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTa                                                                                              >  2:281345/1‑250 (MQ=255)                                                                                                                                   tgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                                      >  1:341916/1‑250 (MQ=255)                                                                                                                                      aatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCACCTTCCCACACGGTAACGAtgacatcg                                                                                  >  2:255707/1‑251 (MQ=255)                                                                                                                                       atgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGGATTGCTACGGAAACCAACTTCCCACACGGTAACGAtgccatcg                                                                                  >  2:395715/1‑250 (MQ=255)                                                                                                                                                  gAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTgg                                                                      >  1:144480/1‑251 (MQ=255)                                                                                                                                                     agcaagcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTgcc                                                                     >  1:400883/1‑247 (MQ=255)                                                                                                                                                      gcaagcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGc                                                                     >  2:142028/1‑248 (MQ=255)                                                                                                                                                           cagcCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCtgc                                                                                                                     >  1:164162/1‑193 (MQ=255)                                                                                                                                                            agcCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTACCCCGATTGCTCGCAAAACACTGAAAGAGCAGTGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCg                                                            >  1:43902/1‑251 (MQ=255)                                                                                                                                                                       aCTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGAAGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATcc                                                 >  1:11926/1‑250 (MQ=255)                                                                                                                                                                            aaTTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAC‑TTGCTCTGAGTCATCAGGCCAAAACCCCGGTAGGCAATACCGCCGCTATCTGTATCTATCATCGCTTTATCCCGATTGCTAGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTCCGCTAACAAACTTCCAACACGGTAACGATGACATCGAAATAACGCTGGCAGAAACCCGTGCGGCCATCGcatcc                                            >  2:174079/1‑247 (MQ=255)                                                                                                                                                                                gATGGACCCGCCCACCCTAAATGACGACGACACCGACGAAAAAGTAA‑TTGCTTTGTGTCACCAGACCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGt                                         <  1:299227/250‑1 (MQ=255)                                                                                                                                                                                   ggACCTGACCACCCTGAATGACGACGACCCCGACGAGAAAGCAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATATGTATCTATCCTCGCTTTATCCCGATTGCACGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCAGCATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAACTCGAGCTGCCAGGAACACTTGCGGCAATCGCCTACGGTGccg                                     >  1:672464/1‑249 (MQ=255)                                                                                                                                                                                       cTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatg                                  >  1:593478/1‑250 (MQ=255)                                                                                                                                                                                        catctccccTGAATGACGACGACACCGACGA‑AAAATAAGTTGCTCTTGGTCATCAGCCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGaa                                <  1:39305/245‑1 (MQ=255)                                                                                                                                                                                             accCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGaa                                >  2:39305/1‑246 (MQ=255)                                                                                                                                                                                                   aaTGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGtggg                     >  1:219543/1‑248 (MQ=255)                                                                                                                                                                                                                     tagaAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGTAAATCCGTATTGCTACGGTAACCAACTTCCAACACGGTAACGATGACCTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGAAGTGGAGTTCCCGTAACGCGCGCt   >  2:509213/2‑251 (MQ=255)                                                                                                                                                                                                                       gaAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTg  >  1:244064/1‑250 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              CGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGA‑TCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTG  >  NC_000913/4617189‑4617652

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