Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,358	C→T	G252G (GGC→GGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,358	0	C	T	98.0%	142.6 / ‑4.0	51	G252G (GGC→GGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (33/16);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

CCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGT  >  NC_000913/4620114‑4620602                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAgg‑t                                                                                                                                                                                                                                                 <  1:78625/250‑1 (MQ=255)   gCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAgg‑tggt                                                                                                                                                                                                                                              >  1:307761/1‑251 (MQ=255)      tatcccGTGTTCCAGTTTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGACCTGTGCGAAACCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGttt                                                                                                                                                                                                                                            <  2:503852/246‑1 (MQ=255)          ccGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTAAAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACAACCTGGCTGTTGAGCCGCCAGCACCGGCCGTGCTGCAGAAACTGGTTGATGGAAA‑CACCGGGCAGG‑GGGTTTCTGTc                                                                                                                                                                                                                                       >  1:52598/1‑251 (MQ=255)                 ccAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑GGGTTTCTGTCGGTaaa                                                                                                                                                                                                                                 >  1:233418/1‑250 (MQ=255)                  cAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCCCGACCTGGCTGTTGGACCCCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGGAAAACAC‑GGTCAGGGTGGTT‑CTTTCGGGAAAAt                                                                                                                                                                                                                               >  2:393229/1‑251 (MQ=255)                      ttgccgtccATGAAGAAACTTTCGGTCTGGATAACCTCTCCGAACTGTGCGAAATCGCCCTTGAAGAGCTGACCAACGGCGGCGACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGc                                                                                                                                                                                                                            <  2:100512/243‑1 (MQ=255)                               acgaaggatccttgggccTGGACAAACTCTACGAATTGTGCGAAATCGCCCGTGACGAGCGGACCACAGGCGGCTTTAATATCGGTGGTGTTATCGCTCGTCCCTTTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACTGCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTAc                                                                                                                                                                                                                  <  2:377518/234‑1 (MQ=255)                                                    aTAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGACGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTTACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGCGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAAG‑GGGTTTATGTCGGTAAAATTGCAGGCATCTAAGCCCACAGCGGTAttcaca                                                                                                                                                                                               >  2:305268/1‑244 (MQ=255)                                                     tAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCa                                                                                                                                                                                               >  1:420452/1‑248 (MQ=255)                                                       aaCTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGGAAAATTGCAGACCTCTACGCCAACTGCGGTATCAcaaagaa                                                                                                                                                                                           >  1:163271/1‑244 (MQ=255)                                                        aCTCTACGAACTGTGCGAAATCGCCCGTGACGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAg                                                                                                                                                                                         <  1:59256/251‑1 (MQ=255)                                                          tctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCCACTGCGGTATCACCAAGAAAGTg                                                                                                                                                                                       >  1:111492/1‑251 (MQ=255)                                                            tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTCCAACCTCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGTCCGTGCTGCCGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGa                                                                                                                                                                                      >  1:189518/1‑250 (MQ=255)                                                              cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGGTGATGAAAAAAAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTTAGCCAACTGCGGTATCACCAAGAAAGTGaaa                                                                                                                                                                                    >  1:409257/1‑250 (MQ=255)                                                              cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCCGCACCGACCGGGCTGCAGAAACTGGGTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGGAAAATTGCCGACATCTACCCCCACTGCGGGATCACCACGAAAGTGaaa                                                                                                                                                                                    >  2:113761/1‑250 (MQ=255)                                                              cGAACTGTGCGAAATCGCCCGTGAAGAGCGGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAg                                                                                                                                                                                   <  2:189518/251‑1 (MQ=255)                                                                accTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCCCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGca                                                                                                                                                                                 <  2:17005/249‑2 (MQ=255)                                                                 aCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACCACATCGGTCGTTTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAGAAAAAC‑GGTCAGG‑TGGTTTCTGTCGGGAAAATTGCAGACATCATCGCCAACTGCGGTATCACCacgaaaaggaaagca                                                                                                                                                                                 >  2:258683/1‑236 (MQ=255)                                                                              gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCCCGTCCGTTTATCGGTGACACAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAGAAAGGTAAAGCAACTGGcacggaagc                                                                                                                                                                   >  2:155146/1‑243 (MQ=255)                                                                                    gAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTACGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGtt                                                                                                                                                             >  1:447211/1‑251 (MQ=255)                                                                                       gagcggaacaagggaggcTAAAAAATCGGTCGGGTTATCGCGAGGCAGTGTAACGGTGAAAAAGACGGTAACGTCCAGCGTACCGGTAACCGGAACGACCGGGCGGTTGAGCCGCCAGCACCGACCGGGCTGCAGAAACTGGGTGAGGAAAAACAC‑GGTCAGG‑GGGTTTCTGTCGGTAAAAATGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGa                                                                                                                                                          <  2:163428/236‑1 (MQ=255)                                                                                           tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                                                                                                                                      >  1:86622/1‑251 (MQ=255)                                                                                                aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGTTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                                                                                                                                 <  1:46538/251‑1 (MQ=255)                                                                                                                             cGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGCGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAac                                                                                                                     >  1:25325/1‑250 (MQ=255)                                                                                                                                    ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                             >  1:33991/1‑251 (MQ=255)                                                                                                                                    ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                             >  1:61734/1‑251 (MQ=255)                                                                                                                                               tgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGACTTCACCAACt                                                                                                  >  2:98955/2‑251 (MQ=255)                                                                                                                                                        cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACt                                                                                         >  1:241862/1‑251 (MQ=255)                                                                                                                                                            aaCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCg                                                                                      <  1:260362/250‑1 (MQ=255)                                                                                                                                                            aaCTTCCAGCGTACCGGTAACCGGCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCg                                                                                      <  2:233418/250‑1 (MQ=255)                                                                                                                                                             aCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAACGAGATGAAAGAAGCGGGAGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGAc                                                                                    >  2:192085/1‑251 (MQ=255)                                                                                                                                                                 ccAGCGAACCGGTAACCGTCACGACATGGCTGTTGTGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGGAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGCAAGTGAAAGCAGCTGGCCGGGACGCGCTGTTTGAGGCCACCATCAAAGAGATGAACGAAGAGGGAGATAACACCATCGGCTTCACCAACTTCGTTGACTTCGActttt                                                                                >  2:183883/1‑248 (MQ=255)                                                                                                                                                                           ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGTCGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCCACTTCGTTGACGTCGACTCTTCCTgggccc                                                                       >  1:268854/1‑247 (MQ=255)                                                                                                                                                                           ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACAATCGTCTTCACCAACTTCGTTGACTTCTACTCTTCCTGGGGcc                                                                       >  1:264961/1‑250 (MQ=255)                                                                                                                                                                                               tGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGtt                                                  >  2:134740/1‑251 (MQ=255)                                                                                                                                                                                                           aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                           >  1:213141/1‑134 (MQ=255)                                                                                                                                                                                                           aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                           <  2:213141/134‑1 (MQ=255)                                                                                                                                                                                                            gCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                     >  1:130684/1‑251 (MQ=255)                                                                                                                                                                                                              accgaccgTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAATCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCAACGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                   <  1:415435/251‑1 (MQ=255)                                                                                                                                                                                                                    cgTGCTGCAGAAACTGGTTGTTGACACACAA‑GGTCACG‑TGGTTTATGTCGGTAAAATTGCCGGCATCTACGCAAACTGCGGTATCACAAGGAAAGCGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAGGCGGGTGCTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTCTGCTGCGGGTATGGAGCTGt                             >  2:301998/1‑251 (MQ=255)                                                                                                                                                                                                                       gctgaAGAAAATGGTTGAAGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATGGCAGACATCGACGCCAACTGCGGTAACACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTGCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCg                          <  2:373389/251‑1 (MQ=255)                                                                                                                                                                                                                         tgcAGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                        >  1:320015/1‑251 (MQ=255)                                                                                                                                                                                                                            aGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                     >  1:370290/1‑251 (MQ=255)                                                                                                                                                                                                                            aGAAACTGGTTGATGAAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                     >  1:492654/1‑251 (MQ=255)                                                                                                                                                                                                                                          gaAAAACTC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGGGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCt       <  1:459015/251‑1 (MQ=255)                                                                                                                                                                                                                                          gaAAAACAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGc        <  2:307761/250‑1 (MQ=255)                                                                                                                                                                                                                                            aaaaCAC‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtt      >  2:457614/1‑249 (MQ=255)                                                                                                                                                                                                                                               acac‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtgttgt  >  2:157239/1‑251 (MQ=255)                                                                                                                                                                                                                                                cac‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGActc                                                                                  >  2:14200/1‑170 (MQ=255)                                                                                                                                                                                                                                                cac‑GGTCAGG‑TGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGActc                                                                                  <  1:14200/170‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACAC‑GGCCAGG‑TGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGT  >  NC_000913/4620114‑4620602

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