Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,096,754	G→A	G25D (GGT→GAT)	rdgB →	dITP/XTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,096,754	0	G	A	100.0%	157.2 / NA	51	G25D (GGT→GAT)	rdgB	dITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/24);  total (27/24)

CGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGA  >  NC_000913/3096514‑3096999                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      cGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATc                                                                                                                                                                                                                                             <  1:392007/251‑1 (MQ=255)                 aagcTGGTGAATTTTCGGGGTAAGCAATGCCGGGTTGCAAAAAGCCAGGTGGTAATTAAAAAAGGCGAATTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGAc                                                                                                                                                                                                                              <  2:452208/247‑1 (MQ=255)                             tttCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGAtt                                                                                                                                                                                                                 >  2:61593/1‑251 (MQ=255)                                  cGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCg                                                                                                                                                                                                              <  1:61593/249‑1 (MQ=255)                                   ggTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTg                                                                                                                                                                                                           >  1:67015/1‑251 (MQ=255)                                         cAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGaaa                                                                                                                                                                                                     <  1:493847/251‑1 (MQ=255)                                              ccGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGc                                                                                                                                                                                                >  1:242585/1‑251 (MQ=255)                                                                      gatggAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGAtt                                                                                                                                                                        <  1:262373/247‑1 (MQ=255)                                                                       aTTGAAAAAGGCGAACTTTGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTc                                                                                                                                                                       >  1:403562/1‑251 (MQ=255)                                                                         gaaaaaaGGCAAATGTGGCCGCCACAAACAAATTAAAAACATTAAGCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGGATCCTATGCAAAAAGGTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCGGCTGAGCGACTTCGATCTTGATATCGGGGCCCAAACAGACCTCGCCTTTGTTTCCGCTGAAGAAACCGGCCTGACCTTTATCGGAAAGGCGATTCTg                                                                                                                                                                     <  2:602786/249‑1 (MQ=255)                                                                          aaaaaaGGCTAACTAGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAAGCGCGGCGTTAATTAATTAGGTATCCTAGGCAAAAAGTTGTCCTCGCAACCGGCAATGGCGGGAAAGTGCGTGAGCTGGCGTCGCTGCGTAGCGACTTCGATCTTGATATCGTGGCCCAAACATACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTAGCTAAAACGCTATTCTGa                                                                                                                                                                    <  2:370331/250‑1 (MQ=255)                                                                           aaaaaGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGa                                                                                                                                                                    >  1:393121/1‑250 (MQ=255)                                                                                         ggCCGCCACAAACAAAGTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCTGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCgg                                                                                                                                                     <  1:383618/251‑1 (MQ=255)                                                                                            cgccACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGTCCTGAcc                                                                                                                                                                                          <  2:399081/211‑1 (MQ=255)                                                                                            cgccACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGAcc                                                                                                                                                                                          >  1:399081/1‑211 (MQ=255)                                                                                              ccACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCaaaa                                                                                                                                                <  1:281199/251‑1 (MQ=255)                                                                                                    acaaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGAcc                                                                                                                                          >  1:123844/1‑251 (MQ=255)                                                                                                     caaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCg                                                                                                                                         <  1:162672/251‑1 (MQ=255)                                                                                                       aaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCt                                                                                                                                       >  1:559722/1‑251 (MQ=255)                                                                                                        aTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTCATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCtt                                                                                                                                      >  1:541542/1‑251 (MQ=255)                                                                                                        aTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCATCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCtt                                                                                                                                      >  1:39119/1‑251 (MQ=255)                                                                                                                taTTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGACAAGGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGc                                                                                                                               <  1:308654/249‑1 (MQ=255)                                                                                                                    aaTCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAAtt                                                                                                                           <  1:269604/250‑1 (MQ=255)                                                                                                                       ccGCAATAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGcc                                                                                                                        <  1:556096/250‑1 (MQ=255)                                                                                                                           aacaaATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGAAATTGCcgacg                                                                                                                    >  1:62644/1‑250 (MQ=255)                                                                                                                               aaTCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCACAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGCAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGActct                                                                                                               >  2:261322/1‑251 (MQ=255)                                                                                                                                         gAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGTCCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCgg                                                                                                     <  1:123219/251‑1 (MQ=255)                                                                                                                                               gcggcgTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGCAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGCAGAtg                                                                                               >  1:19606/1‑251 (MQ=255)                                                                                                                                                        attaattaGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCGTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTggcg                                                                                      <  2:263024/251‑1 (MQ=255)                                                                                                                                                                 gTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGActc                                                                                                                >  2:88298/1‑216 (MQ=255)                                                                                                                                                                 gTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGActc                                                                                                                <  1:88298/216‑1 (MQ=255)                                                                                                                                                                        aTGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGGCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTTGCGGCGCGCCGGGCAtt                                                                        <  1:446951/249‑1 (MQ=255)                                                                                                                                                                             aaaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcg                                                                 >  1:391264/1‑251 (MQ=255)                                                                                                                                                                                aGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCAGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGCCCCAAACAGACCTCGGCGTTGCATCCTCTGAAGAAACCGGCCTAACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGAGGCAAAAATGACAGCTTTCCCGGCAATTGCCGACCACTCTGGCCTGTCGGAAGATGTGCTTGGCTGCCCGCCGGGCATTTACTCcacgc                                                              >  2:611885/1‑251 (MQ=255)                                                                                                                                                                                   tGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGGTGTGCTTGGCGGCGCGCCGGGCCTTTACTCCGCGCGtt                                                           >  2:185538/1‑251 (MQ=255)                                                                                                                                                                                    gTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTCAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTTCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGtt                                                           >  1:162056/1‑250 (MQ=255)                                                                                                                                                                                      ccTCGCAACCGGCAATGTTGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCGCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTAGCGAAAACGCGAGTCTGAAAGCGCGCCATGCGGCAAAAGTGGCCGCGTTACCGGCAATTGCCGACGACTCTGGTCTGGCGTTAGATGGGCTTGGCGGCGCTCCGGGCATTTACTCCGCGCGTTAtt                                                        <  2:610361/251‑1 (MQ=255)                                                                                                                                                                                         cGCAACCGGCAATGTCGGGAAAGTGCGGGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCGTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGGCTGGCGGTAGATGTGCTTGGCGGCGCGCCGTGCATTTACTCCGCGCGTTATTCCg                                                     <  2:79426/251‑1 (MQ=255)                                                                                                                                                                                          gCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGACCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGCCCTTTCTCGAAAACGCGCTTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGCAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGGTCTTCCgg                                                    >  2:71573/1‑251 (MQ=255)                                                                                                                                                                                            aaCCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGGAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTg                                                  >  2:273929/1‑251 (MQ=255)                                                                                                                                                                                             aCCGGCAATGGCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATATTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTAGCGAAAACTCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGa                                                 <  1:400341/251‑1 (MQ=255)                                                                                                                                                                                                  cAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGAGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcg                                            >  2:437567/1‑251 (MQ=255)                                                                                                                                                                                                    aTGTCGGTAAAGTGCGTGGGCTGGCTCCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCTGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCAATTCTGAAACCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGCACGTCGATGTGCTCGGCGGCGCTCCGGGCATTTACTCCGCGCGTTATTCCgctcaataggct                                          >  2:11687/1‑240 (MQ=255)                                                                                                                                                                                                          gTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGCCCGAt                                    >  2:172083/1‑251 (MQ=255)                                                                                                                                                                                                                        cTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGCAAACGCGATTCTGAAAGCGCGCCATGCGGCAAACGTGACCGCTTTACCGGCAATTGCCGCCGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTTTTCCGGTGAAGACGCGACCGATCAAAAGAATCTGc                       >  2:550361/1‑250 (MQ=255)                                                                                                                                                                                                                          ggCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaa                     >  1:303317/1‑250 (MQ=255)                                                                                                                                                                                                                            cGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaaaa                  <  1:185538/251‑1 (MQ=255)                                                                                                                                                                                                                                    cgtAGCGACTTCGATCTTGATATTGTGGCCCAAACAGACCGCGGAGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGAGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCGGGCGGTAGATGTGCTTGGCGGAGCGCCGGGCATTTACTCCGCGTGTTATTCCGGTGAAGACGCGAACGATCAAAAGAATCTGCAAAAACTGCTGGa          <  1:611885/249‑1 (MQ=255)                                                                                                                                                                                                                                     tcaGCGACTTCGATCTTGATATCGTGGCCCAAACAGCCCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGGAAACGCGATTCTGAAAGCGCGCCATGCGGCAAACGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCGTGGCGGCGCGCCGGGCCTTTACTCCGCTCGTCATTCCGGTGAAGACGCGACCGACCACACGAATCTACAAAACctgcctgc          >  2:377306/3‑246 (MQ=255)                                                                                                                                                                                                                                       cgCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGGCTCTGGTCTGGCTGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGaaa        <  1:389215/249‑1 (MQ=255)                                                                                                                                                                                                                                             ttCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGa  >  1:619241/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      CGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGA  >  NC_000913/3096514‑3096999

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