Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,385	G→A	A261A (GCG→GCA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,385	0	G	A	96.2%	165.3 / NA	53	A261A (GCG→GCA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/16);  minor base C (2/0);  total (37/16)
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.87e-01

ATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTG  >  NC_000913/4620144‑4620628                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     acgaaggatccttgggccTGGACAAACTCTACGAATTGTGCGAAATCGCCCGTGACGAGCGGACCACAGGCGGCTTTAATATCGGTGGTGTTATCGCTCGTCCCTTTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACTGCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTAc                                                                                                                                                                                                                                             <  2:377518/234‑1 (MQ=255)                      aTAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGACGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTTACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGCGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAAGGGGTTTATGTCGGTAAAATTGCAGGCATCTAAGCCCACAGCGGTAttcaca                                                                                                                                                                                                                          >  2:305268/1‑244 (MQ=255)                       tAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCa                                                                                                                                                                                                                          >  1:420452/1‑248 (MQ=255)                         aaCTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGGAAAATTGCAGACCTCTACGCCAACTGCGGTATCAcaaagaa                                                                                                                                                                                                                      >  1:163271/1‑244 (MQ=255)                          aCTCTACGAACTGTGCGAAATCGCCCGTGACGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAg                                                                                                                                                                                                                    <  1:59256/251‑1 (MQ=255)                            tctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCCACTGCGGTATCACCAAGAAAGTg                                                                                                                                                                                                                  >  1:111492/1‑251 (MQ=255)                              tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTCCAACCTCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGTCCGTGCTGCCGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGa                                                                                                                                                                                                                 >  1:189518/1‑250 (MQ=255)                                cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGGTGATGAAAAAAACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTTAGCCAACTGCGGTATCACCAAGAAAGTGaaa                                                                                                                                                                                                               >  1:409257/1‑250 (MQ=255)                                cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCCGCACCGACCGGGCTGCAGAAACTGGGTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGGAAAATTGCCGACATCTACCCCCACTGCGGGATCACCACGAAAGTGaaa                                                                                                                                                                                                               >  2:113761/1‑250 (MQ=255)                                cGAACTGTGCGAAATCGCCCGTGAAGAGCGGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAg                                                                                                                                                                                                              <  2:189518/251‑1 (MQ=255)                                  accTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCCCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGca                                                                                                                                                                                                            <  2:17005/249‑2 (MQ=255)                                   aCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACCACATCGGTCGTTTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAGAAAAACGGTCAGGTGGTTTCTGTCGGGAAAATTGCAGACATCATCGCCAACTGCGGTATCACCacgaaaaggaaagca                                                                                                                                                                                                            >  2:258683/1‑236 (MQ=255)                                                gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCCCGTCCGTTTATCGGTGACACAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAGAAAGGTAAAGCAACTGGcacggaagc                                                                                                                                                                                              >  2:155146/1‑243 (MQ=255)                                                      gAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTACGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGtt                                                                                                                                                                                        >  1:447211/1‑251 (MQ=255)                                                         gagcggaacaagggaggcTAAAAAATCGGTCGGGTTATCGCGAGGCAGTGTAACGGTGAAAAAGACGGTAACGTCCAGCGTACCGGTAACCGGAACGACCGGGCGGTTGAGCCGCCAGCACCGACCGGGCTGCAGAAACTGGGTGAGGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAAATGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGa                                                                                                                                                                                     <  2:163428/236‑1 (MQ=255)                                                             tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                                                                                                                                                                 >  1:86622/1‑251 (MQ=255)                                                                  aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGTTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                                                                                                                                                            <  1:46538/251‑1 (MQ=255)                                                                                               cGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGCGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAac                                                                                                                                                >  1:25325/1‑250 (MQ=255)                                                                                                      ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                                                        >  1:61734/1‑251 (MQ=255)                                                                                                      ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                                                        >  1:33991/1‑251 (MQ=255)                                                                                                                 tgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGACTTCACCA‑ACt                                                                                                                            >  2:98955/2‑251 (MQ=255)                                                                                                                          cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACt                                                                                                                   >  1:241862/1‑251 (MQ=255)                                                                                                                              aaCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCg                                                                                                                <  1:260362/250‑1 (MQ=255)                                                                                                                              aaCTTCCAGCGTACCGGTAACCGGCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCg                                                                                                                <  2:233418/250‑1 (MQ=255)                                                                                                                               aCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAACGAGATGAAAGAAGCGGGAGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGAc                                                                                                              >  2:192085/1‑251 (MQ=255)                                                                                                                                   ccAGCGAACCGGTAACCGTCACGACATGGCTGTTGTGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGGAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGCAAGTGAAAGCAGCTGGCCGGGACGCGCTGTTTGAGGCCACCATCAAAGAGATGAACGAAGAGGGAGATAACACCATCGGCTTCACCA‑ACTTCGTTGACTTCGActttt                                                                                                          >  2:183883/1‑248 (MQ=255)                                                                                                                                             ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGTCGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCC‑ACTTCGTTGACGTCGACTCTTCCTgggccc                                                                                                 >  1:268854/1‑247 (MQ=255)                                                                                                                                             ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACAATCGTCTTCACCA‑ACTTCGTTGACTTCTACTCTTCCTGGGGcc                                                                                                 >  1:264961/1‑250 (MQ=255)                                                                                                                                                                 tGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGtt                                                                            >  2:134740/1‑251 (MQ=255)                                                                                                                                                                             aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                                                      >  1:213141/1‑134 (MQ=255)                                                                                                                                                                             aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                                                      <  2:213141/134‑1 (MQ=255)                                                                                                                                                                              gCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                                               >  1:130684/1‑251 (MQ=255)                                                                                                                                                                                accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAATCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCAACGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                                             <  1:415435/251‑1 (MQ=255)                                                                                                                                                                                      cgTGCTGCAGAAACTGGTTGTTGACACACAAGGTCACGTGGTTTATGTCGGTAAAATTGCCGGCATCTACGCAAACTGCGGTATCACAAGGAAAGCGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAGGCGGGTGCTAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTCTGCTGCGGGTATGGAGCTGt                                                       >  2:301998/1‑251 (MQ=255)                                                                                                                                                                                         gctgaAGAAAATGGTTGAAGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATGGCAGACATCGACGCCAACTGCGGTAACACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTGCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCg                                                    <  2:373389/251‑1 (MQ=255)                                                                                                                                                                                           tgcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                  >  1:320015/1‑251 (MQ=255)                                                                                                                                                                                              aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                               >  1:370290/1‑251 (MQ=255)                                                                                                                                                                                              aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                               >  1:492654/1‑251 (MQ=255)                                                                                                                                                                                                            gaAAAACTCGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGGGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCt                                 <  1:459015/251‑1 (MQ=255)                                                                                                                                                                                                            gaAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGc                                  <  2:307761/250‑1 (MQ=255)                                                                                                                                                                                                              aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtt                                >  2:457614/1‑249 (MQ=255)                                                                                                                                                                                                                aacacgttcAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATg                             <  1:432244/243‑1 (MQ=255)                                                                                                                                                                                                                 acacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtgttgt                            >  2:157239/1‑251 (MQ=255)                                                                                                                                                                                                                  cacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGActc                                                                                                            >  2:14200/1‑170 (MQ=255)                                                                                                                                                                                                                  cacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGActc                                                                                                            <  1:14200/170‑1 (MQ=255)                                                                                                                                                                                                                          ggtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCtgctgc                   >  1:351695/1‑251 (MQ=255)                                                                                                                                                                                                                          ggtggatcatgacGGTAAAATTGCAGACATCTCCGACAACTGCGGTATAACCAAGACAGACAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTC‑CCACACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTTTGGGGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCtgctgc                   >  2:493516/13‑251 (MQ=255)                                                                                                                                                                                                                                   tGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGGGCTGTTCGACCGCCTTCTGCCGTGGCTGATGTCTCTTCTGCGCGagtacg          >  2:348502/1‑246 (MQ=255)                                                                                                                                                                                                                                      cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                     >  1:361863/1‑237 (MQ=255)                                                                                                                                                                                                                                      cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                     <  2:361863/237‑1 (MQ=255)                                                                                                                                                                                                                                       ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt      >  2:406045/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt      >  1:319837/1‑251 (MQ=255)                                                                                                                                                                                                                                           aaaTTGCAGACATCTACGCCAACTGCGGTATCACCACGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTAGGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTg  >  2:73568/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     ATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTG  >  NC_000913/4620144‑4620628

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