Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,359	C→T	I377I (ATC→ATT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,359	0	C	T	97.7%	124.9 / NA	43	I377I (ATC→ATT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/26);  minor base G (0/1);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATG  >  NC_000913/4619121‑4619605                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          cgAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTc                                                                                                                                                                                                                                               >  2:479376/1‑251 (MQ=255)  gAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg                                                                                                                                                                                                                                              <  1:228064/251‑1 (MQ=255)     aTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAAACGACTGCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCt                                                                                                                                                                                                                                           <  2:372800/251‑1 (MQ=255)          gcgcagTGCTGGACAACGGTCAAGCGGCGGAAGTCTTTGGTCGTATGGTAGCGGAACAACAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCAGACAGCGCTGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGTTCCCCGCGCGCTGGGGGTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCCTCTGACCCCCTTGATTACaacggagcctttctt                                                                                                                                                                                                                                      >  2:103888/6‑237 (MQ=255)                         aCGGAAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAACAGGCACAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGATGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGTAATGGATACCCGCCCGCTGGGGATGGGAGTGTTTGCAATGGGGGGCGGTCGGCGTCGAGCATCTGACACCATTGATTACAGCGACGGCTTTACTGATATGGCGCGtttg                                                                                                                                                                                                                       >  2:355067/1‑248 (MQ=255)                             tAAAGCGGCCGAAGGCGTTGGTCGAATGGGAACGGCCCAAAAAGGCCCCACCGACTTCGTTGAGAACGCAGCGAAGTATCAGCAGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACACGCGCGATGGGGATGGCAGTGGTTGCCATGGGCGGCGGTCTACTTCCGGCATCTGCCACCATTGTTTACAGCGTCGGCGTTACTGCTATGGCGCGTTTggggg                                                                                                                                                                                                                   >  2:492891/1‑249 (MQ=255)                                              ttGGCCGTCTGGTACCGGTACAAAAGTGCCCGACCGACTTCATGGAGACCTCCGCGAATCACCTGACGTCCGCCATGCTGACGAACGCAGTCTATGCTCATACCGAAGGTTTTGTCAGGGAATCGGTTACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCTTCGTCAGGCATCTGACGCCATTGATTACAGCGTCGGATTTACTGATATGGCGCGTTTGGGCGagcattgagccggtca                                                                                                                                                                                                   <  2:62730/250‑16 (MQ=255)                                               tGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAg                                                                                                                                                                                                  <  1:274465/250‑1 (MQ=255)                                                           aGCGGCACAAAAAGGCCCAACCGGCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCg                                                                                                                                                                                      <  1:379803/250‑1 (MQ=255)                                                           aGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGTCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCg                                                                                                                                                                                      <  1:77781/250‑1 (MQ=255)                                                              ggCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTa                                                                                                                                                                                  <  1:19236/251‑1 (MQ=255)                                                                  caAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACACGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGCCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCATTACTGATATGGCGCGTTTGGGCGGCGAGGTAGACGGGCAGCGTCCGCTGGCGGTTaaccc                                                                                                                                                                              >  2:193677/1‑247 (MQ=255)                                                                       aGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAAACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTGCGCTGGGGGGTATctcaccgg                                                                                                                                                                         >  2:47189/1‑244 (MQ=255)                                                                        ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCAcgcgc                                                                                                                                                                         >  1:476055/1‑249 (MQ=255)                                                                        ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGGTATCCAcgcgg                                                                                                                                                                         >  2:437241/1‑249 (MQ=255)                                                                            caaccgacTTCGTTGAGAACTACGCGAATTATCTGCCGACAGCGATGCGTACGAAAGCAGTCTATGCTGATACCGAATGTTTGGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTTCTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGAc                                                                                                                                                                   <  2:406831/249‑1 (MQ=255)                                                                                        gtGAGAACTACGCGAAGTATCGGCCGACAGCGATGATGACGAAAGCAGTATATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGACAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTCCGCTGGCGGTTATCCACGCT‑AAAGACGAAAACAGCTgg                                                                                                                                                       <  2:434412/250‑1 (MQ=255)                                                                                             aaCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAgg                                                                                                                                                   >  1:214521/1‑250 (MQ=255)                                                                                               cTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCCCG‑AAAGACGAAAAAAGCTGGCAGGa                                                                                                                                                  >  1:186712/1‑249 (MQ=255)                                                                                                 aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAg                                                                                                                                                <  1:495334/249‑1 (MQ=255)                                                                                                 aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAgcg                                                                                                                                              <  2:82640/251‑1 (MQ=255)                                                                                                        gtttcgtcctgcagtgttGCTGACGAAAGCAGTCTAGGCTTATGCAGAATGTTGTGGCAGTGAAATGGAGACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGGGGCGGGCGTCGTCAGGCATCTTCCTCCAGTGATTTCAGCTGCTGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGGCCGCTGGCGGATATCCACGCG‑AAAGACTAAAACAGCTGGCAGGAAgcggcgcaa                                                                                                                                        <  2:441375/233‑4 (MQ=255)                                                                                                              gCCGACAGCGATGCGGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGt                                                                                                                                   <  1:350867/249‑1 (MQ=255)                                                                                                                      agaggCTGACGAAAGCAGTTTCTGCTGAGACCGAAGGTTGTGGCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATGGATTACAGCGTCGGCTTTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGAAGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGCGA‑CAGCGGCa                                                                                                                        <  2:97849/247‑1 (MQ=255)                                                                                                                         tGCTGAAGAAAGCAGTCTATGCTGAAACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGAAGTGGTTGCACTGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTGACCACGCGGAAAG‑CGAAAACAGCTGGGAGGGAGCGGCGGAGACGGGGA‑AAGCAGCAAtt                                                                                                                     >  1:181723/1‑251 (MQ=255)                                                                                                                              aCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACt                                                                                                                <  2:81775/251‑1 (MQ=255)                                                                                                                                        tCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGATATCCAAGCG‑ACAGCCCAAAACCACCGGCAGGCACCGGCGCAAGCGGGGAGAAG‑GGCAATTAAACTTTCCGAtca                                                                                                       >  1:364682/1‑248 (MQ=255)                                                                                                                                             gCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCAcc                                                                                                 <  1:437241/251‑1 (MQ=255)                                                                                                                                                    ccGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGc                                                                                          <  1:502749/251‑1 (MQ=255)                                                                                                                                                        agtttttGTCAGGGAAAAGGATACCCGCGCGCTGGGGAGGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCaca‑c                                                                                     <  2:366105/248‑1 (MQ=255)                                                                                                                                                          ggtttGGCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCa                                                                                   <  1:496009/249‑1 (MQ=255)                                                                                                                                                               gTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTg                                                                               <  1:199447/250‑1 (MQ=255)                                                                                                                                                                       aaTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTGTCCACCCG‑GAAGACGAAAAAAGCTGGGAGGAAAAGGCGAAAACGGGGA‑AAGGCGCAATTAAAATTTCCGATTAAACACCGGaaaagacacaaaattgctattccc                                                                       >  1:503623/1‑230 (MQ=255)                                                                                                                                                                                            gATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGA‑AAGCGGCAAATAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAAATGTCTATCGCCGTATCAGCGAATAACggg                                                    >  1:424322/1‑247 (MQ=255)                                                                                                                                                                                                                 cgtcggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAttt                            <  1:120460/241‑1 (MQ=255)                                                                                                                                                                                                                   gcgggcgtcGTCAGGCATCTGAAACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACT‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGa                          <  1:43828/243‑1 (MQ=255)                                                                                                                                                                                                                    gggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGTTATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTGTCTATCGTCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAGATTTGaa                         <  2:86294/244‑1 (MQ=255)                                                                                                                                                                                                                    cggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGGGGGGA‑AAGAGGCAATTAAACCTGCCGATAAAACACCGGAAAG‑ACAACCAACTGTCTATCGCCGCATCACAGAAACACGGGATACTGATCTGATCATTTAAAtttta                          >  2:73279/8‑248 (MQ=255)                                                                                                                                                                                                                        cgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCa                      <  1:427915/247‑1 (MQ=255)                                                                                                                                                                                                                        cgtcGTCAGGCAGCTGACACCATTGATTACAGCGTCGGCGTTACTGATAGGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACT‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCAc                     <  2:205469/248‑1 (MQ=255)                                                                                                                                                                                                                            gtcgggCACCTGACACCATTGATTACAGCGTCGGCTGTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCCCTGAg                 <  2:53276/247‑1 (MQ=255)                                                                                                                                                                                                                                       gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACa      <  1:479376/251‑1 (MQ=255)                                                                                                                                                                                                                                           ccATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCAGCG‑AAAGACGAAAACAGCCGGCAGGAAGCGGCGGAAGCGGGGG‑ACGAGGCCATTAAACTTGCCGCTAAAACACCCGGAAGCACA‑CCAACTGGCTATTGCCCTATCAGCGAAAAACGGGATACTGATCTGATTCTTTAAATTTTAAGCACTgtggaccgggaaaatatt  >  2:175285/1‑234 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          CGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATG  >  NC_000913/4619121‑4619605

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