Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	774	T→C	T146T (ACT→ACC)	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	774	0	T	C	97.6%	126.1 / ‑4.3	41	T146T (ACT→ACC)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (33/7);  total (34/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

ATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGA  >  NC_000913/533‑1019                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       aTGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAAAGTTCCCGTTATCggt                                                                                                                                                                                                                                               >  2:430324/1‑249 (MQ=255)   gtcgatcGTATTTTGTCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAACCTTGCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCGGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGAATCAATGCTGCGCTGATTTGCCGTGGCGAGAAAATTTCGTTCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATcc                                                                                                                                                                                                                                             <  2:52090/245‑1 (MQ=255)                  gCCGAACGTTTGACGGGACTCGCCGCCGCCCACCCGGGATTCCCGCTGGCGCAATTTAAAACTGTCGTCGATCAGGAATCTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGACAAActgct                                                                                                                                                                                                                              <  2:29612/251‑1 (MQ=255)                   ccGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAActgctg                                                                                                                                                                                                                             >  1:22972/1‑251 (MQ=255)                          tttGACGGGACTCGCCCCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTACTGCATGGCATTAGTTTGTTGGGGCAGTGACCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGCTCCGGTCGAAAAACTGCTGGCAGT‑gg                                                                                                                                                                                                                     >  1:228733/1‑251 (MQ=255)                               cGGGACTAGCCACCGCCCAGCAGGGATTCCCGCTGGCGCACTTGAAAACTTTCGTCGCTCAGGACATTGCCCAAATAAAACATGTCCTGCATGGCATTAGCTTGTTGGGGCAGTGCCCGGATAGCTTCACCGCTGCGCTGATTTGCCGTGGCGAGCAAATGTCGCTCGCCATTATGGCCGGCGCATTAGAAGCGCGCGGTCCCAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGccagggggcat                                                                                                                                                                                                                 >  2:244473/1‑241 (MQ=255)                                              cccAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTAcc                                                                                                                                                                                                 <  1:215490/251‑1 (MQ=255)                                               ccAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCg                                                                                                                                                                                                >  1:173648/1‑251 (MQ=255)                                                   ccGGGCTTCCAGCTGGCGACATTGAAAACTTTCGTCGATCAGGAAGTTGCCCAATTAAAACATGTCATGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAAAGCTGAGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGAGAGGTCAAAACATTACCGTTATTGATCCGGACGAAAAAATGCTGGCAGTTGGGACAT‑ACCATGAATCTACCGTCgg                                                                                                                                                                                            >  2:305385/1‑250 (MQ=255)                                                      ggattCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGAATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCTGTCCGGTCGAAAAACTGCTGGCAGT‑GGGGCCTTACCTCGAATCTACCGGCGatat                                                                                                                                                                                         >  2:445854/4‑251 (MQ=255)                                                       gattCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAACCTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATAtt                                                                                                                                                                                        >  1:358628/3‑251 (MQ=255)                                                        attcccgcgggtcgaaTTGAAAACTTTCGTCGATCAGGAATTGGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCATCGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTg                                                                                                                                                                                       <  2:103520/237‑1 (MQ=255)                                                           cccGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCt                                                                                                                                                                                     <  1:68595/250‑1 (MQ=255)                                                                tGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTcc                                                                                                                                                                               >  1:104052/1‑251 (MQ=255)                                                                   cgcAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCAcc                                                                                                                                                                            >  1:67361/1‑251 (MQ=255)                                                                           aaaaCTTTCGTCGATCAGGAATTTGCCCAAATAAAACACGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCg                                                                                                                                                                                             >  1:111062/1‑226 (MQ=255)                                                                           aaaaCTTTCGTCGATCAGGAATTTGCCCAAATAAAACACGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCg                                                                                                                                                                                             <  2:111062/226‑1 (MQ=255)                                                                                             gAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGCTATTGCTGAGTCCACCCGCCGTATTGCGTCAAGTCGCATTcc                                                                                                                                                  >  1:437118/1‑251 (MQ=255)                                                                                                   gCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCg                                                                                                                                                                                                         >  1:167972/1‑190 (MQ=255)                                                                                                   gCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCg                                                                                                                                                                                                         <  2:167972/190‑1 (MQ=255)                                                                                                          tAAAACATGTCCTGCATGGCATTAGTTTGTTGGGACAGTGCCCAGATAGCCTCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTCGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCCAGTCCCATTCCGGCTGATCACATg                                                                                                                                     >  2:71102/1‑251 (MQ=255)                                                                                                                   tCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGTTCCCATGGTGCTggtg                                                                                                                            >  2:495253/1‑248 (MQ=255)                                                                                                                   tCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTCCCGTGGCGAGAAACTATCGATCGCCATTATGACCGGCGTCTTAGCAGCGCGCGGTCACACCGTTACCGTTATCGATCCGGTCGACAAAATGCTGGCAGT‑GGGGCATTACCTCGAATCTACCATCGCTATTGCTGAATCCACCCGCCGAATTGCGGCCAGTCGCATTACGGCTGATCACATggccctgctg                                                                                                                            >  1:228718/1‑243 (MQ=255)                                                                                                                      tGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGc                                                                                                                          >  1:428158/1‑250 (MQ=255)                                                                                                                        cATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAgg                                                                                                                       >  1:196053/1‑251 (MQ=255)                                                                                                                             cATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGCTATTGCTGAGGCCACCCGCCGTATTGCGGCCAGTCGCATTCCGGCTGATACCATGGTGCTTATGGCCGGTTTc                                                                                                                   >  2:339316/1‑250 (MQ=255)                                                                                                                                   tttGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCAccgccg                                                                                                            >  1:396160/1‑251 (MQ=255)                                                                                                                                               gTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGCTATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGATGATGGAAGGTTTCACCGCCGGTAATGacaacg                                                                                                >  2:408436/1‑246 (MQ=255)                                                                                                                                                    cGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGa                                                                                            >  1:121641/1‑250 (MQ=255)                                                                                                                                                            aTCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGAGTATTAGAAGCGCGCGGTCACAACGATACCGTTATAGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGAATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCACATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTCATGAAAAAGGCGAACtgggggg                                                                                   >  1:506650/1‑247 (MQ=255)                                                                                                                                                               aaCGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTggtggtact                                                                                >  1:297945/1‑248 (MQ=255)                                                                                                                                                                   ctgcgctgATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCgg                                                                                                           >  1:100203/1‑220 (MQ=255)                                                                                                                                                                    tgcgctgATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGTCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGAc                                                                           >  1:43447/1‑251 (MQ=255)                                                                                                                                                                                    ggCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGCTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGCGCTGATGGCAGGGTTCACCGCCGGTAATGAAAAAGGCAGACTGGTGGGACTTGGACGCACCGGTTCCGacc                                                            >  2:482177/1‑249 (MQ=255)                                                                                                                                                                                           aaaTGACGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTAACACCGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCCAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGACAAAGGCGAACTGGTGGGACTTGGACGCAACGGTTCCGACTActccgcg                                                    >  2:17594/1‑247 (MQ=255)                                                                                                                                                                                             aTGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGTTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTActccgcggc                                                  >  1:195314/1‑245 (MQ=255)                                                                                                                                                                                                 cGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGCTCACATGGTGCTGATGGCAGGTTTCCCCGCCGGGAATGAAAAAGGCGAACTGGTGAGACTTGGTCGCAACGGTTCCGACTACTCCGcggcgggg                                              >  2:273019/1‑249 (MQ=255)                                                                                                                                                                                                   aTCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAAATGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGgct                                            >  2:74656/1‑251 (MQ=255)                                                                                                                                                                                                                            aaGCGCGCAGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCTGGTTTCACCGCCGGTAATGAAAAAGGCGCACTGGGGGTCATTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGGCGGCTGCCTGTTTTAGCGCCGGttgt                   >  2:283972/1‑251 (MQ=255)                                                                                                                                                                                                                             aGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCgg                                                                                                                                                >  1:211983/1‑125 (MQ=255)                                                                                                                                                                                                                                             gTTACCGTTATCGATCCGGTCGAAAAACAGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCAGACTACTCCGCGGCGGTGCTGGCTGACTGTTTTCGGGACGATTGTTGCGAGGTTTggacggg  >  1:107071/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       ATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGT‑GGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGA  >  NC_000913/533‑1019

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