Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,874,741	T→C	T227A (ACA→GCA)	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,741	0	T	C	100.0%	110.0 / NA	29	T227A (ACA→GCA)	cysN	sulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/13);  total (16/13)

ATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTC  >  NC_000913/2874615‑2874867                                                                                                                               |                                                                                                                               aTGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTa                                                                                                                    <  2:601345/139‑1 (MQ=255)       gCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACgg                                                                                                              >  2:175874/1‑139 (MQ=255)        cGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGc                                                                                                             >  1:271821/1‑139 (MQ=255)          ttCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCAt                                                                                                           <  1:691714/139‑1 (MQ=255)           tCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATa                                                                                                          >  1:518220/1‑139 (MQ=255)                                                      tGCACCGGGAAGCGCTTTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACAACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCagaga                                                               >  1:467146/1‑139 (MQ=255)                                                      tGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCagaga                                                               >  1:548162/1‑139 (MQ=255)                                                           cGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGc                                                          >  1:311425/1‑139 (MQ=255)                                                            gggAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCATGTTGTCGCCTTCCAGTGCAGAGAGCGGca                                                         <  2:467146/139‑1 (MQ=255)                                                            gggAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGca                                                         <  2:271821/139‑1 (MQ=255)                                                             ggAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGcac                                                        >  1:485824/1‑139 (MQ=255)                                                                    cATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCgg                                                 <  2:548162/139‑1 (MQ=255)                                                                      ttGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGat                                               <  1:175874/139‑1 (MQ=255)                                                                                       acTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                              >  2:647529/1‑139 (MQ=255)                                                                                         tCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGc                            >  1:2599/1‑139 (MQ=255)                                                                                                     cACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCa                                                                   >  1:390523/1‑88 (MQ=255)                                                                                                     cACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCa                                                                   <  2:390523/88‑1 (MQ=255)                                                                                                     cACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAg                >  1:270313/1‑139 (MQ=255)                                                                                                     cACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAg                >  2:317102/1‑139 (MQ=255)                                                                                                     cACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAg                >  2:333596/1‑139 (MQ=255)                                                                                                           ttCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCaaa          <  2:518220/139‑1 (MQ=255)                                                                                                            tCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGat                                               >  2:143701/1‑101 (MQ=255)                                                                                                            tCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGat                                               <  1:143701/101‑1 (MQ=255)                                                                                                                gCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATc     <  1:333596/139‑1 (MQ=255)                                                                                                                   cTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTc  <  1:540361/139‑1 (MQ=255)                                                                                                                    ttCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCaa                  >  1:405537/1‑122 (MQ=255)                                                                                                                    ttCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCaa                  <  2:405537/122‑1 (MQ=255)                                                                                                                      cGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGat                                               >  2:183966/1‑91 (MQ=255)                                                                                                                      cGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGat                                               <  1:183966/91‑1 (MQ=255)                                                                                                                               |                                                                                                                               ATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTC  >  NC_000913/2874615‑2874867

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