Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,543,800	C→T	intergenic (‑272/+32)	cysP ← / ← ucpA	thiosulfate/sulfate ABC transporter periplasmic binding protein CysP/oxidoreductase UcpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,543,800	0	C	T	100.0%	130.8 / NA	38	intergenic (‑272/+32)	cysP/ucpA	thiosulfate/sulfate ABC transporter periplasmic binding protein CysP/oxidoreductase UcpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/17);  total (21/17)

ATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTTATATACCGTTTGGTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTCCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTCGACTGGCGGGCAATGCTTTCC  >  NC_000913/2543566‑2544045                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      atcatATTTATAACGGACGTGAAAGGAGTAACGGTTTTATATACCGTTTGGTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGgga                                                                                                                                                                                                                                       <  2:253717/251‑1 (MQ=255)           taACGGACGTGAAAGGAGTAACGGTTTTATATACCGTTTGGTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACag                                                                                                                                                                                                                               >  2:233551/1‑249 (MQ=255)                gACGTGAAAGGAGTAACGGTTTTATATACCGTTTGGTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTg                                                                                                                                                                                                                          <  1:443840/249‑1 (MQ=255)                gACGTGAAAGGAGTAACGGTTTTATATACCGTTTGGTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTg                                                                                                                                                                                                                          >  2:443840/1‑249 (MQ=255)                   gTGAAAGGAGTAACGGTTTTATATACCGTTTGGTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATc                                                                                                                                                                                                                      >  1:5234/1‑250 (MQ=255)                                                  ggTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGtgtg                                                                                                                                                                                      >  1:20532/1‑251 (MQ=255)                                                     gATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGtgctg                                                                                                                                                                                   <  2:273351/251‑1 (MQ=255)                                                                tGAAAAGGGAATAAGAACTTGGTGCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAATGCGGAGCGTACACATAGTACGTGAGGATTTCGAGGACTGCCCGGGGCCAAAATTACAAATAAAATAGCCTGGTGAACTTAGTTCAAAACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAAGCAGATACCGTCGCTAACCGTCTCCGGCAGTGTGCTGCCGCCatcaa                                                                                                                                                                         <  2:361017/250‑1 (MQ=255)                                                                ggAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCatcaat                                                                                                                                                                        <  2:390432/250‑1 (MQ=255)                                                                       ggAATAAGAACTTGAATCGACCAGCGTAATTGTTGCCGTCAGGTTGGACACGAACAACGCGGCGACACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCAGGGGCCAAAATGAAAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCATTCCCACAACAGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTACACGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTc                                                                                                                                                                 >  2:171419/1‑251 (MQ=255)                                                                              gAACTTGGTTCACCAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTAc                                                                                                                                                          >  2:336658/1‑251 (MQ=255)                                                                                 cTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCcaca                                                                                                                                                                                                                                                 >  1:341287/1‑161 (MQ=255)                                                                                 cTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCcaca                                                                                                                                                                                                                                                 <  2:341287/161‑1 (MQ=255)                                                                                                           agtcagtTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCg                                                                                                                             >  1:453362/1‑251 (MQ=255)                                                                                                                   tGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGgc                                                                                                                     >  1:76613/1‑251 (MQ=255)                                                                                                                    ggACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGgcc                                                                                                                    <  2:167474/251‑1 (MQ=255)                                                                                                                             cAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCg                                                                                                           >  1:394809/1‑251 (MQ=255)                                                                                                                                   ggAGAAACCGGAGCGTACACATAGTACGTGAGGCTTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGCGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACt                                                                                                     >  1:426202/1‑251 (MQ=255)                                                                                                                                       aaaCCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCATGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAAGAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCa                                                                                                 <  2:89490/251‑1 (MQ=255)                                                                                                                                       aaaCCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTcc                                                                                                  >  2:239799/1‑250 (MQ=255)                                                                                                                                       aaaCCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTcc                                                                                                  >  1:254731/1‑250 (MQ=255)                                                                                                                                        aaCCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAg                                                                                                <  2:453362/251‑1 (MQ=255)                                                                                                                                           cGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCgg                                                                                             >  1:458327/1‑251 (MQ=255)                                                                                                                                                                gAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCgg                                                                        >  2:168554/1‑251 (MQ=255)                                                                                                                                                                         gAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCttt                                                               <  1:267658/251‑1 (MQ=255)                                                                                                                                                                          aGCACTGCCCGGTGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCGTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTTGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGTCGACGCATCGGGATTGCTTTc                                                              <  2:426202/251‑1 (MQ=255)                                                                                                                                                                            cACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGc                                                            >  1:298999/1‑251 (MQ=255)                                                                                                                                                                                  ccGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATc                                                                                          >  2:116194/1‑215 (MQ=255)                                                                                                                                                                                  ccGGGGCCAAAATGACAAATAAAATAGCCGGGTGAACTTAGTTCAAGACGAAATCCATCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTtc                                                      <  2:238897/251‑1 (MQ=255)                                                                                                                                                                                           aaaTGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCAc                                             >  2:89722/1‑251 (MQ=255)                                                                                                                                                                                             aTGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCg                                           >  1:403162/1‑251 (MQ=255)                                                                                                                                                                                               gACAAATAAAATAGCATGGTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGAc                                         <  2:394809/251‑1 (MQ=255)                                                                                                                                                                                                                 gTGAACTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTc                       <  1:89722/251‑1 (MQ=255)                                                                                                                                                                                                                     aCTTAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTCGACt                   >  2:313168/1‑251 (MQ=255)                                                                                                                                                                                                                        tAGTTCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTCGACTgg                 <  2:20532/250‑1 (MQ=255)                                                                                                                                                                                                                            tCAAGACGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCTAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTCGACTGGCggg             <  1:154036/250‑1 (MQ=255)                                                                                                                                                                                                                                  cGAAATCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTCGACTGGCGGGCAATGCt      >  1:91543/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTCCTTCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTCGACTGGCGGGCAATGCTTTcc  <  1:239799/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      ATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTTATATACCGTTTGGTGATTTGGAAGTTGAAAAGGGAATAAGAACTTGGTTCAACAGGCGTAATTGTTGCAGTCAGTTTGGACACGAACAGCGCGGAGAAACCGGAGCGTACACATAGTACGTGAGGATTTCGAGCACTGCCCGGGGCCAAAATGACAAATAAAATAGCCTGGTGAACTTAGTTCAAGACGAAATCCTCCCCACAAATGCAGGGAGGAAACAGAGGTGAATCAGATACCGACGCTAACCGTCTCCGGCAGTGTGCTGCCGCCATCAATCACATTCTGTGTACCGGTTAAATAGCTGGATTCATCCGATGCGAGGAAGGCCGCCAGTTCGCCGACTTCCAGCGGATCGGCGAGGCGACGCATCGGGATTGCTTTCGCCATTTCAGTCAGCACCGACTCTGGATCTTCCGGGTTCGACTGGCGGGCAATGCTTTCC  >  NC_000913/2543566‑2544045

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