Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,126,316	G→A	E54K (GAA→AAA)	yceH →	DUF480 domain‑containing protein YceH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,126,316	0	G	A	100.0%	117.3 / NA	38	E54K (GAA→AAA)	yceH	DUF480 domain‑containing protein YceH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/17);  total (21/17)

GTTGATTGATGGACAATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGG  >  NC_000913/1126074‑1126556                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 gTTGATTGATGGACAATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGt                                                                                                                                                                                                                                              <  1:424203/247‑1 (MQ=255)       tgatGGACAATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGTTCATGAATCTGAGTGAATCCAAAGTGCAGGAAc                                                                                                                                                                                                                                      >  1:431374/1‑249 (MQ=255)       tgatGGACAATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAAc                                                                                                                                                                                                                                      >  1:528406/1‑249 (MQ=255)           ggACAATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTg                                                                                                                                                                                                                                 <  1:388077/250‑1 (MQ=255)                               aCTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAg                                                                                                                                                                                                                                               >  2:396804/1‑216 (MQ=255)                               aCTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAg                                                                                                                                                                                                                                               <  1:396804/216‑1 (MQ=255)                                                       cTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAACAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGtttt                                                                                                                                                                                    >  1:286275/1‑251 (MQ=255)                                                        tGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAACCGTCATTATCTACGCACAGTGAGCGGTTTTg                                                                                                                                                                                   >  1:435243/1‑251 (MQ=255)                                                         ggACGCCCGGCCGCTAAGGAGAATCTCATGAAATAGCAACTTACTGCCCTTGAAGCTCTCGTGATTGGTTGC‑TTGCTGTAAAATCATGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATAAGAAAACGAACCGTGAACCGGTCATTAATCTGAGTGAATCCAAAGTGCAGGAACAGCTTGATAATCTGGTCAAACGTCATTATCTACGCACAGTTAGCGGTTTTg                                                                                                                                                                                   <  2:625167/250‑1 (MQ=255)                                                                 ggCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATAACCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGATGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCggg                                                                                                                                                                          >  1:664866/1‑251 (MQ=255)                                                                  gCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGCGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGt                                                                                                                                                                         >  2:13349/1‑251 (MQ=255)                                                                                       aaaTATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTCTC‑TTGCTGGAAAAGCAGGTGACCCCGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCATAAAACGATCCGTGATCCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCAGGTCAAACGTCATTATCTACGCCCAGTGAGCGGTTATGGTAATCGGGTCACCAAATATGAGCAACggtt                                                                                                                                                    <  2:57644/251‑4 (MQ=255)                                                                                               aCTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAATCGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAAtt                                                                                                                                            <  1:78897/251‑1 (MQ=255)                                                                                                             aagctatcgTGATAGGATGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTATCAGTCAATGGTGTAGTCACGGCCTGTACTCAGAAAACGAACCGTGAACAGGCCATGACTCTGAGTGAATCCAAAGTGCAGGAACAGGTGGATAATCTGGTAAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGAAATCGGGTCACCCAATATGCGCAACGTTTTTGTAATTAAGAAttggg                                                                                                                                  >  2:243315/8‑244 (MQ=255)                                                                                                             aaGCTCGCTTGAGTGGTTGCCTTAC‑GGAAAAGCTGGTGACCACGCCGGAGAAATACCCGCGCTCAGTCAATGGGGTAGTCCCGGCCTGTAATCACAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGTTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGATCGGTTTTGGTAATCGGGTCACCAAATATGTGCAACGTTTTGGTCATTCATATTTTGGCGAt                                                                                                                              <  2:285524/251‑1 (MQ=255)                                                                                                              aGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATc                                                                                                                             <  1:13349/251‑1 (MQ=255)                                                                                                                       gATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTg                                                                                                                     >  1:493061/1‑250 (MQ=255)                                                                                                                           ggTTGC‑TTGCTAGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAgcgc                                                                                                                >  2:353657/1‑251 (MQ=255)                                                                                                                                                      cccGCCGGAGCAATACTCGCTCTCAGTCAATGGGGTAGGTACGGCCAGTAAGCAGAAAACGAACCGTTAACATGTCAGGAATCGTAGGTAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCGACGCACAGTGATCGGTTTTGGTAATCGGGTCACCAAGTATGAGCAACGTTTTTGTATTTCAGAATTGGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCTCCACgt                                                                                      <  2:322219/249‑1 (MQ=255)                                                                                                                                                        cGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACgttg                                                                                    >  1:235451/1‑251 (MQ=255)                                                                                                                                                           cGGAGCAATACCCGCTCTCAGTCAATGGTGTATTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTa                                                                                 <  2:528406/251‑1 (MQ=255)                                                                                                                                                                   tACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCgtgg                                                                         >  1:412939/1‑251 (MQ=255)                                                                                                                                                                    aCCCTCTCTCAGTACATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCTTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAAGCGGGTCACCAAATATGAGCAACGTTTTGGTAATTCAGAATTTGGCTATCTGAAACTGAGCGCAGCGGAAGTGGCTTTAATCACCACGTTGTTATTGCgtgg                                                                         <  2:435243/250‑1 (MQ=255)                                                                                                                                                                              aGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCAAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGAACAGTGAGCGGATTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGCTCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCACAGACGc                                                              >  2:429706/1‑251 (MQ=255)                                                                                                                                                                                  aaTGGTGTAGTCACGGCCTGTAACCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAgg                                                          >  1:585175/1‑251 (MQ=255)                                                                                                                                                                                         tAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTg                                                   >  1:92596/1‑251 (MQ=255)                                                                                                                                                                                         tAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTg                                                   >  1:13930/1‑251 (MQ=255)                                                                                                                                                                                         cagtcccagCCTGTAAGAAGAAAACGAACCGTGAACCGGTCATGAGTCTGAGTGAATCCAAAGTGCAGGAACATCTGGATAAACTGGTCAAACGTCATTATCTACGCATATTGAGCGGTGTTGGTAATCGGGTCACCAAATATGATCAACGTTTTTGTAATTCAGAATGTGGCGAGCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTAGTGCTTGGTGCCCAGGCGCCAGGTGAACt                                                    <  2:673925/242‑1 (MQ=255)                                                                                                                                                                                          aGTCACGGCCTGTAATCAGAAAACGAACCGGGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTAGCTACGCACAGTTAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAAGCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTg                                                   <  1:564505/250‑1 (MQ=255)                                                                                                                                                                                               cGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAgc                                             >  1:37290/1‑251 (MQ=255)                                                                                                                                                                                                          cAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTcgcgcg                                  <  1:37561/251‑1 (MQ=255)                                                                                                                                                                                                           aGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGccgcgtcgcg                                    >  2:334366/1‑248 (MQ=255)                                                                                                                                                                                                              aaaCGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATg                              <  1:652834/251‑1 (MQ=255)                                                                                                                                                                                                                 cGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTAt                           >  1:291466/1‑251 (MQ=255)                                                                                                                                                                                                                   aaCCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGa                         <  1:302033/251‑1 (MQ=255)                                                                                                                                                                                                                     ccGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAAt                       <  1:329053/251‑1 (MQ=255)                                                                                                                                                                                                                             cGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGat               >  1:362192/1‑251 (MQ=255)                                                                                                                                                                                                                                            gTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTgg  >  1:238697/1‑249 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 GTTGATTGATGGACAATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGC‑TTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGG  >  NC_000913/1126074‑1126556

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