Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	674,518	G→A	A89V (GCG→GTG)	leuS ←	leucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	674,518	0	G	A	94.8%	180.2 / NA	59	A89V (GCG→GTG)	leuS	leucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (36/19);  minor base C (1/2);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTT  >  NC_000913/674287‑674761                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    aGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGAcc                                                                                                                                                                                                                                  >  2:225780/1‑251 (MQ=255)                  cGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCg                                                                                                                                                                                                                  >  1:581385/1‑250 (MQ=255)                    ggCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTAACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATc                                                                                                                                                                                                               >  1:762306/1‑251 (MQ=255)                     gCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGATTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCNGCAGGCAGCCCAAACGCGTCCCAGCCGATCg                                                                                                                                                                                                              >  2:363247/1‑251 (MQ=255)                       aCCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGATCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCCGGCAGACCAAACGCGTCCCAGCCGGTCGGc                                                                                                                                                                                                            >  1:673927/1‑251 (MQ=255)                                             cttATATACCAGGCCTTTTTTATACAGCTCGGTGAAGACTTTCTTTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCCGTCATAACCAAAGCCCCGCATTTTGAGCTGGTTTTTCATCT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGTGCGGTGTTGTTTTTCACCGTCACGCCTTACGCAGGCCGACCAAACGCGTCCCCNCCGCTCGGCTTCAGGACGTTTTTTCCCaccc                                                                                                                                                                                      >  2:389778/1‑248 (MQ=255)                                                     ccAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTgg                                                                                                                                                                              >  2:197744/1‑251 (MQ=255)                                                       aGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGt                                                                                                                                                                             >  1:72218/1‑250 (MQ=255)                                                        ggCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGt                                                                                                                                                                             >  1:576433/1‑249 (MQ=255)                                                          ccTTTTTAAAACAGCTGGGTAAAGAATTTCTGTTCCCACAGGTAGTTTTCCGGGTTCCAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGGCCGGGTTTTTCATAA‑ACCCGATGTTGTCGTCCGTCCACGGGGCCGGGGCGTTGTTGTTTTTCACCGCCCCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCg                                                                                                                                                                         <  1:31651/251‑1 (MQ=255)                                                                ttATACAGCTCGGTGAAGAATTTCTGTTCCCACCGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGGTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGGCCCAGCCGATCGGCTGCAGGAAGTTTTTGCCCAGCATCCGCTGGGAGCGGGGGAt                                                                                                                                                                   >  2:450980/1‑251 (MQ=255)                                                                 tataCAGCTCGGTGAAGAATTTCTGTTCCAAACGGTAGTATTCAGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCACAGCATATTGAGCTGGGTTTTCATAT‑ACGCGATGTTTTCGTACGTCCACGGTGCCGGAGCGGAGTTGTTTTTCACCCCCCCGCCTTCCGCAGGCAGACCAAACGCGTACCAGCCGATCGGCTGAAGGTCTTTTTTGACCCGCATACACTGGTCGCGGGCAATc                                                                                                                                                                  >  2:479334/1‑251 (MQ=255)                                                                      aGCTGGGTGAAGATTTTCTTTCCCCAAGGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCTTTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTc                                                                                                                                                             <  1:523355/251‑1 (MQ=255)                                                                          gggtgaaataTTTCTGTTCCCAGCGGTATTATTCCGGGTTACGGTTTGCCGGCTCGCGGCCCCAGTCATAACCAAAGCCCAGCAGTTTGGGCTGGTTTTCCAAAT‑CCGCGTGGTTGTCGTACTCCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCg                                                                                                                                                         <  2:536008/242‑1 (MQ=255)                                                                          cGGTGAAGAATTTCTGTTCCCACCGGTAGTATTCCGGCGTACAGTTTGCCAGCCCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGTTGTTGTCGTACGTCCACGTTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCAcc                                                                                                                                                          <  1:489726/250‑1 (MQ=255)                                                                               aagaaTTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAACGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCCGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGCTCACGTCACCGCTggg                                                                                                                                                    >  2:618378/1‑250 (MQ=255)                                                                               aagaaTTTCTGTCCCAAACGGTAGTATCCCGGCTTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGGGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGgt                                                                                                                                                    <  1:197744/251‑1 (MQ=255)                                                                                 aaaTTTCTGTCCCCAACGGTAGTATCCCGGCGTCCAGGTTGCCAGCTCGCGGCTCCTGTCAAACCCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑CCGCGATTTTGTCGTACGTCCACGGTGCCGGAGCGGTTTTGTTTTCCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGgtgt                                                                                                                                                  <  1:328506/250‑1 (MQ=255)                                                                                  aaTTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATCCGCTGGTAGCGGGCGATCACGTCACCGATGgtga                                                                                                                                                  >  1:690954/1‑249 (MQ=255)                                                                                      tCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCAGGTAGCGGGCGATCACGTCACCGATGGTGTAGtt                                                                                                                                              >  2:748913/1‑250 (MQ=255)                                                                                            cccAACGGTAGTTTTCGGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACAAAGGCCGAGCATTTTGAGCTGTTTTTTCATAT‑ACGGGATGTTGTCGTGCGTCCACGGTGCGGGAGGGGTTTTGTTTTTCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTtacgtacg                                                                                                                                       <  1:363247/251‑1 (MQ=255)                                                                                              caacggttatgccccGGCGTACGGGTTGCTAGACCGGGGCTCATGTCAAAACCAAAGCCCAGCGTTGTGGTCTG‑TTTTTTAGATTGCGCGTTGTTGTCGTACTCCCCCGGTGCCGGGGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTg                                                                                                                                     <  2:135780/238‑1 (MQ=255)                                                                                                     agtaTTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCAtg                                                                                                                              >  1:401161/1‑251 (MQ=255)                                                                                                      gtaTTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCAtgg                                                                                                                             >  1:488709/1‑250 (MQ=255)                                                                                                            cGGCGTACAGGTTCCCACCTCGCGGCCCCATTCATAACCAAAGCCCGGCATTTTGAGCTGGTTTTTCATAT‑CCGCGATGTTGTCGTACGTCCACGGTGCCGGGGGGGTGTTGTTTTTCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTc                                                                                                                       <  1:536290/251‑1 (MQ=255)                                                                                                                   cAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGATTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGAcc                                                                                                                   >  1:103510/1‑248 (MQ=255)                                                                                                                               ccgcgGCTCCAGTCATAACCAAAGCCCCGCTTTTGGAGCTGGTTTTTCATAT‑CCGCGATGTTGTCGTCCGTCCAGGGTGCGGGAGCGTTGTTGTTTTTCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAg                                                                                                    <  1:748913/250‑1 (MQ=255)                                                                                                                                        cAGTCATAACAAAAGCCCAGCTTTTGGGGCGGTTTTTTCAGAT‑ACGGGGTTTTGTGGTCCCCCCACGTGGCGGGAGGGGTTTTTTTTTCCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAg                                                                                           <  2:576433/251‑1 (MQ=255)                                                                                                                                            cATAACCAAAGCCAAGCATTTGGAGCTGGTTTTTCATAT‑CCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAg                                                                                       <  1:225780/251‑1 (MQ=255)                                                                                                                                                aCCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAAt                                                                                   >  1:492006/1‑251 (MQ=255)                                                                                                                                                    aaGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCCGAAGGATAGGGAAGCATAGACAGGCAGTAATACtt                                                                               >  2:211283/1‑251 (MQ=255)                                                                                                                                                     aGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCCGGACGTTTTTGCCCAGCATACGCTGGAAGCGGGCGCTCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGAGTAGTCGACCAGAAGNATAGGGAAGCATAGTCAGGCAGTCATACTtc                                                                              >  2:371492/1‑251 (MQ=255)                                                                                                                                                      gCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTtct                                                                             >  1:226454/1‑251 (MQ=255)                                                                                                                                                         cAGCTTTTGGAGTTGGTTTTTCAATT‑CCGCGATGTTGTCGTCCGTCCAGGGTGCCGGGGCGGTGTTGTTTTCCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCtt                                                                          <  1:704506/251‑1 (MQ=255)                                                                                                                                                          aGCATTTTGAGCGGGTTTTTCATAT‑GCGCGATGTGGTCGTCCGCCCCCGGTGCCGGAGCGGTGTTGTTTTCCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCttt                                                                         <  1:450980/251‑1 (MQ=255)                                                                                                                                                               tttGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGctct                                                                    >  1:197475/1‑251 (MQ=255)                                                                                                                                                                  gAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGt                                                                 >  1:62215/1‑251 (MQ=255)                                                                                                                                                                  gAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGt                                                                 >  1:354524/1‑251 (MQ=255)                                                                                                                                                                  gAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGGGTAGTTACGTACGTGGCCCATGTGGAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGt                                                                 >  2:574422/1‑251 (MQ=255)                                                                                                                                                                         tttttCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATCCGCTGGTAGCGGGCGATCACGTCACCGGTGGTGTAGTTACGTACGTGGCCCCTGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACTGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGt                                                          >  2:58110/1‑251 (MQ=255)                                                                                                                                                                              cATAT‑ACGGGATTTTTTCGTACGCCCACGGTGCCGGAGCGGTTTTGTTTTCCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACt                                                      <  1:223352/250‑1 (MQ=255)                                                                                                                                                                               atat‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTc                                                    >  1:608454/1‑251 (MQ=255)                                                                                                                                                                                 at‑aCGCGATGTTGTCGTACGTCCAGGGTGCCGGAGCGGTGTTGTTTTCCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCa                                                   <  1:618378/250‑1 (MQ=255)                                                                                                                                                                                 at‑aCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCa                                                   >  1:267577/1‑250 (MQ=255)                                                                                                                                                                                       cgATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGCCAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAAtgtg                                             >  2:269183/1‑251 (MQ=255)                                                                                                                                                                                       cgATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAAtgtg                                             >  1:607921/1‑251 (MQ=255)                                                                                                                                                                                       cgATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAAtgtg                                             >  1:502851/1‑251 (MQ=255)                                                                                                                                                                                        gATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGCTAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAAtgtgg                                            >  1:46220/1‑250 (MQ=255)                                                                                                                                                                                                 gTACGTCCACAGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCAt                                   >  1:671163/1‑251 (MQ=255)                                                                                                                                                                                                  ccctcccacgttgccggttgtgtgtTTTTTTTGACCGCCACGCCTTCCGCATGCAGACCAAACGCGTCCCAGCCGATCGGTTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGAGGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGAGAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCGCTCTTCGTCGGTTACTTCTAATGTGCTCTTCTCATc                                  <  1:538217/230‑1 (MQ=255)                                                                                                                                                                                                    cgcccacgggtgcggggCGGTTTTGTTTTTCCCCGCCCCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTGGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATccc                                <  1:69166/235‑1 (MQ=255)                                                                                                                                                                                                                  gAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACt                  >  1:322572/1‑251 (MQ=255)                                                                                                                                                                                                                    gCGGTTTTGTTTTTCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACttt                <  1:520996/251‑1 (MQ=255)                                                                                                                                                                                                                       ggtttgttTTTCACCGCCACTCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGa             <  1:479334/248‑1 (MQ=255)                                                                                                                                                                                                                       gggttgttTTCCCCCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCGTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGa             <  2:607921/249‑1 (MQ=255)                                                                                                                                                                                                                         gttgttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGc                                                                                         >  1:563397/1‑173 (MQ=255)                                                                                                                                                                                                                         gttgttTTCCACCGCCACGCCTTCCGCAGGGAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGAt            <  1:195188/250‑1 (MQ=255)                                                                                                                                                                                                                            gttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGCAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCATCGGTTACTTCAAATGTGCGCTTCTCCTCCCAATGAAGCTGTACTTTGGATTCt         >  2:40059/1‑250 (MQ=255)                                                                                                                                                                                                                                ttccccgccaCGCCTTCCCCAGGCAGACCAACCGCGTCCCAGCCGATCGGCTGCAGGTCGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTTTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTAtct     <  2:525878/246‑1 (MQ=255)                                                                                                                                                                                                                                  cACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCtt  >  1:457823/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    AGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATAT‑ACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTT  >  NC_000913/674287‑674761

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