Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,369,793	G→A	E35E (GAG→GAA)	pspE →	thiosulfate:cyanide sulfurtransferase (rhodanese)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,369,793	0	G	A	100.0%	125.7 / NA	37	E35E (GAG→GAA)	pspE	thiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/15);  total (22/15)

TACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGC  >  NC_000913/1369705‑1369889                                                                                         |                                                                                                 tACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAggg                                                                                       <  1:1042867/100‑1 (MQ=255)         ctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCa                                                                               <  1:1182559/100‑1 (MQ=255)               cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCGAGCGCAGTATCAGCAAGAACACGTTCAGGGGGCCATGAATATTc                                                                         >  2:1403335/1‑100 (MQ=255)               cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTc                                                                         >  1:893886/1‑100 (MQ=255)                      tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                  >  1:212680/1‑100 (MQ=255)                         cACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaagga                                                               <  1:532520/100‑1 (MQ=255)                          aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaa                                                              >  1:1537609/1‑100 (MQ=255)                             gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                           <  2:778640/100‑1 (MQ=255)                             gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                           >  2:1556117/1‑100 (MQ=255)                             gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                           <  1:1483673/100‑1 (MQ=255)                              cccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGa                                                          <  2:687192/100‑1 (MQ=255)                                cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaaa                                                        >  1:971903/1‑100 (MQ=255)                                     tCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGGATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcg                                                   <  1:1372796/100‑1 (MQ=255)                                       gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa                                                 >  2:1074776/1‑100 (MQ=255)                                         cgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAtt                                               >  1:424600/1‑100 (MQ=255)                                            tGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGcc                                            <  2:1277768/100‑1 (MQ=255)                                               acacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAcc                                         >  2:680059/1‑100 (MQ=255)                                                    ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt                                    >  1:1177456/1‑100 (MQ=255)                                                     gATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGtt                                   >  2:568228/1‑100 (MQ=255)                                                        cGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg                                <  1:1403343/100‑1 (MQ=255)                                                           tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                             <  2:1537599/100‑1 (MQ=255)                                                            gTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGTTGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTTCCACCGCCGTTCCGGATa                            >  2:546005/1‑100 (MQ=255)                                                               cGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGGGCATTGCCACCGCCGTTCCGGATaaaa                         >  1:615273/1‑100 (MQ=255)                                                                  gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      <  1:1333860/100‑1 (MQ=255)                                                                  gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      >  1:521002/1‑100 (MQ=255)                                                                  gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      >  2:655704/1‑100 (MQ=255)                                                                  gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      <  2:1160204/100‑1 (MQ=255)                                                                    tCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGac                    >  1:1158130/1‑100 (MQ=255)                                                                     ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                   >  1:1057124/1‑100 (MQ=255)                                                                     ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                   >  2:1566292/1‑100 (MQ=255)                                                                          gCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTg              <  2:174780/100‑1 (MQ=255)                                                                               aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgt         >  1:1106058/1‑100 (MQ=255)                                                                               aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg          >  2:1147696/1‑99 (MQ=255)                                                                                   gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt     >  2:1397196/1‑100 (MQ=255)                                                                                   gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt     <  1:577611/100‑1 (MQ=255)                                                                                    cAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt     >  1:1255134/1‑99 (MQ=255)                                                                                     aaGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTg   <  2:1336797/100‑1 (MQ=255)                                                                                      agaaCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc  >  1:1240599/1‑100 (MQ=255)                                                                                         |                                                                                                 TACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGC  >  NC_000913/1369705‑1369889

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