Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,662,523	A→G	intergenic (‑392/+60)	iscR ← / ← trmJ	isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,662,523	0	A	G	100.0%	129.5 / NA	38	intergenic (‑392/+60)	iscR/trmJ	isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/25);  total (13/25)

CGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGG  >  NC_000913/2662378‑2662663                                                                                                                                                  |                                                                                                                                             cgcgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGtt                                                                                                                                         <  1:660832/151‑1 (MQ=255)  gcggcgCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGtt                                                                                                                                         <  1:61536/150‑1 (MQ=255)     tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTg                                                                                                                                      <  1:242395/150‑1 (MQ=255)                    gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCTTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGa                                                                                                                      <  1:216679/151‑1 (MQ=255)                       gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGAt                                                                                                                     <  2:284533/149‑1 (MQ=255)                                    tGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGTACCCCTTTTCCGACCCGCGGTTGTGTTCAGTGTACGGCGGAGGAGACTCAtttc                                                                                                         >  1:169529/1‑146 (MQ=14)                                     gAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                      <  2:236011/150‑1 (MQ=255)                                       aCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                                    >  2:231440/1‑150 (MQ=255)                                         gCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGGCTGATAAGACGCATTACGCGTc                                                                                                    >  1:71013/1‑148 (MQ=255)                                             tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                              <  1:302600/150‑1 (MQ=38)                                              aTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                             <  2:555157/150‑1 (MQ=38)                                                        cGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                                                                     <  1:532612/148‑1 (MQ=17)                                                            tCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAt                                                                               <  1:282335/150‑1 (MQ=18)                                                                    cTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTg                                                                      >  2:49520/1‑151 (MQ=38)                                                                      gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                     <  1:595918/150‑1 (MQ=255)                                                                         ggCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                 >  2:166909/1‑151 (MQ=255)                                                                          gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt                                                                <  2:1055/151‑1 (MQ=255)                                                                           ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                                                               <  1:84268/151‑1 (MQ=255)                                                                              gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt                                                              <  2:212483/149‑1 (MQ=255)                                                                                 aaGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTct                                                         >  2:434051/1‑151 (MQ=255)                                                                                   gACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgc                                                       >  2:8572/1‑151 (MQ=255)                                                                                       cATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtgctgc                                                    >  2:509199/1‑150 (MQ=255)                                                                                           aCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAAt                                               <  1:370155/151‑1 (MQ=255)                                                                                           aCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAAt                                               <  2:351987/151‑1 (MQ=255)                                                                                                            aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                               >  2:653539/1‑150 (MQ=255)                                                                                                              ggCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCgcgc                            <  1:172465/151‑1 (MQ=255)                                                                                                               gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCa                           <  1:560249/151‑1 (MQ=255)                                                                                                               gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCa                           <  2:215007/151‑1 (MQ=255)                                                                                                               gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCa                           <  2:14162/151‑1 (MQ=255)                                                                                                                  gTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAgg                         >  1:643784/1‑150 (MQ=255)                                                                                                                    cGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGat                       >  2:617920/1‑150 (MQ=255)                                                                                                                      gATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATAtt                    >  1:277931/1‑151 (MQ=255)                                                                                                                          cggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                <  1:369878/151‑1 (MQ=255)                                                                                                                            gcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc              <  1:72175/151‑1 (MQ=255)                                                                                                                             cgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTcc             <  1:651270/151‑1 (MQ=255)                                                                                                                                gAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg          >  1:394304/1‑151 (MQ=255)                                                                                                                                   cGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCtt       <  2:609020/151‑1 (MQ=255)                                                                                                                                         aTCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCgg  <  2:75326/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             CGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGG  >  NC_000913/2662378‑2662663

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