Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	158,615	C→A	S171I (AGC→ATC)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	158,615	0	C	A	97.0%	208.9 / NA	68	S171I (AGC→ATC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (33/32);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTTAAGCGATTCTTCAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTC  >  NC_000913/158369‑158859                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      aGTTTAAGCGATTCTGCAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGGGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTa                                                                                                                                                                                                                                                   <  2:136899/251‑1 (MQ=255) aGTTGAAGCGATTCTGCAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACTGACGGTAAAATCCGCTACGATGTa                                                                                                                                                                                                                                                   <  2:101010/251‑1 (MQ=255)    ttAAGCGATTCTTCAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAAt                                                                                                                                                                                                                                                 >  1:352403/1‑250 (MQ=255)      aaGCGATTCTTCAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATAc                                                                                                                                                                                                                                               >  1:30944/1‑250 (MQ=255)               ttcAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTACGGGCTGCTGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTgt                                                                                                                                                                                                                                     >  1:567387/1‑250 (MQ=255)                 cAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGAt                                                                                                                                                                                                                                    <  1:326901/250‑1 (MQ=255)                    aCAGGCGTGCCGGCAGGATATCGTTCAGCAGGGCAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGACGGCAAAACGTACCGCGCGCAGCATACGTACCGGCTCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATACCGACGTCCTTCAGATCCTTCATGCCGCCAACGTCATCACGCACGGTAAAATCCGCTACGCTGAAATACAGGCTGATGGTAGTg                                                                                                                                                                                                                                >  2:116402/1‑251 (MQ=255)                      gggCGTGCCGGTGGGATATCGGTCAGCAGGTTAGCGAGGCGAGGGATCGGTTCTGCGGTTGCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGTCGTCCTTCAGATCCTGCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCAACGATGTAATACAGGCTGTTGATAGTGaa                                                                                                                                                                                                                              <  2:270114/250‑1 (MQ=255)                            gCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATcgcg                                                                                                                                                                                                                        <  1:333475/251‑1 (MQ=255)                                     atatCGGTCAGCAGGGTAGCGAGGCGAGGGATCGGGTCTGCGGTTGCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCGTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAAGCCGTGACGAGGTAATACAGTCTGTTGATAGTGAAATCTCGGCGCTGG‑Gc                                                                                                                                                                                                              <  2:29589/251‑1 (MQ=255)                                              agcagGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCttcttc                                                                                                                                                                                                     <  2:175867/251‑1 (MQ=255)                                                  gggTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATGTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGAGCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGAGCCTTCATGCCGCGAACGTAATCTCGGACGGTAAAATCCGCTACGATGTAATACAGGCTTTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATg                                                                                                                                                                                                 <  2:392749/251‑1 (MQ=255)                                                      aGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGc                                                                                                                                                                                             >  2:207811/1‑251 (MQ=255)                                                         gaggcgaggGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCGTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGGTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGa                                                                                                                                                                                          <  1:242231/251‑1 (MQ=255)                                                               aggGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCTCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAAtgt                                                                                                                                                                                    >  1:337587/1‑251 (MQ=255)                                                                     cGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCCAAACGTACCGCGCGCAGCATCCGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTTATAGTGAAATCGCTGCGCTGG‑GCGTCTTCTTCGATGGGGACGAAAATGTTGTcgc                                                                                                                                                                              >  2:318711/1‑251 (MQ=255)                                                                      ggTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgcg                                                                                                                                                                             <  1:440917/251‑1 (MQ=255)                                                                           tGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc                                                                                                                                                                         >  1:312888/1‑250 (MQ=255)                                                                                tttCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATg                                                                                                                                                                   >  1:108059/1‑251 (MQ=255)                                                                                 ttCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATg                                                                                                                                                                   <  1:80607/250‑1 (MQ=255)                                                                                  tCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGcc                                                                                                                                                                 >  1:15439/1‑251 (MQ=255)                                                                                   ccGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTTCCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACTGACGGTAAAATCCGCTACGATGTAATACAGGCTGTGGATAGTGAAATCGCGGCGCTGG‑GCGTCTGCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCg                                                                                                                                                                <  2:93510/251‑1 (MQ=255)                                                                                        cTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTg                                                                                                                                                           >  1:33478/1‑251 (MQ=255)                                                                                          gATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGcc                                                                                                                                                         >  2:10087/1‑251 (MQ=255)                                                                                              cgcATACCCAATTTTGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTGTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCcgcg                                                                                                                                                     <  1:440040/251‑1 (MQ=255)                                                                                                   accaaaTTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGAGAACGGCTGCCTGCAGATCCTTCATGCCGCCAGCGTAAACACGGACGGTAAAATCCGCTACGAGGTAATACATGCTGTTGAGAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGGCGAATATGTTGTGGCGCAGTAACATGCCGGTTTGCGCGCGTTgg                                                                                                                                                 <  2:366400/246‑1 (MQ=255)                                                                                                     ccAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAg                                                                                                                                              <  1:326646/251‑1 (MQ=255)                                                                                                       aaTTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGt                                                                                                                                            >  1:414695/1‑251 (MQ=255)                                                                                                                 cAAAACGTACAGCGCGCAGCATACGTACCGGATCTTCACGGTAGAGCGTTTCCGGGTTACACCTCAGACGGATACCGCCGTCCTTCAGCTCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTACGATGGAGCCGAAACTGTTTTCGCGCAGCACCATGCCGTTTTGCGCGCGTCgggcgtgcttgcggtcg                                                                                                                                  >  2:282349/1‑237 (MQ=255)                                                                                                                    aaCGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTg                                                                                                                               >  1:297943/1‑251 (MQ=255)                                                                                                                      cGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGAc                                                                                                                             <  1:56784/251‑1 (MQ=255)                                                                                                                        tACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTTG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGt                                                                                                                           >  2:325685/1‑251 (MQ=255)                                                                                                                             cgcgcAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGGTGGGAGGTCGTGCGGTCGCTGACGTTACCt                                                                                                                      >  2:511830/1‑251 (MQ=255)                                                                                                                                gcagcaTACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCg                                                                                                                   >  2:168148/1‑251 (MQ=255)                                                                                                                                  agcaTACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGCAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtg                                                                                                                 >  1:535258/1‑251 (MQ=255)                                                                                                                                           aCCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCCATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACTGACGGTAAAAGCCGCTATGGTGTAATACAGGCTGTTGATAGTGAAATCGCGTCGCTGG‑GCGTCTTCTTCGAGGGAGCCGAATAGGGTGTCGCGCAGCAACATGCCGGTTTGGCCGCGTTGGGATGTCGTGCGTTCGCTGACGTTACTGTCGTGGTGTCCAc                                                                                                         <  2:87589/250‑1 (MQ=255)                                                                                                                                              ggATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                                                                                      >  2:427257/1‑250 (MQ=255)                                                                                                                                               gATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                                                                                      >  1:457936/1‑249 (MQ=255)                                                                                                                                                  cTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGAAAAATCCGCTACGATGGAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCCGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTc                                                                                                 >  2:198935/1‑251 (MQ=255)                                                                                                                                                    tCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCGCGGAAGGTc                                                                                                 <  1:470291/249‑1 (MQ=255)                                                                                                                                                    tCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGATACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTc                                                                                                 <  1:212760/249‑1 (MQ=255)                                                                                                                                                     cACGGTAGCGCCTTTCCGGGTTACCAATCAGACGGATAACGCCGTACTTCAGATCCTTAATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTACTACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCCGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGCCGTTACCTTCGTGGTGTCCACGGAAGGTCGcc                                                                                              >  2:311920/1‑250 (MQ=255)                                                                                                                                                      aCGGTAGCGCGTTTCCGGGTTACCAATGAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTTAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAAGGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGTTCGCTGACGTTACCTTCGTGGTTTCCACGGAAGGTCGCaa                                                                                             <  2:558808/251‑1 (MQ=255)                                                                                                                                                      aCGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCa                                                                                              >  2:206637/1‑250 (MQ=255)                                                                                                                                                                gTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGGTGGGAGGTCGTGCGGTCGCTGCCGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAt                                                                                   >  2:149841/1‑251 (MQ=255)                                                                                                                                                                   tCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTc                                                                                                                                           >  1:289500/1‑192 (MQ=255)                                                                                                                                                                   tCCGGGTTACCAAGCAGACGGATAACGCCGGCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAAGCGCGGCGCTGG‑GCGTCGTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTGGCCCGCGTTGGGAGGTc                                                                                                                                           <  2:289500/192‑1 (MQ=255)                                                                                                                                                                    ccGGGTTACCAATCAGACGGATAACGCCGTCCTGCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTGCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGATGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtatc                                                                                <  2:187072/250‑4 (MQ=255)                                                                                                                                                                        gTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGCTAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGc                                                                           >  1:193546/1‑251 (MQ=255)                                                                                                                                                                          tACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCGCGGACGGTAAAAGCCTCTACGATGTAATGCAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTGTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAATGTCGCAACTTCGATAATCTCCGGGcc                                                                          <  2:173010/250‑1 (MQ=255)                                                                                                                                                                           aCCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTGCTTCGATGGAGCCGAAAATGTTGTCGCGTAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGGTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCa                                                                         <  1:263367/250‑1 (MQ=255)                                                                                                                                                                                    aCGGATAAAGCCGTCCTTCAGATCCTTCATGCCTCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAAGACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACAt                                                               <  2:213874/251‑1 (MQ=255)                                                                                                                                                                                      ggATAACGCTGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTGCGATGTAAGACAGGCTGTTGATAGTGAAATCGAGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGa                                                             <  2:435284/251‑1 (MQ=255)                                                                                                                                                                                      ggATAACGCCGTCCTTCAGATCCTTCATGCCGCAAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGCTAGTGAAATCGCGGCGCTGG‑GCGACTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGCCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCCAACATTACATgg                                                             >  1:391365/1‑250 (MQ=255)                                                                                                                                                                                         tAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGcc                                                          >  1:89060/1‑251 (MQ=255)                                                                                                                                                                                          aaCGCCGTCCTTCAGATCCGTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCTGGCCAAACATTACATGAGCCa                                                         <  1:406296/251‑1 (MQ=255)                                                                                                                                                                                                  ccTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGCTGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGaa                                                  >  1:300699/1‑250 (MQ=255)                                                                                                                                                                                                     tCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGTAAAc                                                <  2:352403/249‑1 (MQ=255)                                                                                                                                                                                                                aTGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACccacc                                     >  1:113427/1‑249 (MQ=255)                                                                                                                                                                                                                         aCGTAATCACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGtt                          <  1:168148/251‑1 (MQ=255)                                                                                                                                                                                                                                cACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGaa                    <  1:10087/250‑1 (MQ=255)                                                                                                                                                                                                                                cACGGACGGTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAAc                   >  2:208944/1‑251 (MQ=255)                                                                                                                                                                                                                                 acggacggTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACa                  <  1:80695/251‑1 (MQ=255)                                                                                                                                                                                                                                    gacggTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGtt               >  2:28996/1‑251 (MQ=255)                                                                                                                                                                                                                                      cggTAAAATCCGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTg             <  1:536519/251‑1 (MQ=255)                                                                                                                                                                                                                                        gTAAAATCCGCTACTATGTAATTCAGGCTGTTGAGAGTGAAATCGCGGCG‑TGGTGCGTCTTCGTCGATGGAGCCGAAAATGTTGTCGCGCTGCAACAGGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGGGGTGTTCACGGAAGGTCGCAACGTCGAGAAGCTCCGGGCGAAACATTTCATGATTGAGGCGTAAACGGCGGCCTATCAGGCGTCATTGACGCAACAGTTTgcg           <  2:585824/251‑1 (MQ=255)                                                                                                                                                                                                                                               ccGCTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGCTGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGc    >  1:469610/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cTACGATGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTc  >  1:286095/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      AGTTTAAGCGATTCTTCAAACAGGCGTGCCGGTGGGATATCGTTCAGCAGGGTAGCGAGGCGAGGGATCGGTTCTGCGGTTTCCGGGCTGATGCGCATACCCAATTTGGCGGCAAAACGTACCGCGCGCAGCATACGTACCGGATCTTCACGGTAGCGCGTTTCCGGGTTACCAATCAGACGGATAACGCCGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGG‑GCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTC  >  NC_000913/158369‑158859

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