Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,874,431	G→A	A330V (GCG→GTG)	cysN ←	sulfate adenylyltransferase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,874,431	0	G	A	100.0%	166.3 / NA	49	A330V (GCG→GTG)	cysN	sulfate adenylyltransferase subunit 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/22);  total (27/22)

TAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATT  >  NC_000913/2874302‑2874563                                                                                                                                  |                                                                                                                                     tAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTgc                                                                                                                            >  2:562195/1‑140 (MQ=255)     gAATGCCATCAACACGCGCGCGCGTCTTCTTCCCGGAATTTTTGATGTCGTAACTCTGCCCGGGAAAAAGCGGTTGTTCCGCCATCCATACCACACCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTgccgcc                                                                                                                        <  2:492396/140‑1 (MQ=255)          ccATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCcagcag                                                                                                                   >  1:21674/1‑140 (MQ=255)             tCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCcagcagcag                                                                                                                >  1:456724/1‑140 (MQ=255)              cAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGa                                                                                                               >  1:818770/1‑140 (MQ=255)               aaCCCGCGCGCGCGTCTTCTTCCCGGCATTTTTGATGGCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCCCCATCCATACCACACCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGAt                                                                                                              <  2:386837/140‑1 (MQ=255)                  ccgcgcgcgggTTTTCTTACCGGAATTTTTGTTGTCGTACCCCTCCCGGGGAGAAGGGGGCTGCTCCCCCACCCTCACCACATCCCCGGCGGCGCTCTCCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGc                                                                                                           <  2:456724/139‑1 (MQ=255)                   cgcgcgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGcc                                                                                                          <  2:21674/140‑1 (MQ=255)                    gcgcgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATCCCACACCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCa                                                                                                         <  1:858752/140‑1 (MQ=255)                    gcgcgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCa                                                                                                         <  1:562195/140‑1 (MQ=255)                     cgcgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCAc                                                                                                        >  1:480625/1‑140 (MQ=255)                     cgcgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCAc                                                                                                        >  1:720906/1‑140 (MQ=255)                                    aCCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGAtct                                                                                         >  1:872357/1‑140 (MQ=255)                                       ggCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCAt                                                                                      >  1:375281/1‑140 (MQ=255)                                           aTTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATCCCACACCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGt                                                                                  <  2:771239/140‑1 (MQ=255)                                                 cattCGTAACTCCGCCCTGGAGAAAGCGGCTGTCCCGCCACCCTTCCCCCACCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCAc                                                                            <  2:790269/137‑1 (MQ=255)                                                     cGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAgg                                                                        >  2:487058/1‑140 (MQ=255)                                                      gTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAggg                                                                       >  2:258636/1‑140 (MQ=255)                                                       tAACTCTGCCCTGGAGAAAGCCGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGt                                                                      <  1:754527/140‑1 (MQ=255)                                                        aaCTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTg                                                                     <  1:487058/140‑1 (MQ=255)                                                            ctGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCg                                                                 >  1:237727/1‑140 (MQ=255)                                                              gCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCt                                                               >  1:746879/1‑140 (MQ=255)                                                              gCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCt                                                               >  1:628916/1‑140 (MQ=255)                                                               cccTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCtt                                                              >  1:10246/1‑140 (MQ=255)                                                               cccTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCtg                                                              >  1:790269/1‑139 (MQ=255)                                                                 cgggAGAAGGCGGCTGTCCCGCCATCCATCCCACACCCCCCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTtct                                                            <  2:638335/138‑1 (MQ=255)                                                                  tGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCCGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTtctc                                                           >  2:11069/1‑140 (MQ=255)                                                                         aGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCaaa                                                    >  2:759775/1‑140 (MQ=255)                                                                               tGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttc                                              <  1:759775/140‑1 (MQ=255)                                                                                        ccccatACCACCTCCCCCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATc                                     <  2:36581/138‑1 (MQ=255)                                                                                         tccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCa                                    >  2:498672/1‑140 (MQ=255)                                                                                                cacaTCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCaa                             <  2:375281/140‑1 (MQ=255)                                                                                                  caTCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCaaaa                           >  2:768797/1‑140 (MQ=255)                                                                                                    tCCACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGt                         >  1:144989/1‑140 (MQ=255)                                                                                                    ccccccGACGCGCTCTGCGCCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGt                         <  1:607830/136‑1 (MQ=255)                                                                                                     ccACCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCATGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTc                        >  1:292642/1‑140 (MQ=255)                                                                                                       aCCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCAc                      >  2:359011/1‑140 (MQ=255)                                                                                                       aCCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCAc                      >  2:1992/1‑140 (MQ=255)                                                                                                       aCCGACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCAc                      >  2:553207/1‑140 (MQ=255)                                                                                                          gACGCGCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGAt                   >  2:698099/1‑140 (MQ=255)                                                                                                            cgcgCTCTGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATcc                 >  1:58473/1‑140 (MQ=255)                                                                                                               gcTCTGCACCGCCGGTAATACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcg              <  1:803008/140‑1 (MQ=255)                                                                                                                 tccgCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgcg            <  1:11069/137‑1 (MQ=255)                                                                                                                   tGCACCGCCGGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGAc          >  1:16981/1‑140 (MQ=255)                                                                                                                       ccgccgGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGttt      <  2:480625/140‑1 (MQ=255)                                                                                                                       ccgccgGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGttt      <  2:720906/140‑1 (MQ=255)                                                                                                                       ccgccgGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGttt      <  2:10246/140‑1 (MQ=255)                                                                                                                         gccgGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGa    <  2:16981/140‑1 (MQ=255)                                                                                                                         gccgGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGGCACGATCCGCGCGACGTTTGa    <  1:553207/140‑1 (MQ=255)                                                                                                                           cgGTAACACTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGAtt  <  1:768797/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     TAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATT  >  NC_000913/2874302‑2874563

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: ‑