Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,652,671	G→T	D60Y (GAT→TAT)	sseA →	3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,652,671	0	G	T	100.0%	141.1 / NA	42	D60Y (GAT→TAT)	sseA	3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (25/17);  total (25/17)

CATATTGATGACCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGC  >  NC_000913/2652536‑2652804                                                                                                                                        |                                                                                                                                      cATATTGATGACCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATca                                                                                                                                   <  1:21549/140‑1 (MQ=255)           aCCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCg                                                                                                                         >  2:478486/1‑140 (MQ=255)            cccGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGc                                                                                                                        >  2:735555/1‑140 (MQ=255)               ggAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTc                                                                                                                     >  2:426986/1‑140 (MQ=255)                gAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTcc                                                                                                                    >  2:309878/1‑140 (MQ=255)                             tCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCTGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCTGCACATGCTGCcg                                                                                                       >  2:239690/1‑140 (MQ=255)                             tCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCcg                                                                                                       >  1:139312/1‑140 (MQ=255)                                 tGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGcc                                                                                                   >  1:631463/1‑140 (MQ=255)                                     cGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGa                                                                                               >  1:485890/1‑140 (MQ=255)                                       cATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGaaa                                                                                             >  1:414310/1‑140 (MQ=255)                                                     gACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCg                                                                               <  2:414310/140‑1 (MQ=255)                                                     gACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCg                                                                               <  1:426986/140‑1 (MQ=255)                                                          gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCg                                                                          >  2:205648/1‑140 (MQ=255)                                                          gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCg                                                                          <  1:735555/140‑1 (MQ=255)                                                           aggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCg                                                                               >  2:548531/1‑134 (MQ=255)                                                           aggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCg                                                                               <  1:548531/134‑1 (MQ=255)                                                                    aCGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATg                                                                            <  1:214173/128‑1 (MQ=255)                                                                    aCGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATg                                                                            >  2:214173/1‑128 (MQ=255)                                                                     cGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCTGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCTGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCg                                                               <  1:239690/140‑1 (MQ=255)                                                                          cTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAAc                                                          >  2:111325/1‑140 (MQ=255)                                                                           tCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAAcc                                                         <  2:373208/140‑1 (MQ=255)                                                                                aGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGAt                                                    <  2:631463/140‑1 (MQ=255)                                                                                aGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGAt                                                    >  2:175688/1‑140 (MQ=255)                                                                                 gTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATa                                                   <  1:175688/140‑1 (MQ=255)                                                                                  tATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAAcc                                                         >  1:34261/1‑133 (MQ=255)                                                                                  tATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAAcc                                                         <  2:34261/133‑1 (MQ=255)                                                                                    tCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGc                                                <  1:111325/140‑1 (MQ=255)                                                                                       gAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCAcc                                             <  1:205648/140‑1 (MQ=255)                                                                                         aTGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTg                                           >  1:75418/1‑140 (MQ=255)                                                                                         aTGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTg                                           >  2:79609/1‑140 (MQ=255)                                                                                          tGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGa                                          <  2:139312/140‑1 (MQ=255)                                                                                              cATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGt                                      >  1:652746/1‑140 (MQ=255)                                                                                                  ttCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCt                                    <  2:173997/138‑1 (MQ=255)                                                                                                  ttCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCt                                    >  1:173997/1‑138 (MQ=255)                                                                                                       ggCGCAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGa                             >  2:512915/1‑140 (MQ=255)                                                                                                        gcgcAGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGGCTATGACGaa                            >  1:363243/1‑140 (MQ=255)                                                                                                            aGTGTTTTTTGATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTa                        >  1:634845/1‑140 (MQ=255)                                                                                                                      gATATCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTtctc              <  1:512915/140‑1 (MQ=255)                                                                                                                        tatCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAg            >  2:428634/1‑140 (MQ=255)                                                                                                                          tCGAAGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGcc          >  2:397115/1‑140 (MQ=255)                                                                                                                             aaGCGCTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCa       <  1:309878/140‑1 (MQ=255)                                                                                                                                  cTTTCTTATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGc  >  2:122663/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                      CATATTGATGACCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCCGCTTCCGCACATGCTGCCGCGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGC  >  NC_000913/2652536‑2652804

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