Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,625,546	C→T	L218L (CTG→TTG)	serB →	phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,625,546	0	C	T	98.4%	185.3 / NA	61	L218L (CTG→TTG)	serB	phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (35/25);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGAAA  >  NC_000913/4625302‑4625786                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gggCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTCATTTCTTTGCTGGATATCTGCGCGCCAAGTTGCGcc                                                                                                                                                                                                                                            >  2:386321/1‑251 (MQ=255)            gggcgaaataaacggacgggggggggcgcGGCAAACTCGATTTTACCGCCAGCCTGCGCCGCCGCGTGGCGACGCTGAAAGGCGCTGACGCAAATATCCTGCAACGGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGCCCGCCGTg                                                                                                                                                                                                                                 <  2:390732/226‑1 (MQ=255)               cGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGt                                                                                                                                                                                                                               >  1:325151/1‑250 (MQ=255)                ggAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGttg                                                                                                                                                                                                                             >  2:400614/1‑249 (MQ=255)                              gggCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTgg                                                                                                                                                                                                                 >  1:366724/1‑249 (MQ=255)                                gCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATAGCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGAt                                                                                                                                                                                                             <  1:391503/251‑1 (MQ=255)                                        cggcgAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACg                                                                                                                                                                                                     >  1:82008/1‑251 (MQ=255)                                           ggaactggtttttCCCGCCAGCCTGCGCAGCGGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTAAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACGTTCGTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCAGGGACGGTa                                                                                                                                                                                                  <  2:247665/242‑1 (MQ=255)                                              atcctaTTTTCCCGCCAGCCGGCGCAGCCGCGGGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAt                                                                                                                                                                                               <  2:133707/246‑1 (MQ=255)                                                  gATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAc                                                                                                                                                                                           >  1:93001/1‑251 (MQ=255)                                                      ttcccGCCAGCCTGCGCAGCCGCGTGGCGCAGCGGAAAGGCGCTGACGCCAATATCCTGCAACAGTTGCTTGAAAACCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGTTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATTTCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCgg                                                                                                                                                                                        <  2:468331/247‑1 (MQ=255)                                                      ttACCGCCAGCCTGAGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTACCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGACAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGATCTGGAGATCATGGACGGTAAATTTACCGGc                                                                                                                                                                                       >  1:19903/1‑251 (MQ=255)                                                       tACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGc                                                                                                                                                                                       >  2:251398/1‑250 (MQ=255)                                                         ccGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAt                                                                                                                                                                                    >  1:386981/1‑251 (MQ=255)                                                               gCCTGCGAGGCCGCGTGGCGATGCTAAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCTGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGATGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATc                                                                                                                                                                              <  1:189948/251‑1 (MQ=255)                                                                ccTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCg                                                                                                                                                                             <  1:251398/251‑1 (MQ=255)                                                                 atGCGCAGCAGCGTGGCGACGCTGAAAGACGCAGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAtg                                                                                                                                                                                   >  1:367517/2‑244 (MQ=255)                                                                                 ggACGCGGAAAGGGCCTGACGCCAATCTTCGGCAACAGTGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGTTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTTCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAcgcg                                                                                                                                                            <  2:346208/250‑1 (MQ=255)                                                                                   aCGCTGAAAGGCGCTGCCGCCAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTAAATATCTGCGCGACAAGTTGAGCCTGACCGCCGTTGGAGCCAATTAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGCTCGGCGGCATCGTAGAcgccc                                                                                                                                                           >  2:240315/1‑248 (MQ=255)                                                                                                       cAATATCCGGCAAAAGGGGCGTGAAAAGCTGCCGATGATGCCAGGCTTAACGCAACTGGGTCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACGTTCTTTGCTGAATATCTGCGCGCCAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGAACGGCGACAGCGTAGACGCGCAGTACAAAGCGAAAActccga                                                                                                                                      <  2:19903/251‑4 (MQ=255)                                                                                                         atatTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGAct                                                                                                                                    <  1:59292/251‑1 (MQ=255)                                                                                                          tatTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGActc                                                                                                                                   >  1:413668/1‑251 (MQ=255)                                                                                                          tatTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGActc                                                                                                                                   >  1:355792/1‑251 (MQ=255)                                                                                                                 caacagtggCGTGAAAATTCGCCGCTGTTGCCAGGCTTACCCCAACTGGTGTCCATGCTGGAAACGCTGGTCTGGAAAGTGGAGATTGCCTCCGGCGGCTTTACTTTTTTTGCTGATTATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACGCTGActcgcctc                                                                                                                              <  2:354239/241‑1 (MQ=255)                                                                                                                 cAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAcc                                                                                                                                                                                          <  1:366753/189‑1 (MQ=255)                                                                                                                 cAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAcc                                                                                                                                                                                          >  2:366753/1‑189 (MQ=255)                                                                                                                 cAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGAAGACGCGCAGTCCAAAGCGAAAACTCTGActcgcctcgc                                                                                                                            >  2:122049/1‑251 (MQ=255)                                                                                                                   aCAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCgcgc                                                                                                                          >  1:220517/1‑251 (MQ=255)                                                                                                                     aGGTGCGTGAAATTCGGCCGCGGATGCCAGGCTTAACGCAACTGGGGCTCAAGCTGGAAACGCTGGGCTGGAAAGGGGCGATTGCCTCCGGCGGCTTTACGTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAg                                                                                                                        <  2:366724/251‑1 (MQ=255)                                                                                                                     aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTa                                                                                                                                                                  >  1:240424/1‑209 (MQ=255)                                                                                                                     aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTa                                                                                                                                                                  <  2:240424/209‑1 (MQ=255)                                                                                                                     aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCgcgc                                                                                                                          >  2:127366/1‑249 (MQ=255)                                                                                                                          cGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGt                                                                                                                    >  1:156500/1‑250 (MQ=255)                                                                                                                           gTGAAAATCTGCTGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGCCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTAt                                                                                                                  <  1:319566/251‑1 (MQ=255)                                                                                                                           gTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTAt                                                                                                                  <  1:439629/251‑1 (MQ=255)                                                                                                                             aaaaaTCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAATTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACTCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATg                                                                                                                 <  1:147676/249‑1 (MQ=255)                                                                                                                                   cTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTCCAAAGCGAAAACTCTGACTCGCCTCGCGCCGGCGTATGAAATccc                                                                                                          >  2:205953/1‑251 (MQ=255)                                                                                                                                   cTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                                          >  1:418433/1‑251 (MQ=255)                                                                                                                                     gCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                                          <  2:141371/249‑1 (MQ=255)                                                                                                                                     gCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATc                                                                                                            <  1:127366/247‑1 (MQ=255)                                                                                                                                      ccGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTACGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTACATTTACCGGAAATGTGATCGGCGACCTCGAAGACGCGCGGTACAAAGCGAAAAATCTGACTCGCCTCGCGCCGGAGTATGAAATCCCGCt                                                                                                       >  2:465212/1‑251 (MQ=255)                                                                                                                                                cAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGGGTCTGAAATCCCGCTGGCGCAGAcc                                                                                             >  2:155909/1‑251 (MQ=255)                                                                                                                                                    cTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtgg                                                                                         >  1:463132/1‑251 (MQ=255)                                                                                                                                                            aaCTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattggc                                                                                 <  1:205226/251‑1 (MQ=255)                                                                                                                                                                 gTGCTCAAGCTAGAAACGCTGGGCTGGAAAGTGTCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCACTGAACTGGAGATCATGGACGGAAAATTTACCGGCACTGAGATCGGCGACATCGTCAACGCGCCGTACAAAGCGGAGACTCTGCCTCGCCTCACGCAGGAGTATGACATCCCGCGGGCGCAGCCCAgggcgagtggcgaacg                                                                            >  2:87470/1‑248 (MQ=255)                                                                                                                                                                         gCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATg                                                                    <  1:400614/251‑1 (MQ=255)                                                                                                                                                                          cTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGTCGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTAAAAAGCGAAAACTCTGACTCGCCCCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGAGATGGCGCCAATg                                                                    >  2:84266/1‑250 (MQ=255)                                                                                                                                                                                                  gcgatgtcctccgtcggtgttattttctgtGCTGAATATCGGCGCGACAAGTTGCGCTTGGCCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGTTAAATTTACCGGCAATGTGATCGGCGACATCTTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAgg                                           <  1:468337/228‑1 (MQ=255)                                                                                                                                                                                                   cGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAGTATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAggg                                          >  1:131912/1‑251 (MQ=255)                                                                                                                                                                                                        gCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAACGCGGCAGGGCTggg                                     >  1:476399/1‑251 (MQ=255)                                                                                                                                                                                                           tCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGAt                                  >  1:29197/1‑251 (MQ=255)                                                                                                                                                                                                             cggcggCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGGTTg                                >  1:286884/1‑251 (MQ=255)                                                                                                                                                                                                                cggcTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCt                             <  1:122049/251‑1 (MQ=255)                                                                                                                                                                                                                      tatttgttttGCTAAATATTTGCGCGAAAAGTTTCGCTTGACCGCCGGGGTAGCCAAGGAACTGGAGATCATGGCCGGTAAATTTACCGGCAATGTGACCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACCCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCAt                        <  2:355792/243‑1 (MQ=255)                                                                                                                                                                                                                         tttctttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgcc                     <  1:155909/250‑1 (MQ=255)                                                                                                                                                                                                                                gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCCTgccaacccaa              >  2:73950/1‑250 (MQ=255)                                                                                                                                                                                                                                gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagccaa              >  1:71427/1‑250 (MQ=255)                                                                                                                                                                                                                                     atatatGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAACGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGAACTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACAATGCCAAGCCCAAAGtg         >  1:20932/6‑250 (MQ=255)                                                                                                                                                                                                                                      tatcTGCGCGACAAGTTGCGCCTGAACGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGAAATGTGATCGGCGACATCGTCGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCACGCGCAGGCGTATGAAATCCCGCTGGCGCAGCACGTGGCGATTGGCGCTGGAGCCAATGACCGGCCGATGATCAAAGCGCCGGGGGTAGGGAGTGCCTACCCTTCCAAGACAAAAGtgc        >  2:340473/4‑249 (MQ=255)                                                                                                                                                                                                                                         cTGCGCGACAAGTTGCGCCTGACCGCCGTGGTATCCAATGCACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGAGAAAACTCTGACTCGCCTCGCGCAGGGGTATGACATCCCGCGGGCGAAGACCGTGGCGCTTGGAGCAGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGACGGGGGTTGCCTaacagtcccagcccaaaggaaatg     >  2:269177/1‑227 (MQ=255)                                                                                                                                                                                                                                             gcgACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGaaa  >  1:6505/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGAAA  >  NC_000913/4625302‑4625786

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