Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,022	T→C	E167E (GAA→GAG)	ettA ←	energy‑dependent translational throttle protein EttA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,022	0	T	C	100.0%	198.5 / NA	64	E167E (GAA→GAG)	ettA	energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/36);  total (28/36)

TCCAGCCAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCA  >  NC_000913/4629777‑4630262                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 cccagccagGAGGAGTAGTTGCCTTCCCACGGAATACCTTCGCCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTcacc                                                                                                                                                                                                                                             <  2:155804/250‑1 (MQ=255)     gccagGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCg                                                                                                                                                                                                                                         <  1:199109/251‑1 (MQ=255)        aggaggAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCgg                                                                                                                                                                                                                                        >  1:258189/1‑249 (MQ=255)         ggaggaGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGgaga                                                                                                                                                                                                                                     <  2:103532/251‑1 (MQ=255)              agtagtTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGTGAGGtg                                                                                                                                                                                                                                 >  2:154512/1‑249 (MQ=255)               gtagtTGCCTTCCCACGGAATACCTCCACCGCGGTCAAGTTCGAGGATCCACCCGGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTGCCAGCCAGGCCGCGGATTCGGCATCCAGGTGGTTGGGCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACTACGCTACGCGACGACTCTCACCACCGGAGAGGTTAg                                                                                                                                                                                                                               <  2:28930/251‑1 (MQ=255)                  gttgccTTCCCACGGAATTCCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGa                                                                                                                                                                                                                            <  2:2936/247‑1 (MQ=255)                         tCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTcg                                                                                                                                                                                                                      <  2:227399/250‑1 (MQ=255)                          cccACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTcgc                                                                                                                                                                                                                     <  1:52964/250‑1 (MQ=255)                           ccACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGt                                                                                                                                                                                                                   >  1:105588/1‑251 (MQ=255)                             aCGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGTTTCGTCGAGCAGCAGCATGTCTGGTTTTGCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTcc                                                                                                                                                                                                                 <  1:258537/251‑1 (MQ=255)                              cGGAAGACCTTCACCGCGGTCAAGTTCGAGGAGCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCTCGGTGCCCTCGAAGTCGGGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCTAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGTTTAGCGAGTTTCGCGTccc                                                                                                                                                                                                                <  2:231651/251‑1 (MQ=255)                                 aaTACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGt                                                                                                                                                                                                             <  2:132892/251‑1 (MQ=255)                                    aCCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTcc                                                                                                                                                                                                           <  1:140105/250‑1 (MQ=255)                                       ttCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCCGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGc                                                                                                                                                                                                        >  1:201186/1‑250 (MQ=255)                                               ggTCAAGTTCGAGGATCCCGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGATAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAgc                                                                                                                                                                                               >  2:170100/1‑251 (MQ=255)                                                        cGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATccgcc                                                                                                                                                                                      >  1:167533/1‑251 (MQ=255)                                                         gAGGATCCAGCTTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATccgccg                                                                                                                                                                                     <  1:158239/251‑1 (MQ=255)                                                         gAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATccgccg                                                                                                                                                                                     >  2:75930/1‑251 (MQ=255)                                                            gATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCAc                                                                                                                                                                                  <  1:87209/251‑1 (MQ=255)                                                                 aGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGGTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGgaga                                                                                                                                                                                                                                     <  1:130511/195‑1 (MQ=255)                                                                 aGCCTGCAACGTTATCGAGGAAGTAACGGTAGTGGGTAATCGACAACACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGCCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCCCAACGCTACGCGACGACGCTCACCACCGgagc                                                                                                                                                                                                                                     >  2:130511/1‑194 (MQ=255)                                                                       cAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGt                                                                                                                                                                       <  1:28960/251‑1 (MQ=255)                                                                         aCGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTAc                                                                                                                                                                     <  1:30457/251‑1 (MQ=255)                                                                                     aaGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCACCCAGGCCACGGATTCGGCCTCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCCTGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATg                                                                                                                                                         >  2:118493/1‑251 (MQ=255)                                                                                           cGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACcgtc                                                                                                                                                   <  1:49688/251‑1 (MQ=255)                                                                                                gTGGGTCATAGCCACCACGGTGCCCTCGAAGTAGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGAACGTCGTg                                                                                                                                                >  2:51042/1‑249 (MQ=255)                                                                                                   ggTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCt                                                                                                                                           >  2:157241/1‑251 (MQ=255)                                                                                                       aTCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCAACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAAt                                                                                                                                       >  1:13070/1‑251 (MQ=255)                                                                                                          gCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGGCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGAt                                                                                                                                    >  2:171618/1‑251 (MQ=255)                                                                                                                    tGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTACAGCAGCCGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGCGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCCGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTcccgc                                                                                                                          >  2:49412/1‑248 (MQ=255)                                                                                                                     gccctCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTGGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTCTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCTCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGa                                                                                                                          <  1:43097/246‑1 (MQ=255)                                                                                                                     gccctCGAAGTCGGGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGAc                                                                                                                         <  2:143847/247‑1 (MQ=255)                                                                                                                      ccctCGAAGTCTTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCCGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGCGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGACGTGAGCCTGCATGATCTCTTCCAGACg                                                                                                                        >  2:37440/4‑251 (MQ=255)                                                                                                                      ccctCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACg                                                                                                                        <  1:241194/248‑1 (MQ=255)                                                                                                                      cccacGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCTCGGATTCGGCATCCAGTTGGTTGGTCGGTTCGTTGAGCAGCAGCATGTCTGGTTTTTCCATCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGATGTTAGCTATTTTCGCGTCCGATTCCGGCAGACGTATCTCTTCCGCCGCACGCTCCAGCTGTGCGTTCAGGTTAGGACTGTCTTGAGCCTGTATAATCGCTTCCAGACg                                                                                                                        <  2:122173/247‑1 (MQ=255)                                                                                                                       cctCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCCTGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACgg                                                                                                                       >  2:26288/3‑251 (MQ=255)                                                                                                                          cGAAGTCGTGCCGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGCATGATCTCTTCCAGACGGCCt                                                                                                                    >  1:194062/1‑251 (MQ=255)                                                                                                                           gAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCAGCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCtt                                                                                                                   <  2:105588/251‑1 (MQ=255)                                                                                                                            aaGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTg                                                                                                                  <  1:154512/251‑1 (MQ=255)                                                                                                                                cGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCTCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCGCTTCCAGACGGCCTTGTTCa                                                                                                              <  1:215567/251‑1 (MQ=255)                                                                                                                                         aGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCCGCGGCCAGc                                                                                                     >  1:53834/1‑251 (MQ=255)                                                                                                                                                cagccagGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAgcg                                                                                                                                                                                              <  1:198084/155‑1 (MQ=255)                                                                                                                                                cagccagGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAgcg                                                                                                                                                                                              >  2:198084/1‑155 (MQ=255)                                                                                                                                                          aCGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGAAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCACCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATc                                                                                    >  1:245877/1‑251 (MQ=255)                                                                                                                                                            ggATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGCATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCg                                                                                  >  1:5277/1‑251 (MQ=255)                                                                                                                                                                 cGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAACGGCCAGCTGGTCAAAATCGGCATCCGGATCg                                                                             <  2:48910/251‑1 (MQ=255)                                                                                                                                                                         ggtgggttgggCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACACCGCTACGCGACGACGCTCACCACCGGAGAGTTTAGCGAGTTTCGCGTCCCAGTCCGGCACACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGTCCGTCGTGAGCCTGAATGATCTCTGCCAAACGGCCTCGTTCAGCGGCCAGCATTTCAAAATCGGCATCCGGATCGGCTTACAg                                                                     <  2:244164/247‑1 (MQ=255)                                                                                                                                                                                tggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGGTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATac                                                              >  2:20169/1‑251 (MQ=255)                                                                                                                                                                                        tcgtcgAGCAGCAGCATGTCTGGTTGTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCAt                                                       <  1:201585/250‑1 (MQ=255)                                                                                                                                                                                          gtcgAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCa                                                    <  1:64142/251‑1 (MQ=255)                                                                                                                                                                                            cgAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAg                                                   <  1:203093/250‑1 (MQ=255)                                                                                                                                                                                             gAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCAGATCGGCGTACAGCGCATACACTTCATCCAGGc                                                 >  2:188428/1‑251 (MQ=255)                                                                                                                                                                                             gAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCAACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCCTCCAGGc                                                 >  2:252366/1‑251 (MQ=255)                                                                                                                                                                                             gAGCAGCAGCATGTCTGGTTTTCCCAGCAGCAGGCGGCAAAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTGTCGCGGCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAGATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGc                                                 <  2:201186/251‑1 (MQ=255)                                                                                                                                                                                              agcagaagcATGTCTGGTTTTTCCAGCAGCAGGCGGCACCACGCTACGCGACGCCGCTCCCCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGAAGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCCGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCCTACAATTCCTCCAGGCg                                                >  2:38751/1‑251 (MQ=255)                                                                                                                                                                                                  gcagcaTGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTc                                            >  2:171732/1‑251 (MQ=255)                                                                                                                                                                                                  cgcgcctGTCTGGTTTTTCCAGCAGCAGGCAGCCCACCGCTACGCGACGACGCCCCCCACCGGAGAGGTTAGCGATTTTCGCGCCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGTTTATGACCGCCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCCTCCAGGCGttgc                                            <  2:78619/245‑3 (MQ=255)                                                                                                                                                                                                   cagcaTGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTGTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCGCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCa                                           <  1:119954/251‑1 (MQ=255)                                                                                                                                                                                                            tGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGAAGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAAACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACCTCATCCAGGCGCTTCAGGGCGTTa                                   >  2:213717/1‑250 (MQ=255)                                                                                                                                                                                                                       gcagcagGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCATTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTGATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGaa                       <  1:136378/251‑1 (MQ=255)                                                                                                                                                                                                                          gcagGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAAcc                    <  2:5766/251‑1 (MQ=255)                                                                                                                                                                                                                                   acaACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCCGAAACCGCttcttc            >  2:56877/1‑250 (MQ=255)                                                                                                                                                                                                                                            cgcgAAGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTca  <  2:48648/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 TCCAGCCAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCA  >  NC_000913/4629777‑4630262

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