Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,369,805	G→A	G39G (GGG→GGA)	pspE →	thiosulfate sulfurtransferase PspE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,369,805	0	G	A	100.0%	54.3 / NA	18	G39G (GGG→GGA)	pspE	thiosulfate sulfurtransferase PspE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/8);  total (10/8)

CGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAG  >  NC_000913/1369682‑1369934                                                                                                                            |                                                                                                                                  cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAat                                                                                                                          >  1:44234/1‑133 (MQ=255)  gCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAata                                                                                                                         <  2:44234/133‑1 (MQ=255)      aCATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTcc                                                                                                                     <  1:35790/133‑1 (MQ=255)           tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTga                                                                                                                <  1:62950/133‑1 (MQ=255)              aaaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCTGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTTAGGGCGCCATCAATATTCCCCTgacgg                                                                                                             >  1:514781/1‑133 (MQ=255)                     gCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGaa                                                                                                      >  2:475515/1‑133 (MQ=255)                         aCTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAagag                                                                                                  <  1:345967/133‑1 (MQ=255)                            tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                                                               <  1:449555/133‑1 (MQ=255)                                         gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgccg                                                                                  >  1:410604/1‑133 (MQ=255)                                                           ttCGCCGCCGAATACTGGTTCGCTGTTAGTATTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                                                                <  1:225396/133‑1 (MQ=255)                                                                                    ttCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCgg                                       >  2:332173/1‑133 (MQ=255)                                                                                       gttttCCAGAGCAGATTCAACAAGTGGCCGTTAGGGGAGCCCTAAATATCCCCCGGAAGGAAGGGAAAGAGCGCTTTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACg                                    <  2:209270/130‑1 (MQ=255)                                                                                       gtgtTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACg                                    <  2:410604/133‑1 (MQ=255)                                                                                           tCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGAGCCATAAATATCCCCCGAAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAg                                <  2:468746/133‑1 (MQ=255)                                                                                                  cAGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATCCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGgc                         <  2:293905/133‑1 (MQ=255)                                                                                                          gCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGAAATGCCGGACGGCAGTCAGGGGAAGCaaaa                 >  1:147940/1‑133 (MQ=255)                                                                                                              gagCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAgaga             >  2:266696/1‑133 (MQ=255)                                                                                                               agCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTTTATTGCAATGCCTGACTCCAGTCAGGGCAATCAAAAGAGAt            >  2:151355/1‑133 (MQ=255)                                                                                                               agCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGAt            >  1:449538/1‑133 (MQ=255)                                                                                                               agCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGAt            >  2:364895/1‑133 (MQ=255)                                                                                                               agCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCAAGTAAGTGCAAGCAAAAGAGAt            >  1:169828/1‑133 (MQ=255)                                                                                                                         gggagCCTTCAATTGTCCCTTAAAGAAAGTAAAAGACTGCATGCCAACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATCCCGGACGCCATTCAGGCCAAGCAAAAGAGATCCTTaccccg  <  2:514781/133‑6 (MQ=255)                                                                                                                            |                                                                                                                                  CGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAG  >  NC_000913/1369682‑1369934

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