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%	140.2 / NA	42	R676P (CGT→CCT)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (28/14);  total (28/14)

CACCATGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCC  >  NC_000913/3824182‑3824664                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               caccaTGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCg                                                                                                                                                                                                                                          >  1:107188/1‑251 (MQ=255)               ccGTAATATCCGTGGCTACCAGAAAGAGCCAGAGACGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAACTCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCg                                                                                                                                                                                                                             >  2:123671/1‑250 (MQ=255)                        ccgtttcTACCAGAAATAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGCGTTCATCACCGAAATCAAGGTGGAGATGTGCAATCATCAGGGTGCGTTGTCAATCCTTACGGCGTCAATTAACACCACGACTTCGAATATTTAAAGTTTGAATACGGTAGAGAAAGATGGTCGGGTCTACAGCGTCTTTATTCGTTTGACCGCTCCTTACCGTGTGCATCTGGCCGATATCATgc                                                                                                                                                                                                                   <  2:251275/245‑1 (MQ=255)                        ccGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATg                                                                                                                                                                                                                    <  1:188247/250‑1 (MQ=255)                         cGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATgcg                                                                                                                                                                                                                  >  2:97594/1‑251 (MQ=255)                                 ccAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATcc                                                                                                                                                                                                          <  2:360979/251‑1 (MQ=255)                                           ccAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGcc                                                                                                                                                                                                >  1:171394/1‑251 (MQ=255)                                               agaAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGAc                                                                                                                                                                                            <  2:36825/251‑1 (MQ=255)                                               agaAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGAc                                                                                                                                                                                            >  1:205846/1‑251 (MQ=255)                                               agaAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGAc                                                                                                                                                                                            <  1:76252/251‑1 (MQ=255)                                                  aGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTg                                                                                                                                                                                         >  2:90632/1‑251 (MQ=255)                                                  aGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTg                                                                                                                                                                                         >  2:110248/1‑251 (MQ=255)                                                      tATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGAAATTAACACCACGACTTCGAATAATCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACATGAtt                                                                                                                                                                                      <  2:259908/250‑1 (MQ=255)                                                       aTGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGATGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGCGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTa                                                                                                                                                                                     >  2:181847/1‑250 (MQ=255)                                                         ggCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTaa                                                                                                                                                                                    >  2:137271/1‑249 (MQ=255)                                                             gtgGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCa                                                                                                                                                                               <  1:139730/250‑1 (MQ=255)                                                                         aaaGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACccgaaaccgaa                                                                                                                                                                   >  2:37689/1‑250 (MQ=255)                                                                             agaCGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGGCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGGGATGCCAGACGTGATTAAAGTCACCCTAAACCGAAATTaa                                                                                                                                                              <  2:27089/251‑1 (MQ=255)                                                                                   cgcAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGtttt                                                                                                                                                        >  2:79109/1‑251 (MQ=255)                                                                                   cgcAGGAGTTCATCACCGAAATCAAGGAGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGtttt                                                                                                                                                        >  2:30349/1‑251 (MQ=255)                                                                                     cAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTAt                                                                                                                                                      >  1:55328/1‑251 (MQ=255)                                                                                      tggAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATg                                                                                                                                                     >  2:119470/2‑251 (MQ=255)                                                                                              atcaCCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAac                                                                                                                                             >  1:340143/1‑251 (MQ=255)                                                                                                    gAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGtt                                                                                                                                        >  1:125892/1‑250 (MQ=255)                                                                                                      aaTCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATg                                                                                                                                     >  1:313677/1‑251 (MQ=255)                                                                                                             gTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACAGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCAt                                                                                                                              <  1:91464/251‑1 (MQ=255)                                                                                                               ggAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCt                                                                                                                            >  1:141091/1‑251 (MQ=255)                                                                                                                gagaTGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTg                                                                                                                           <  2:205846/251‑1 (MQ=255)                                                                                                                  gaTGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCg                                                                                                                         <  1:123061/251‑1 (MQ=255)                                                                                                                                          ggCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTg                                                                                                 <  1:69373/251‑1 (MQ=255)                                                                                                                                                         ggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGc                                                                                   >  2:280028/1‑250 (MQ=255)                                                                                                                                                             aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATgg                                                                               >  1:171006/1‑250 (MQ=255)                                                                                                                                                             aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATgg                                                                               >  2:247/1‑250 (MQ=255)                                                                                                                                                                      acGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTcc                                                                      >  1:196327/1‑250 (MQ=255)                                                                                                                                                                           ttCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAAc                                                                >  1:107704/1‑251 (MQ=255)                                                                                                                                                                           gtCGAATAGTCAAAGGTTGAATACGGAAGAGAAAGAGGGGCGCGTCGACAGCGCCGTGAGTCGTCGGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGAAGGCTGATCTGACCGTCTGCATGGAGCAGGTACACAAAc                                                                <  2:34/250‑1 (MQ=255)                                                                                                                                                                                          tttGAATACGGAAGCGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACAGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGCCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTg                                                 >  2:136723/1‑251 (MQ=255)                                                                                                                                                                                                    gAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCg                                       >  2:89769/1‑251 (MQ=255)                                                                                                                                                                                                            aGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGc                                 >  2:161740/1‑249 (MQ=255)                                                                                                                                                                                                                ggTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGc                           >  1:282331/1‑251 (MQ=255)                                                                                                                                                                                                                ggTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATAAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATg                            <  2:239357/250‑1 (MQ=255)                                                                                                                                                                                                                                         gTCTGACCGCTCCTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGcc  >  2:368748/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               CACCATGAATCCTGCCGTAATATCCGTGGCTACCAGAAAGAGCCAGAGAAGTTTATGGCTGTGGAATGGGATAAAGAGACGGCGCAGGAGTTCATCACCGAAATCAAGGTGGAGATGTTCAATCATCAGGGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCC  >  NC_000913/3824182‑3824664

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