Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	229,552	C→T	S129F (TCC→TTC)	dkgB →	methylglyoxal reductase DkgB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	229,552	0	C	T	97.8%	164.4 / NA	45	S129F (TCC→TTC)	dkgB	methylglyoxal reductase DkgB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/14);  minor base G (0/1);  total (30/15)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGCAGAAAGTGGCGTGCCACGTCATGAACTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGT  >  NC_000913/229313‑229794                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   tGCAGAACGTGGCGTGCAACGTCCTGAACTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAACGCCAACTTGATCCCACGACTGAAAGAGAGCCTGCAAAAATTGCGAACCGATTATGTTGATCTGACGCTACTCCACTGGCCGTCACCAAACGATGAAGTATCTGTTGAAGAGTTTATGCAGGAGCTGCTGGAAACCAACAAACAAGGGCTGACGCGAGAGATCGGTATTTTCAACTTCActg                                                                                                                                                                                                                                         >  2:190021/1‑249 (MQ=255)  gCAGAAAGTGGCGTGCCACGTCATGAACTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGAt                                                                                                                                                                                                                                        >  1:287975/1‑251 (MQ=255)                       cATGAACTCTACATCACCACTAAAATCTTGATTGAAAATCGCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGa                                                                                                                                                                                                                    <  2:187645/250‑1 (MQ=255)                          gAACTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATtgc                                                                                                                                                                                                                >  1:65630/1‑251 (MQ=255)                           aaCTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATtgct                                                                                                                                                                                                               >  2:7276/1‑251 (MQ=255)                           aaCTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATtgct                                                                                                                                                                                                               <  1:2045/251‑1 (MQ=255)                            aCTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAACATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATtgctg                                                                                                                                                                                                              >  2:401740/1‑251 (MQ=255)                                     cagcattAAAAGCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTTAAAGAGAGCCGGTAAAAATTGCGTACCGATTATGGTGATCTGTGGCTAATCCACTGTCCGTCACCAAACGATGAATTCTCTGTTTAAGAGGTGATGCAGGCGCTGCTGGAAGGCAAAAAACAAGTTCTGAGGCGTGAGATCGTTATGTGCAGCTTCACGATCCCGTGGATGGAAAAAGCGATTGCTGCTGTTGGTg                                                                                                                                                                                                     <  1:25051/245‑1 (MQ=255)                                              aaTCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAAc                                                                                                                                                                                             >  2:196559/1‑250 (MQ=255)                                                 cTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCg                                                                                                                                                                                         >  1:46441/1‑251 (MQ=255)                                                              ctcAGCACAGACAAATTGATCCCAAGTCTGAAAGAGGGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAAcc                                                                                                                                                                                >  2:147582/1‑247 (MQ=255)                                                              ctcAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAAcc                                                                                                                                                                                >  2:34866/1‑247 (MQ=255)                                                                cAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGAt                                                                                                                                                                            <  1:131132/249‑1 (MQ=255)                                                                     aataCAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAActc                                                                                                                                                                     >  1:259935/4‑251 (MQ=255)                                                                      aGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAActct                                                                                                                                                                    >  1:246923/1‑251 (MQ=255)                                                                       gACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAACACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAActctc                                                                                                                                                                   >  1:21316/1‑251 (MQ=255)                                                                         cAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACtctctc                                                                                                                                                                 >  1:123282/1‑251 (MQ=255)                                                                          aaaTTGATCCCAAGTCGGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAAACCACTGGCCGTCACCAAACGATGAAGTCGCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTcc                                                                                                                                                                <  2:246923/251‑1 (MQ=255)                                                                                                   gCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGATGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGAGTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGtt                                                                                                                                       >  1:154330/1‑251 (MQ=255)                                                                                                          aaaaTTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTggg                                                                                                                                >  1:407120/1‑251 (MQ=255)                                                                                                                gCGTACCGAGTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAAc                                                                                                                          <  2:406199/251‑1 (MQ=255)                                                                                                                  gTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAg                                                                                                                        <  1:7276/251‑1 (MQ=255)                                                                                                                                     gACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACt                                                                                                      <  2:65630/250‑1 (MQ=255)                                                                                                                                          tAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCtat                                                                                                >  2:111085/1‑251 (MQ=255)                                                                                                                                           aaTCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGCTCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATATGCAAAACCGTAACGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCtat                                                                                                >  1:223287/1‑250 (MQ=255)                                                                                                                                              ccACTTGCCGTCACCAAACGAGGAAGTCTCTGTTGAAGAGTTTATTCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGTTGCTGAAAACATCGCTACGAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAGCAGCACGGCATCCATATTACTTCCTATATGa                                                                                            <  1:165385/251‑1 (MQ=255)                                                                                                                                                         cACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTAt                                                                                 <  1:175860/251‑1 (MQ=255)                                                                                                                                                                cGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAgg                                                                          >  1:138341/1‑251 (MQ=255)                                                                                                                                                                    gAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCAGCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                      <  1:374474/251‑1 (MQ=255)                                                                                                                                                                        tctcTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCCTATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaa                                                                   >  2:97610/1‑250 (MQ=255)                                                                                                                                                                           ctGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaa                                                                   >  1:253320/1‑247 (MQ=255)                                                                                                                                                                              ttGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATg                                                              >  2:176216/1‑249 (MQ=255)                                                                                                                                                                                    agTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAtt                                                      <  2:237569/251‑1 (MQ=255)                                                                                                                                                                                     gTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTg                                                     <  1:111085/251‑1 (MQ=255)                                                                                                                                                                                             aGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGTTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATc                                             >  1:28443/1‑251 (MQ=255)                                                                                                                                                                                              ggCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGt                                                                                                                                                >  1:271775/1‑151 (MQ=255)                                                                                                                                                                                              ggCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGt                                                                                                                                                <  2:271775/151‑1 (MQ=255)                                                                                                                                                                                                       ggAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                    <  1:356020/250‑1 (MQ=255)                                                                                                                                                                                                                 aaaCAGGGGCTGACGCGTGAGATAGGTATTTTCAACTTCACGATCCAGTTGATGGAAAACGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGACCTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGCAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGAc                         >  2:71947/1‑251 (MQ=255)                                                                                                                                                                                                                 aaaCAAGTGCTGACGCGGGCGATCGGGATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATATGCAAAAACGTAAAGAGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTCTGGTAAGTCCATGAAAGATGAGGTTATTGCACGTATCGCAGCTAACAACCATGCGAc                         >  2:248987/1‑251 (MQ=255)                                                                                                                                                                                                                   aCAAGGGCTGACGCGTGAGCTCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTc                       >  1:358859/1‑251 (MQ=255)                                                                                                                                                                                                                         gCTGACGCGTGAGATCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGcac                 >  1:379419/1‑251 (MQ=255)                                                                                                                                                                                                                                    agaTCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTg      >  2:393554/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTCGGTATTTTCAACTTCACGATCCCGTTGATGGAAAAAGCGCTTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGc    >  1:128816/1‑251 (MQ=255)                                                                                                                                                                                                                                        cGGTATTTTCAAATTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTACAGTGGATGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTACGGCCCTGAAAGATCAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGt  >  1:58378/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   TGCAGAAAGTGGCGTGCCACGTCATGAACTCTACATCACCACTAAAATCTGGATTGAAAATCTCAGCAAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGGTGCTGAAAACATCGCTACTAACCAGATTGAACTCTCTCCTTATCTGCAAAACCGTAAAGTGGTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGT  >  NC_000913/229313‑229794

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