Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,338	C→T	R370R (CGC→CGT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,338	0	C	T	96.4%	162.1 / NA	55	R370R (CGC→CGT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/32);  minor base G (0/2);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATT  >  NC_000913/4619097‑4619578                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    ggAAAGACGACGACGAAGCGCTCGGGAAAGTGCAAGCAGTGCTGGACAACGGTAAAGCGGAAGAAGTCGTTGGTCGGATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTAGGCTGATACCGAATGTTTTGTCAGTGAAATGGATACCCGCGCGCCGTGGATGTCAGTGGTCGCAATGTGCGGCGGTCGTCGTCAGGCa                                                                                                                                                                                                                                            <  2:71509/250‑1 (MQ=255)    aaGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCt                                                                                                                                                                                                                                         <  1:365441/251‑1 (MQ=255)                 aGCGCGCGCAAAATTGCAAGCTATTCTGGACAACGGTAAAGTGGCACAAGTCTGTGGTCGTATGGTATCGGCTCAAAAAGGCCTAACCGACTGCGTTGAGAACCACGCGACGTATCTGCCGACAGCGATGCGGACGAAAGCAGTCTATGCTTATCCCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGAtt                                                                                                                                                                                                                            <  2:111905/251‑1 (MQ=255)                     cgcgcgAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAg                                                                                                                                                                                                                        <  2:478238/251‑1 (MQ=255)                             aTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGATGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTAAAGCGTCGGCt                                                                                                                                                                                                                >  1:427992/1‑251 (MQ=255)                               tgcaagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCtt                                                                                                                                                                                                               <  2:280967/242‑1 (MQ=255)                                   agcagTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTg                                                                                                                                                                                                          >  1:337401/5‑251 (MQ=255)                                            ggACAACGGTAAAGCGGCAGAAGTCGTTGGTCGTAGGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAGACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCGGACACCATTGATTACAGCGTCGGCTTTACTGATATGgcg                                                                                                                                                                                                  <  2:237042/250‑1 (MQ=255)                                                   ggTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTggg                                                                                                                                                                                          >  1:99819/1‑251 (MQ=255)                                                        aGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGAcc                                                                                                                                                                                     >  1:375919/1‑251 (MQ=255)                                                                 aGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACg                                                                                                                                                                            >  1:41726/1‑251 (MQ=255)                                                                                           aaaaaGGCCCAACCGTCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCat                                                                                                                                                  <  2:93133/251‑2 (MQ=255)                                                                                             aaaGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTCGACGGTccgc                                                                                                                                                                      >  2:212870/1‑226 (MQ=255)                                                                                             aaaGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGc                                                                                                                                                                      <  1:212870/229‑1 (MQ=255)                                                                                             aaaGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTACGCTGGCGGGTATGCAccg                                                                                                                                                >  2:42325/1‑249 (MQ=255)                                                                                                  cccaaccgacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCCGCGTCCGCTGGCGGtttttccccgcgaa                                                                                                                                            >  2:178473/5‑238 (MQ=255)                                                                                                  cccaaccgacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAg                                                                                                                                           <  1:407745/247‑1 (MQ=255)                                                                                                          acTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCCGCGTCCGCTGGCGGTTTTCCCCGCGAAAGACgcaac                                                                                                                                    >  2:84678/1‑246 (MQ=255)                                                                                                          acTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAaac                                                                                                                                   <  1:178473/251‑1 (MQ=255)                                                                                                             tCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAaacag                                                                                                                                 <  2:337401/250‑2 (MQ=255)                                                                                                                    gaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAgg                                                                                                                         <  1:172383/251‑1 (MQ=255)                                                                                                                             gAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCTCGCGCAAGACGaaaaacgccggcaggaaacggcgg                                                                                                                >  1:10128/1‑231 (MQ=255)                                                                                                                               aGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGaaa                                                                                                              <  1:264142/251‑1 (MQ=255)                                                                                                                                gTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCGCTAAAGACGAAAACAGCTGGCAGGAAGCGGCGaa                                                                                                               >  1:307633/1‑249 (MQ=255)                                                                                                                                         gACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGaa                                                                                                               <  1:273906/240‑1 (MQ=255)                                                                                                                                               gATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCCCGCGAAAAACGACAACACCTGGGAGGGAACGGCGAAAGCGGGG‑AAAACGG‑Ca                                                                                             >  1:303129/1‑250 (MQ=255)                                                                                                                                                          aaGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGcc                                                                                <  1:325780/251‑1 (MQ=255)                                                                                                                                                          aaGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGACAGCGTCCGCTGGCGGTTATACAAGGGAAAGAAGAAAACAGCGGGCAGGAGGCGGAGAAAACGGGG‑AAAGCGGGCAATTAAA‑‑TTGcc                                                                                >  1:246717/1‑251 (MQ=255)                                                                                                                                                               gTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCTCGCTGGGGATGGCATTGGGTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCATGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATaa                                                                           <  2:398955/251‑1 (MQ=255)                                                                                                                                                                    ggcggATACCGAAGGTTGGGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGTCGGCAGGCATCGGACACCATTGATTACATCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCAc                                                                      <  2:363682/247‑1 (MQ=255)                                                                                                                                                                    ggCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCAc                                                                      <  2:162842/250‑1 (MQ=255)                                                                                                                                                                        gATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCg                                                                    <  1:483270/249‑1 (MQ=255)                                                                                                                                                                         aTACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATCCACGCGAAAGACGGAAACAGCTCGCAAGAAGCGGCGACAGCGGTGGAAAG‑GG‑CAATTAAAGCTT‑CCGTTAAAACACCGGa                                                                  >  1:242313/1‑250 (MQ=255)                                                                                                                                                                          ggcagAAGGTTTTGTCTGTGAAATGGATACCCGCGCGCTGGTGATGGCAGGGGGTGCAATGGACGGCGGTCGTCGTCAGGCATCTGACACCATTGATGACAGCGTCGGCGTTACTGATATGGCGCGTGTGGGAGATCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCATGAAGCGGCGAAAGCGGTG‑AAAGTGG‑CAATTAAA‑CTTGCCGATAAAGCGCCGGaca                                                                <  2:428508/247‑3 (MQ=255)                                                                                                                                                                           aCCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGaaa                                                                <  1:93898/250‑1 (MQ=255)                                                                                                                                                                             cGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTTTCCCCCCGAAAGACGGAAACCGCTGGCCGGGAGCGGCGAAAGCGGGG‑GAAGCGG‑CAAATAAA‑CTTGCCGATAAAGCCCCgggaaacc                                                             >  1:311984/1‑245 (MQ=255)                                                                                                                                                                             cGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCGGGCGGGTATCCCCGCGAAAGAAGAAAACAGCCGGGCGGAAGCGGCGGAAGGGGTG‑AAAGGGG‑GAATAAAA‑CTTTACGATAAAACACCGGaaaaa                                                              >  2:425870/1‑248 (MQ=255)                                                                                                                                                                                  gTTTTGTCAGTGAAATGGATACCCGCGCGTTGGGGATGTCAGTGGGTGCAATGGTCGGCGGGCGGTGGCAGTCATCTGACACCAGTGATTACAGCGTCGGCTTACCTGATACGGCGCGTGTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGCTATTCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCtcacta                                                        <  2:133399/251‑3 (MQ=255)                                                                                                                                                                                    tttGGCAGGGAAAGGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCaa                                                       <  2:388590/250‑1 (MQ=255)                                                                                                                                                                                          aGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGc                                                                                                                            <  2:471459/178‑1 (MQ=255)                                                                                                                                                                                           gTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAGAAAAAAGCTGGCAGGAAGCGGCGAAAGAGGGG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAACCCACCAACTGTccat                                               >  1:80098/1‑248 (MQ=255)                                                                                                                                                                                             gAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGG‑GAAGCGG‑CAATTAAA‑CTTTCCGATAAAGCAACGGGAAGCACCCCAACCGTCtttcc                                             >  1:242568/1‑247 (MQ=255)                                                                                                                                                                                              aaaTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCCACCGTCTATccc                                            >  1:6631/1‑249 (MQ=255)                                                                                                                                                                                                aTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCg                                          <  2:451083/251‑1 (MQ=255)                                                                                                                                                                                                  ggATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTa                                        <  1:446711/251‑1 (MQ=255)                                                                                                                                                                                                       cccGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCTTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGc                                   <  2:435103/251‑1 (MQ=255)                                                                                                                                                                                                           cgggcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGACAGGCATCTGACACCATTGATTACAGCGTCTGCTTGACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGGCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGGCTATCGCCGTATCAGCGAAt                               <  2:90088/248‑1 (MQ=255)                                                                                                                                                                                                             cgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGcc                                           <  2:185159/237‑1 (MQ=255)                                                                                                                                                                                                                 ggggATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCAGTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCTAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                           <  2:459529/249‑1 (MQ=255)                                                                                                                                                                                                                     aTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCCCGGAAGAACAAAACAGCTGGGAGGAAGCGGCGGAAACGGGG‑AAAAGGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAAAAAACCCACTGGCTATCCCCGTATCaccggcaaaagggatac                     >  1:311959/1‑235 (MQ=255)                                                                                                                                                                                                                          aGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatc                <  2:6631/251‑1 (MQ=255)                                                                                                                                                                                                                            tGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCCCGCGAAAGACGAAAACAGCCGGCCGGAAGCGGCGAAAGCGGGG‑AAAACGG‑CCATTAAA‑CTTGCCCATAAAGCACCGGGAAGCACACCAACCGTCCATCGCCGGATCAGCGGATAACGGGATAccgatctg              >  1:388420/1‑251 (MQ=255)                                                                                                                                                                                                                                tGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCAGTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGa                                                                                                                                          <  1:24359/126‑1 (MQ=255)                                                                                                                                                                                                                                    gtggtgggcggTCGTCGTCCGGCATCTGACACCATTGAATACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAAGAACTTTATACTGATCTGATCATTTa      <  2:244998/245‑1 (MQ=255)                                                                                                                                                                                                                                        gcggcgGTCGGCGGCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAtt  <  2:388420/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    CGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATT  >  NC_000913/4619097‑4619578

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