Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,142,598	G→A	L117L (CTC→CTT)	udk ←	uridine/cytidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,142,598	0	G	A	98.6%	261.9 / NA	74	L117L (CTC→CTT)	udk	uridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (37/35);  minor base C (1/0);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACC  >  NC_000913/2142352‑2142840                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   gATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCaaga                                                                                                                                                                                                                                                 >  2:160857/1‑250 (MQ=255)  aTGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAat                                                                                                                                                                                                                                               >  1:459430/1‑251 (MQ=255)   tGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCACCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAatg                                                                                                                                                                                                                                              >  1:162476/1‑251 (MQ=255)     cGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGGTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAatgat                                                                                                                                                                                                                                            <  1:656054/251‑1 (MQ=255)       gTTTTTCCCGCCGCGCGGCACTATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGAGGCGGCGCATCAGGCAGATAGCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCGTCAAGAAtgatga                                                                                                                                                                                                                                          <  2:747646/251‑1 (MQ=255)           ttCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACctt                                                                                                                                                                                                                                      >  1:445390/1‑251 (MQ=255)            tCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGACTTGCAGGCACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGCATGATGACcttc                                                                                                                                                                                                                                     >  1:103001/1‑251 (MQ=255)            tCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACctt                                                                                                                                                                                                                                      <  1:7847/250‑1 (MQ=255)                         cACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCg                                                                                                                                                                                                                        <  2:459430/251‑1 (MQ=255)                            gATAATGTCCGCATATTGTTTAGAAGGATCAATGAATTGCTGGAACATCGGGCGCAAGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTGCTTCGGCTCAACCGTCa                                                                                                                                                                                                                     <  2:178099/251‑1 (MQ=255)                             aTAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCAc                                                                                                                                                                                                                    <  1:708194/251‑1 (MQ=255)                                aTGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGt                                                                                                                                                                                                                 <  1:5917/251‑1 (MQ=255)                                       cATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctttc                                                                                                                                                                                                          >  1:326781/1‑251 (MQ=255)                                       cATATTGGTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctttc                                                                                                                                                                                                          <  2:243226/251‑1 (MQ=255)                                           ttGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAACTGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAc                                                                                                                                                                                                      >  2:651538/1‑251 (MQ=255)                                             gTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAcg                                                                                                                                                                                                     <  1:160857/250‑1 (MQ=255)                                               ttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGt                                                                                                                                                                                                  >  2:223999/1‑251 (MQ=255)                                                     ggCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTc                                                                                                                                                                                            >  2:330883/1‑251 (MQ=255)                                                         cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCaa                                                                                                                                                                                          <  1:454670/249‑1 (MQ=255)                                                          aaTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACAt                                                                                                                                                                                       >  2:350689/1‑251 (MQ=255)                                                               aTTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTg                                                                                                                                                                                  >  1:499011/1‑251 (MQ=255)                                                                         cATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGc                                                                                                                                                                         >  1:703312/1‑250 (MQ=255)                                                                            cGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGt                                                                                                                                                                     >  1:111372/1‑251 (MQ=255)                                                                                cgctcGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCGTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAAt                                                                                                                                                                 <  1:473251/247‑1 (MQ=255)                                                                                 gcACGGTTTTTTGATATTGCGCCATCACGGAAGCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAAtt                                                                                                                                                                <  1:507587/251‑1 (MQ=255)                                                                                  cACGGTTTTTTTATATGGCTCCATCACTGAAGCCATTGAACGCCCACGCTCGGTAACGTCACGCTTGATGCGGCGCATCAGGCAGATTTCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTg                                                                                                                                                               <  2:490807/251‑1 (MQ=255)                                                                                    cGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGcc                                                                                                                                                             <  2:326781/251‑1 (MQ=255)                                                                                    cGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCGCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGcc                                                                                                                                                             <  1:208626/251‑1 (MQ=255)                                                                                         ttttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGcc                                                                                                                                                        >  2:700293/1‑251 (MQ=255)                                                                                            tGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCcgcg                                                                                                                                                     >  1:227133/1‑251 (MQ=255)                                                                                             gATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGt                                                                                                                                                    <  2:410254/251‑1 (MQ=255)                                                                                                    cgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACg                                                                                                                                             >  2:588853/1‑251 (MQ=255)                                                                                                            cTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTaa                                                                                                                                      >  1:225556/1‑250 (MQ=255)                                                                                                              gAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTaaa                                                                                                                                     <  1:104155/249‑1 (MQ=255)                                                                                                                 tCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTc                                                                                                                                >  2:504137/1‑251 (MQ=255)                                                                                                                    aTTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCa                                                                                                                               >  1:142032/1‑249 (MQ=255)                                                                                                                      tGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAagc                                                                                                                            <  2:225556/250‑1 (MQ=255)                                                                                                                           gCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGa                                                                                                                      <  1:327693/251‑1 (MQ=255)                                                                                                                              cACGCTCGTTCACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAAAAAAATGCCTTCCAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACAGGCAGGTCAATTGCCGAGCCGCGTTTCACCGCTTGTAAATGCCCAAGCCGCAGCCtg                                                                                                                   >  2:753960/1‑251 (MQ=255)                                                                                                                              cACGCGCGTTAACGGCACGCTGGATTCGGGACACCAGGCAGATATCCCGCGGGGTCACAACGAAAATGGAGAAGTTAAGTTCGTCTCGCAATCGCGCATTCGTCAGCAACAAAATGCCTTCAAGAATGATTACCTTCTTCGGCTCAACCGTCCCCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGTAGGTAAATAGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACtg                                                                                                                   <  2:269265/251‑1 (MQ=255)                                                                                                                                 gCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACtgtg                                                                                                                 >  1:481862/1‑250 (MQ=255)                                                                                                                                    cGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGAt                                                                                                               <  1:22368/249‑1 (MQ=255)                                                                                                                                     gTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACTCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTCTGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGAAGCAGACTGTGATccc                                                                                                            >  1:515374/1‑250 (MQ=255)                                                                                                                                     gTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATcc                                                                                                             >  1:102495/1‑250 (MQ=255)                                                                                                                                               cGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAAGGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGGCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTggg                                                                                                  <  2:742537/251‑1 (MQ=255)                                                                                                                                                 cTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTACACCGGCAGGTCACTTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTgtgtt                                                                                                >  2:296815/1‑247 (MQ=255)                                                                                                                                                   tgaagcggcgCATCAGGCAGATAGCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGAt                                                                                                <  2:241455/245‑1 (MQ=255)                                                                                                                                                     aTGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTc                                                                                            <  1:350689/251‑1 (MQ=255)                                                                                                                                                                cAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTc                                                                                  <  1:149980/250‑1 (MQ=255)                                                                                                                                                                       atatCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAAcgcg                                                                          >  1:401092/1‑251 (MQ=255)                                                                                                                                                                         atCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGtt                                                                        <  2:162476/251‑1 (MQ=255)                                                                                                                                                                            cAGCGGGGTATTAACGAAATTGGAGAAGTTAAGTTCGTCACTCAAACGCGCATCCGTCAGCAACAAAATGCTTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCGTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGttctt                                                                     <  2:445390/251‑1 (MQ=255)                                                                                                                                                                            cAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACTCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGGGATCCAGCGCGCTGGGAGGGTGGTAGTTGGTCTTAACGCGttctt                                                                     <  2:103001/251‑1 (MQ=255)                                                                                                                                                                                ggggTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCAt                                                                 >  1:433145/1‑251 (MQ=255)                                                                                                                                                                                         aCGAAAATGGAGAAGGTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAAAAAACTGCCTTCAAGAATGAGGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTAGGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATg                                                        <  2:190552/251‑1 (MQ=255)                                                                                                                                                                                           gAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAAAGCGTTCTTCCATCGCCAGATGGc                                                      >  2:418557/1‑251 (MQ=255)                                                                                                                                                                                              aaTGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCaa                                                                                                                              >  2:289278/1‑176 (MQ=255)                                                                                                                                                                                              aaTGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCaa                                                                                                                              <  1:289278/176‑1 (MQ=255)                                                                                                                                                                                              aaTGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCttt                                                   <  2:433145/251‑1 (MQ=255)                                                                                                                                                                                                  gagaAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGGCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATc                                               >  2:697867/1‑251 (MQ=255)                                                                                                                                                                                                          aagttCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGGTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATTGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTaa                                        <  1:753960/250‑1 (MQ=255)                                                                                                                                                                                                           agttCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATa                                      <  1:418557/251‑1 (MQ=255)                                                                                                                                                                                                               cgtaccGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCGTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGCCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGGGAGCCATCGCGCTGGGATGGTCGTATTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAg                                  <  2:703312/246‑1 (MQ=255)                                                                                                                                                                                                               cGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGc                                    >  2:95443/1‑249 (MQ=255)                                                                                                                                                                                                                   aCGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGCTGGCTTTGATCTTTGTAATAGCAGTCtt                              >  2:439785/1‑251 (MQ=255)                                                                                                                                                                                                                          cgcgCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGCCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCggggat                       >  1:564422/1‑248 (MQ=255)                                                                                                                                                                                                                                    tCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAAGCGTATGTTCAACATAGCTGTAAACCGGCAGGGCAATTACCGAGCCGCGTTTCAAAGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTGGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAATCTTCGAGAATTACGCAGAt              >  1:470342/1‑250 (MQ=255)                                                                                                                                                                                                                                     cAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGAtgt            <  1:700293/251‑1 (MQ=255)                                                                                                                                                                                                                                      aGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGCCAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGAtgtg           >  2:487583/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGCTgtgt          >  1:180557/1‑251 (MQ=255)                                                                                                                                                                                                                                        caacaaAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGtt         >  1:679794/1‑251 (MQ=255)                                                                                                                                                                                                                                        caacaaAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCAGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAAcgcg                                                                          >  1:291959/1‑186 (MQ=255)                                                                                                                                                                                                                                               aTGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCAcc  >  1:564614/1‑251 (MQ=255)                                                                                                                                                                                                                                                  ccTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAagcagc                                                                                                                         <  1:164950/129‑1 (MQ=255)                                                                                                                                                                                                                                                  ccTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAagcagc                                                                                                                         >  2:164950/1‑129 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACC  >  NC_000913/2142352‑2142840

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