Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,393	C→T	L389L (CTG→TTG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,393	0	C	T	96.5%	181.0 / NA	58	L389L (CTG→TTG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/34);  minor base G (0/2);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCA  >  NC_000913/4619147‑4619635                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      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)                                                                                                                                                                                      gcggtgggcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGTGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTAGCAGCGAATAACGGTATACTGATCTGAGCATTTAAAtt                                                           <  2:311984/240‑1 (MQ=255)                                                                                                                                                                                      gcggcgGTCGGCGGCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAtt                                                           <  2:388420/251‑1 (MQ=255)                                                                                                                                                                                            gtcgtctgcAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc                                                     <  2:198413/243‑1 (MQ=255)                                                                                                                                                                                                     aGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCCCGGAAGACGAAAAAAGCTTGCCGGGAGCGGCGAAAGCGGGG‑AAAGGGG‑CGAATAAA‑ATTTCCGATAAAGCACCGGAAAGCACACCAAATGTCTTTCGCCGATTCCGCGAATAACGGGTTACTTATCTGATACTTTTAATTTGTAGCaatgtgaac                                            >  2:57398/1‑243 (MQ=255)                                                                                                                                                                                                     aGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTAc                                            <  1:57398/251‑1 (MQ=255)                                                                                                                                                                                                             gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGGGGACGGAGAAc                                     >  1:153723/1‑250 (MQ=255)                                                                                                                                                                                                             gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                                        <  2:6987/247‑1 (MQ=255)                                                                                                                                                                                                                      gATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACg                           <  1:125511/251‑1 (MQ=255)                                                                                                                                                                                                                         tACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTg                         <  1:248331/250‑1 (MQ=255)                                                                                                                                                                                                                                  ggCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTtattat               <  2:153723/251‑1 (MQ=255)                                                                                                                                                                                                                                        aCTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATACACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACAACAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAAAGTGCATTTATTATGGTGCt         >  1:21981/1‑251 (MQ=255)                                                                                                                                                                                                                                               gggCGCGTTTGGGCGACCAGTGAGACGGTCAGCGTCCGCTGGCGGATAGCCTCGCGAAAGCCGAAAACAGCTGGCAGGAGGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAGTTGAAGCACTGAGTTCGGAGAACATATGAAACGTGCATTTCTTATGGTGCTGGACTCa  <  2:69636/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCA  >  NC_000913/4619147‑4619635

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