Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,335	A→T	G369G (GGA→GGT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,335	0	A	T	96.4%	166.1 / NA	56	G369G (GGA→GGT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (21/32);  minor base G (0/2);  total (22/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

GGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATT  >  NC_000913/4619087‑4619578                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       ggCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGAAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCAAAAAAAGGCCCAACCGACTTCGTGGAGAACTACGCGAAGAATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGAAACCGAAGGTTTTGTAAGGGAAATGGATACACGCGAGCTGGGGTGGGCAGAGGTAGCAATGGGGGGCGGACg                                                                                                                                                                                                                                                      >  2:424745/1‑251 (MQ=255)           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)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       GGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATT  >  NC_000913/4619087‑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: ‑