Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,162	C→T	H255Y (CAC→TAC)	spoT →	bifunctional (p)ppGpp synthase/hydrolase SpoT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,823,162	0	C	T	100.0%	228.4 / NA	62	H255Y (CAC→TAC)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/29);  total (33/29)

TTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTT  >  NC_000913/3822916‑3823404                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   gtAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGAGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTGTGCAGGAAGCGGGAATAGCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCATCGTTTTTACTc                                                                                                                                                                                                                                                <  2:299642/250‑1 (MQ=255)       aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCCAAGAGCAGCGTTTTTACTCGAt                                                                                                                                                                                                                                             >  2:342667/1‑248 (MQ=255)       aTCCACCACACTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCTTGCCGCGTCAGTGGTCGCGAGAAGCAGCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGAt                                                                                                                                                                                                                                             <  1:342667/248‑1 (MQ=255)         ccaccaCATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTCTCCCCACCGTTATCGCGTAATCAACGAAGTGGAGACAGCCGCGCGCGGAAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCCGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATgg                                                                                                                                                                                                                                        >  1:201192/1‑251 (MQ=255)         ccaccaCATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATg                                                                                                                                                                                                                                         <  2:8530/250‑1 (MQ=255)                   aaaaCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTGCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTAct                                                                                                                                                                                                                               <  2:95059/250‑2 (MQ=255)                      aCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCttt                                                                                                                                                                                                                           >  2:3946/1‑251 (MQ=255)                        cGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTcc                                                                                                                                                                                                                         >  1:102715/1‑251 (MQ=255)                           aCTCGAAGAGCTGGGTTTAGAGGCACTGTCTCCCAACCGTTATCGCGTAATCAAAGCCGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCcgcg                                                                                                                                                                                                                      >  2:157906/1‑251 (MQ=255)                                 agagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCg                                                                                                                                                                                                                >  2:311608/1‑251 (MQ=255)                                 agagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTGTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCg                                                                                                                                                                                                                <  2:336707/251‑1 (MQ=255)                                   agCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTc                                                                                                                                                                                                              >  2:178325/1‑251 (MQ=255)                                                ggCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGGCAATGATTCTGaca                                                                                                                                                                                                 <  2:114918/251‑1 (MQ=255)                                                      gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGt                                                                                                                                                                                            <  2:267905/250‑1 (MQ=255)                                                            cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcgcg                                                                                                                                                                                     >  2:291842/1‑251 (MQ=255)                                                               ccGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGc                                                                                                                                                                                  <  2:340722/251‑1 (MQ=255)                                                               ccGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGc                                                                                                                                                                                  <  1:256199/251‑1 (MQ=255)                                                                      cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGccc                                                                                                                                                                           >  2:299393/1‑250 (MQ=255)                                                                           aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTGATCGCGTGCTGGGCCAGATGc                                                                                                                                                                      <  2:344534/251‑1 (MQ=255)                                                                               aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGcac                                                                                                                                                                    >  2:108699/1‑249 (MQ=255)                                                                                           aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGcc                                                                                                                                                      <  2:120443/251‑1 (MQ=255)                                                                                                   gcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGAGCATGGACATCTACGCTTTTCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCgg                                                                                                                                              <  2:171137/251‑1 (MQ=255)                                                                                                         cacccGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgcg                                                                                                                                        >  2:53189/4‑251 (MQ=255)                                                                                                              gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                                                                                   <  1:121676/251‑1 (MQ=255)                                                                                                              gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTAGCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                                                                                   <  1:147488/251‑1 (MQ=255)                                                                                                                aaatagaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                                                                                 <  2:83841/247‑1 (MQ=255)                                                                                                                    agaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                                             >  1:9772/1‑251 (MQ=255)                                                                                                                          tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGcc                                                                                                                       >  2:312314/1‑251 (MQ=255)                                                                                                                              gaagaTTGTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTc                                                                                                                   >  2:207936/1‑251 (MQ=255)                                                                                                                              gaagaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCAGGGCCGCGTGAAAGACTATATCGCCATTc                                                                                                                   >  2:231907/1‑251 (MQ=255)                                                                                                                                   ttctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCa                                                                                                                 <  2:68099/248‑1 (MQ=255)                                                                                                                                   ttctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCa                                                                                                                 >  1:68099/1‑248 (MQ=255)                                                                                                                                     ctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGc                                                                                                            >  2:301968/1‑251 (MQ=255)                                                                                                                                           gAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACAGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACAACTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                                        <  2:245949/249‑1 (MQ=255)                                                                                                                                               tCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTAt                                                                                                  >  1:252601/1‑251 (MQ=255)                                                                                                                                                cGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGCGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATCTCGCCATTCCAAAAGCGAACGGCTATc                                                                                                 >  2:238472/1‑251 (MQ=255)                                                                                                                                                               ggAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGGGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGGACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCt                                                                                  <  2:183528/251‑1 (MQ=255)                                                                                                                                                                 aaGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                                                <  1:263597/251‑1 (MQ=255)                                                                                                                                                                   gCGGGAATACCTTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTTGTTGCACACCTCgat                                                                              <  1:284955/251‑1 (MQ=255)                                                                                                                                                                   gCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgat                                                                              <  1:284945/251‑1 (MQ=255)                                                                                                                                                                        aaTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                                         <  1:299393/251‑1 (MQ=255)                                                                                                                                                                           aCCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGcc                                                                      <  1:12211/251‑1 (MQ=255)                                                                                                                                                                            ccGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGccc                                                                     >  1:306910/1‑251 (MQ=255)                                                                                                                                                                              gTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGc                                                                   >  1:92128/1‑251 (MQ=255)                                                                                                                                                                                 ccGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                                                >  1:276509/1‑251 (MQ=255)                                                                                                                                                                                 acgcgTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                                                >  1:276506/2‑251 (MQ=255)                                                                                                                                                                                        aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGTTGCTCAAAGAGCCGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                                         <  2:129995/251‑1 (MQ=255)                                                                                                                                                                                         gTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCg                                                        <  1:3946/251‑1 (MQ=255)                                                                                                                                                                                                       cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTcc                                              >  2:287739/1‑247 (MQ=255)                                                                                                                                                                                                           tttATTCGATTTACTGCACAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCCCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGCCTATATCGCCATTCCAAACGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGCGGTCCAGCTCCGt                                      >  1:265122/1‑251 (MQ=255)                                                                                                                                                                                                            ttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                                      >  2:235036/1‑250 (MQ=255)                                                                                                                                                                                                                 cGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGaa                                <  1:290164/251‑1 (MQ=255)                                                                                                                                                                                                                        cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgg                         >  1:144761/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaTGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAg                    <  2:144761/251‑1 (MQ=255)                                                                                                                                                                                                                               aTGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc                      >  1:148083/1‑247 (MQ=255)                                                                                                                                                                                                                                         aaGAGCAGAGTTTTTACTCGATCATGGACATCTAAGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAACGCGCACGGCTATCAGTCTTTGCACACCTCGCTGATCGGCCCGCACGGTGTGCCGGTTGAGGTACAGATCCGTACCGCAGATATGGACCAGATGGCGgagagg        >  2:216438/1‑249 (MQ=255)                                                                                                                                                                                                                                          agagCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGaa                                >  1:271728/1‑226 (MQ=255)                                                                                                                                                                                                                                          agagCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGCTATGGACCAGATGGCGGAGATgg       >  2:119111/1‑251 (MQ=255)                                                                                                                                                                                                                                          agagCAGCGTTTTTACTCGATCATGGACATCTACGCTTCCCGCGTGATCGTCAATGATTCTGCCACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGaa                                <  2:271728/226‑1 (MQ=255)                                                                                                                                                                                                                                           gagCAGCTTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATggg      <  2:32806/251‑1 (MQ=255)                                                                                                                                                                                                                                             gcagcGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGgtg    >  2:94876/1‑251 (MQ=255)                                                                                                                                                                                                                                               agcGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGtt  >  1:350226/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   TTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTT  >  NC_000913/3822916‑3823404

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