Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867	C→T	S177S (AGC→AGT)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867	0	C	T	97.3%	101.2 / NA	37	S177S (AGC→AGT)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/3);  minor base G (0/1);  total (33/4)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTA  >  NC_000913/627‑1112                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      aTTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATACGGTCGAAAAACTGCTGGCAGAGTGGCATTACCGCGAATCTACCGTCGATATTGCTGAGTCCCCCCGCCGTATTGAGGCCAGTGGCATTCCgg                                                                                                                                                                                                                                             >  2:325789/1‑251 (MQ=255)     gCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCCGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTTCGGCAAGTCGCATTCCGGCtt                                                                                                                                                                                                                                          >  2:78575/1‑249 (MQ=255)      cccAAATAAACCATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTCTTAGAAGCGCGCGGTCAAAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCCGTGGGGCATTACCTCGAATCTACCGTCGGTATTGCTGAGTCCACCCGCCGTATTTGGGCAAGTCGCATTCCGGCTGAt                                                                                                                                                                                                                                        >  2:83639/1‑251 (MQ=255)                      ccTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGCATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATgg                                                                                                                                                                                                                        >  2:214416/1‑251 (MQ=255)                       cTGCATGGCATTAGTTTGTTGGGGCCGTGCCCGGATAGCATCACCGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCAATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTTCCGTTAACGATCCGGTCGAAACACGGCTGGCAGTGGGGCCTTTCCACGAATCTACCGTCGATTATGCTGAATCCCCCCGCCGAATTGCGGCAAGTCGCATTCCGACTGGTCACATGGGACTGATgg                                                                                                                                                                                                                        >  2:398358/1‑250 (MQ=255)                                             ggCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGa                                                                                                                                                                                                  >  1:549299/1‑250 (MQ=255)                                                     ccGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCt                                                                                                                                                                                                                                           >  1:513121/1‑201 (MQ=255)                                                     ccGGATAGCATAAACGCTGCGCTGAATTGCCGTGGCGAGAAAATGTCGAGCGCAATACTGGCCGGCGTATTAGAAGAGCGCGGTCTCAACGTTCCCGTTATCGATCCGGTCGAAAGACTGCTGGCGGTGGGGCATTGCCGCGACTCTACCGTCGATATTGCTGAGTCCACCCGACGTATTTCGGCAAGTCGCTTTCCGGCTGTTCACATGGGGATGGTGGCCGGTTTAACCGCCGGTAATGAAAAAGGCgg                                                                                                                                                                                         >  2:233340/1‑250 (MQ=255)                                                      cGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGa                                                                                                                                                                                         >  1:215441/1‑250 (MQ=255)                                                               tCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAACGGCGGACTggtggt                                                                                                                                                                                >  1:676358/1‑250 (MQ=255)                                                                   cgctgcgcggATTTGCCGTGGCGAGAAAATGTCGCTCGACATTATGGCCGGCGGATTAGAAGCGCGCGGTCTCGACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCAGTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCACGTCGCATTCAGGCTTGTCACATGGAGCTGAAGGCAGGTTTACCCCCCGGTAGTGAAAAAgaagaaccggggggaattg                                                                                                                                                                           >  2:494519/1‑233 (MQ=255)                                                                        cgctgATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGATATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCCTTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGGGATGGTGGCAGGTTTCACCCCCGGTAATGAAAAAGGCGAACTGGGGGAACTTGGACGc                                                                                                                                                                      >  2:128778/1‑251 (MQ=255)                                                                          ctgATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGGGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCa                                                                                                                                                                     >  1:411671/1‑250 (MQ=255)                                                                            gATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGGACTTGGACGCaccg                                                                                                                                                                  >  1:524860/1‑248 (MQ=255)                                                                                        cGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCCACGGTTCCGacacct                                                                                                                                                      >  2:669351/1‑247 (MQ=255)                                                                                               aTGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACAGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGCCGCAACGGTTCCGACTActccgcggc                                                                                                                                               >  1:43059/1‑245 (MQ=255)                                                                                                             tATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGGGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGGAATGAAAAAGGCGAACTGGCGGGACTTGGGCGCAACGGTTCCGACTActccgaggcggggggtggtgacc                                                                                                                                 >  2:581259/1‑231 (MQ=255)                                                                                                                 gCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGAGGTGGCTGCCTGttt                                                                                                                             >  2:441901/1‑251 (MQ=255)                                                                                                                    ggCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTa                                                                                                                            >  1:428024/1‑249 (MQ=255)                                                                                                                      cGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACtg                                                                                                                                                                                     >  1:644604/1‑190 (MQ=255)                                                                                                                         aTTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGGCGCCTGTTTTCGCGccc                                                                                                                     >  2:94603/1‑250 (MQ=255)                                                                                                                               aGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCGTTTTACGCCCCGgttgtt                                                                                                               >  2:88358/1‑251 (MQ=255)                                                                                                                                                     gTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGCTATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTTCGAGATTTGGGAGGACGTTgg                                                                                         >  1:583021/1‑250 (MQ=255)                                                                                                                                                        aTCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGCTTGTTGAGAGATTTGGACGGACGTTGACgg                                                                                      >  1:55842/1‑251 (MQ=255)                                                                                                                                                          cGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGAGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAACAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGGTGGCTGCCTGTTTACGCGCCGATTGTTGCGGGATTTGGAGCGACGTTGACgggg                                                                                    >  1:669483/1‑251 (MQ=255)                                                                                                                                                                            gctgGCAGTGGGGCATTACCTCGAATATACCGTCGATATTGCTGAGTCCACCCGCAGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGGATTGGACGGACGGTGACGGGGGCGATACCTGCGACCcgc                                                                  >  1:250280/1‑251 (MQ=255)                                                                                                                                                                              tgGCAGTGGGGAATTACCTCGAATCTACCGTCGATATTGCTGAGACCACCCGCCGTAATGAGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACAACTCAGCGACGGGGGCGGCGGCCTGTTTACGAGCCGATTGGTGAGAGGATTGGACGGGCGTTGAGgggggggatacaggcgagccgcgg                                                                >  2:642952/1‑231 (MQ=255)                                                                                                                                                                               ggCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGACGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGAAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGGCTGGCTGCCTGTTTACGGGACGATTTTTGCGAGATTTGGGCGGACGTTGACGGGGTCTTTACCTGCGACCCGCGTc                                                               >  1:2111/1‑251 (MQ=255)                                                                                                                                                                                                      cTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCCGGTTTCACCGCCGGTAATGAAAAAGGCGACCTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCCGTTTCCGCGCCGGTTGTTGCGAGTGTTGGGCGGAAGTTTAAGGGGGCTATGCCTGCGACCCGCGTCAGGTTACCGATGCGAGGgtgttg                                        >  2:368464/1‑251 (MQ=255)                                                                                                                                                                                                      cTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGGTTTGGACGGACGTTGACGGGGTCTAAACCTGCGACCCGCGGCAGGGGCCCGATGCGAGGttgttg                                        >  1:450362/1‑251 (MQ=255)                                                                                                                                                                                                       ccggtgtttcttggttgcccccccccGCCTTATGGCGGCAAGTGGCATTCCGGTTTATCACAGGGTGCTGAGGGCGGGTTTCCCCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGa                                       <  1:581259/228‑1 (MQ=255)                                                                                                                                                                                                         cggttaattggcatctcccccccggccgTTTTGGGGAAAGGCGCATTCCGGTTGATCCAATGGGGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAg                                     <  1:142722/227‑1 (MQ=255)                                                                                                                                                                                                                       gAGTCCACCCGCCGTATTGAGGCAAGTCGCATTCCGGATGATCACCTGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGGCTGGCTGCCTGTTTACGGCCCGGATGTTGCGCGGTTTGGACGGACGTTGGCGGGGTTTATACCTTCAACCCCCGTCAGGGGCCCGGTGGGAGGGTGTTGAAggcggtggctccacc                       >  2:290775/1‑237 (MQ=255)                                                                                                                                                                                                                            cACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGGGAGGTTTGGACGGAGGTTGCCGGGGTCTATACCTGCGCCCCGCGGCAGGTGCCCGCAGCGAGGGTGTTGAAGTCGATGTCCtccccgga                    >  2:292341/1‑242 (MQ=255)                                                                                                                                                                                                                                  ccgTTTTGCGGCAAGTCGCATTCCGGCGGATCACATGTGGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATgg            <  1:590559/251‑1 (MQ=255)                                                                                                                                                                                                                                            gCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTAACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGGGGTTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGGTGCGGGGTTGTTGAAGTCGATGTCCTACCAggcagtagagtgggagggcgt  >  2:596651/1‑232 (MQ=255)                                                                                                                                                                                                                                            gCAAGTCGCATCCGGGCTGATCACAGGGTGCGGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGc                 <  1:596651/236‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      ATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTA  >  NC_000913/627‑1112

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