Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,824,426	G→C	R676P (CGT→CCT)	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,824,426	0	G	C	100.0%	143.3 / NA	41	R676P (CGT→CCT)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (28/13);  total (28/13)

ATCCACCATGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCT  >  NC_000913/3824179‑3824668                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   aTCCACCATGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGa                                                                                                                                                                                                                                                 <  2:266287/251‑1 (MQ=255)    caccaTGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCg                                                                                                                                                                                                                                              >  1:144261/1‑251 (MQ=255)    caccaTGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCg                                                                                                                                                                                                                                              >  2:169747/1‑251 (MQ=255)              ccTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTg                                                                                                                                                                                                                                    >  1:210782/1‑251 (MQ=255)               cTGCCGTAATATCCGTGGCTACCAGAAAGCGCCAGAGAAGCTTATGGCTGTGGAATTGGATAAAGCGACGGCGCAGGAGT‑CACTCACCGACATCAAGGTGGAGATGTGCAATCATCAGGGTGCGCTGGCAAACCTGCCGGAGGCGATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAACGCTCCTGAcgtagtggatctgg                                                                                                                                                                                                                                   >  2:121088/1‑238 (MQ=255)                   cGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGCATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCgc                                                                                                                                                                                                                                >  1:239608/1‑249 (MQ=255)                         atCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTTGCGAATATCAt                                                                                                                                                                                                                         <  2:101249/251‑1 (MQ=255)                           ccGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCAAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAAAACCACGACTTCGAATATTCAACGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCCTgc                                                                                                                                                                                                                       >  2:137615/1‑251 (MQ=255)                                        aaaGAGCCAGAGACGTTTATGGCTGTGGACTGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGt                                                                                                                                                                                                          >  2:151788/1‑251 (MQ=255)                                               cAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGcc                                                                                                                                                                                                    >  2:50293/1‑250 (MQ=255)                                                  agaAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGAc                                                                                                                                                                                                <  1:55915/251‑1 (MQ=255)                                                   gaAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCCGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTCTTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGACTATCATGCGCAAACTCCGCGTGATGCCAGACg                                                                                                                                                                                               >  2:236212/1‑251 (MQ=255)                                                         tATGGCTGTGGAATGGGCTAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGAtt                                                                                                                                                                                          >  1:227917/1‑250 (MQ=255)                                                         tATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGAtt                                                                                                                                                                                          <  1:177580/250‑1 (MQ=255)                                                             gCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTCCAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCCTGCGCAAAATCCGCGTGATGCCAGACGTGATTaa                                                                                                                                                                                        >  2:190619/1‑248 (MQ=255)                                                               tgtgGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGAGGAGCTGTTCAATCATCAGGGTGCGCTGGCACACCTGTCGGCGGCAATTAACACCACGACTTCGAATCTTCCAAGTTTGACTACGGACGAGAAAGATGGTCGCGTCTCCAGCGCCTTTATTCGTCCGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCACAATCCGCGTGATGCCTGACGTGATTAAAGTcc                                                                                                                                                                                   >  1:23981/1‑250 (MQ=255)                                                                            aaaGAGACGGCGCAGGAGTTCA‑TCAACGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGACAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGGGATGCCAGACGTGATTAAAGTCACccgacaccgaa                                                                                                                                                                       >  2:121818/1‑250 (MQ=255)                                                                               gagaCGGCGCAGGCGTTCA‑TCACCGAAATCGCGGTGGCGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCCCGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAAtt                                                                                                                                                                    >  2:111422/1‑250 (MQ=255)                                                                                             gTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAAcc                                                                                                                                                     <  1:70442/251‑1 (MQ=255)                                                                                                    caCCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAacac                                                                                                                                               <  1:169747/251‑1 (MQ=255)                                                                                                         aaaTCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTAt                                                                                                                                          <  1:65611/251‑1 (MQ=255)                                                                                                                 gTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAAAACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAAACGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCAt                                                                                                                                  >  1:250579/1‑251 (MQ=255)                                                                                                                  ttgagaTGTTCAATCCCCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGCGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGCCCGTGTGCATCTGGCGACTATCATGCGCAAACTCCGCGTGATGCCAGCCGAGATTAAAGTCACCCGAAACCGACATTAATGTTTTCTGAACCCAACACGTTATGCACGCATc                                                                                                                                 >  2:168342/3‑251 (MQ=255)                                                                                                                         gTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATCTTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGaa                                                                                                                           >  2:17596/1‑250 (MQ=255)                                                                                                                                         gcgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCa                                                                                                          <  1:264348/251‑1 (MQ=255)                                                                                                                                                   aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTg                                                                                                >  1:251846/1‑251 (MQ=255)                                                                                                                                                              gcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCAt                                                                                     >  2:27784/1‑251 (MQ=255)                                                                                                                                                                   tAACACCACGACTTAGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCAAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGc                                                                                >  2:165407/1‑251 (MQ=255)                                                                                                                                                                    aaCACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGc                                                                                <  1:205989/250‑1 (MQ=255)                                                                                                                                                                              cTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGGCCACaa                                                                      >  2:84473/1‑250 (MQ=255)                                                                                                                                                                                    atatTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTCTATTCGTCTGACAGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGAGATGCAAGACGTGATTACAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAAAGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCAt                                                               >  1:264498/1‑251 (MQ=255)                                                                                                                                                                                     tatTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCAt                                                               >  1:271121/1‑250 (MQ=255)                                                                                                                                                                                            aGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGAGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTc                                                       >  1:50568/1‑251 (MQ=255)                                                                                                                                                                                             gTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCt                                                      >  2:104265/1‑251 (MQ=255)                                                                                                                                                                                             gTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCt                                                      >  1:70064/1‑251 (MQ=255)                                                                                                                                                                                                 gAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCg                                                   >  1:101082/1‑250 (MQ=255)                                                                                                                                                                                                               aaGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATATCGTTTCTGCGATTATTCGTACCGCa                                    <  1:189346/251‑1 (MQ=255)                                                                                                                                                                                                               aaGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCa                                    <  2:36324/251‑1 (MQ=255)                                                                                                                                                                                                                  aTGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGAt                                 <  1:34237/251‑1 (MQ=255)                                                                                                                                                                                                                         cgTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTg                          >  2:194839/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCGGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCt  <  1:166470/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   ATCCACCATGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCA‑TCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCT  >  NC_000913/3824179‑3824668

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