Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,204,395	T→A	intergenic (+106/‑201)	yehK → / → yehL	uncharacterized protein YehK/putative AAA(+) MoxR family ATPase YehL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,204,395	0	T	A	100.0%	111.5 / NA	38	intergenic (+106/‑201)	yehK/yehL	uncharacterized protein YehK/putative AAA(+) MoxR family ATPase YehL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (18/20);  total (18/20)

CCCGCTAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTG  >  NC_000913/2204151‑2204638                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    cccGCTAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAatcgttc                                                                                                                                                                                                                                               <  1:445904/251‑1 (MQ=255)              gatgCGTTGACTTATATACGAGATGAGCGTATTGATGAGCGGAAATTCCGGGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCGTTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACCCGAt                                                                                                                                                                                                                                  >  2:254579/1‑251 (MQ=255)                   gTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGCTGTTGc                                                                                                                                                                                                                             >  1:390418/1‑251 (MQ=255)                           atataAGAGAGGAGGGTATTGAGGAGCGGTAATTACGAGGTGAATCTCCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTTGTCCGATCCGGCGCATCGTCCCAGTGGCCTGGATTATAAATCTCATGATCTGAATTGCAAAGGGGGTCAGCCGTGGTTATAAAGCGGTCGTTGTATCACGATTTTGCTCAGGCt                                                                                                                                                                                                                      <  2:387879/250‑1 (MQ=255)                                agagATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTgtggtg                                                                                                                                                                                                                <  1:155969/251‑1 (MQ=255)                                  agaTGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGtt                                                                                                                                                                                                              >  1:303200/1‑251 (MQ=255)                                            aTTGATGAGCGTAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATAGATGCAACCAAGAAATAGAATAATATTCTCTAAATAATTCAAATTGGTCCGATCCGGCTCAACGTCCCAATGGCCCTGATTAGAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGTCTGTTGTGTTAGGCCTGaa                                                                                                                                                                                                     <  2:388603/250‑1 (MQ=255)                                                     cGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGtaaat                                                                                                                                                                                           >  2:161830/1‑251 (MQ=255)                                                       gAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGtaaataa                                                                                                                                                                                         <  2:189026/251‑1 (MQ=255)                                                          aTTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAg                                                                                                                                                                                      <  2:187800/251‑1 (MQ=255)                                                                       acACCTGAAACTGGAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTACCAATGGCCTGGATTATCAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAACAGAACCTGCCTCAcc                                                                                                                                                                         >  2:312068/2‑251 (MQ=255)                                                                                       tAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttatt                                                                                                                                                         >  1:343853/1‑251 (MQ=255)                                                                                                             aCAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAg                                                                                                                                     <  1:187888/249‑1 (MQ=255)                                                                                                                        tCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTcc                                                                                                                                       >  1:7427/1‑236 (MQ=255)                                                                                                                        tCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTcc                                                                                                                                       <  2:7427/236‑1 (MQ=255)                                                                                                                               aGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATAGTTTTCGCTCAaacaac                                                                                                                 <  2:372952/251‑1 (MQ=255)                                                                                                                                gAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACt                                                                                                                >  1:568728/1‑251 (MQ=255)                                                                                                                                       aataatATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATgcgc                                                                                                         >  1:143196/1‑251 (MQ=255)                                                                                                                                         atatatACTCTAAATAATTCAAATTGGCCCAATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCAGTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGGAAATAAAAGAACCTGCCTCACCAGCAGGTTTTGTTATTTACTGTGATCTGCTATCCAGATATTGGTCGCTCATACAACTAATGCGCCa                                                                                                       <  2:62241/249‑1 (MQ=255)                                                                                                                                                ctctAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACAttt                                                                                                >  1:31408/1‑251 (MQ=255)                                                                                                                                                  ctAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGTATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACAttt                                                                                                <  1:8700/249‑1 (MQ=255)                                                                                                                                                  ctAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACAttt                                                                                                <  1:8670/249‑1 (MQ=255)                                                                                                                                                  ctAAATAATTCAAATTGGTCCGATCCGGCGCAAAGTCCCAATGGACTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGACGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTAACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGAGCCAAACATTTAt                                                                                              >  1:627066/1‑251 (MQ=255)                                                                                                                                                            cAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTaa                                                                                     >  1:49523/1‑250 (MQ=255)                                                                                                                                                              aaTTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACAGTTATTGCGCGTAAAat                                                                                  <  2:598788/251‑1 (MQ=255)                                                                                                                                                                        aTCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATAGCGTTTAtt                                                                        <  1:161830/251‑1 (MQ=255)                                                                                                                                                                                  aCGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATAc                                                              <  1:384099/251‑1 (MQ=255)                                                                                                                                                                                      accAATGGCCTGGATTATAACTCTCATTCTCTTAATTGCAACGGGGTCCAGCCGAGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATGTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCAAAACATTTATTGCCAGTAAAATATCGTTTATTTCATTAATCCAttt                                                          >  2:145150/2‑251 (MQ=255)                                                                                                                                                                                               cTGGATTATAAATCTCATTATCTTCATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGCTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAACTAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAAAATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATg                                                  >  2:328163/1‑250 (MQ=255)                                                                                                                                                                                                      ataAATCTCAGTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAAGCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTCTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTTAAACAACTAATGCGCCAAACAGTTATTGCGCGTAAAAGATCGTGTATTTCATTAATACATTTCAGGGATGAATATATg                                          <  2:291226/251‑1 (MQ=255)                                                                                                                                                                                                            ctcATTATCTTAATTGCACCGGGGTACAGCCGTGGTTTTAAATCGTTCGTTGGTACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGTATATAGTAATTAAATGAACCTGCCTCACCAGCAGGTTTTTTTTTTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTACTGCGCACAACATTTATTTCGCGTATATTATCGTTTATTTCCTTATTCCATTTCTGGGTTGATTATATGTCACCg                                    >  2:569409/1‑251 (MQ=255)                                                                                                                                                                                                              taTTATCTTAAGTGCAACTGGGTCGAGCCGTGGTTATAAATCGTTCGTTGGAACACGATTGTGCACAGGCTTTGGGGTTAGGCCTTAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCAc                                      <  2:456209/246‑1 (MQ=255)                                                                                                                                                                                                              cATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGc                                   >  1:130573/1‑250 (MQ=255)                                                                                                                                                                                                                       aaTTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAatc                          <  1:233630/250‑1 (MQ=255)                                                                                                                                                                                                                       aaTTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTGTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAatc                          <  2:234250/250‑1 (MQ=255)                                                                                                                                                                                                                                 gggTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGt               <  1:604574/251‑1 (MQ=255)                                                                                                                                                                                                                                    tCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCg            >  2:619895/1‑251 (MQ=255)                                                                                                                                                                                                                                      cAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGcc          >  1:369977/1‑251 (MQ=255)                                                                                                                                                                                                                                              gTTATAAAGCGTTCATTGGAACACGATGGTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGATCCTGCCTCACCAGCAGGTTTTTATATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCTCctgctgctg  <  2:470731/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    CCCGCTAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTG  >  NC_000913/2204151‑2204638

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