Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,591,603	G→A	*119* (TAG→TAA)	yffB →	putative ArsC family reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,591,603	0	G	A	100.0%	103.3 / NA	34	*119* (TAG→TAA)	yffB	putative ArsC family reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/17);  total (17/17)

ATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTG  >  NC_000913/2591470‑2591727                                                                                                                                      |                                                                                                                             aTGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTAt                                                                                                                         >  2:458741/1‑139 (MQ=255)                   cATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGa                                                                                                       <  1:86789/139‑1 (MQ=255)                                aaTGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                   >  2:926706/1‑114 (MQ=255)                                aaTGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                   <  1:926706/114‑1 (MQ=255)                                        cAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                  <  2:1021070/139‑1 (MQ=255)                                                    gTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCt                                                                      <  1:458741/139‑1 (MQ=255)                                                     tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTTGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCtt                                                                                       <  1:1340116/121‑1 (MQ=255)                                                     tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCtt                                                                                       >  2:1340116/1‑121 (MQ=255)                                                       cATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTccc                                                                       >  2:582344/1‑135 (MQ=255)                                                       cATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTccc                                                                       <  1:582344/135‑1 (MQ=255)                                                             tctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtga                                                             <  1:1172278/139‑1 (MQ=255)                                                                gCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtgatga                                                          >  1:970710/1‑139 (MQ=255)                                                                    gCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCa                                                      <  1:601834/139‑1 (MQ=255)                                                                      ccGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAgg                                                    >  1:1051791/1‑139 (MQ=255)                                                                                  tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCttt                                        >  2:67498/1‑139 (MQ=255)                                                                                  tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCttt                                        >  1:1243860/1‑139 (MQ=255)                                                                                      gctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTttgttg                                    >  1:26062/1‑139 (MQ=255)                                                                                               cAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                   <  1:1219513/51‑1 (MQ=255)                                                                                               cAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                   >  2:1219513/1‑51 (MQ=255)                                                                                                    aTTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCa                      >  2:442690/1‑139 (MQ=255)                                                                                                         aGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGa                 <  1:67498/139‑1 (MQ=255)                                                                                                           ttATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATc               <  2:1243860/139‑1 (MQ=255)                                                                                                              tCAGTAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACacaa                                                                                             >  1:168964/1‑58 (MQ=255)                                                                                                              tCAGTAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACacaa                                                                                             <  2:168964/58‑1 (MQ=255)                                                                                                                aGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGAt                <  1:659581/133‑1 (MQ=255)                                                                                                                aGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGAt                >  2:659581/1‑133 (MQ=255)                                                                                                                aGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGttt          >  1:910150/1‑139 (MQ=255)                                                                                                                        ttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTg  <  1:442690/139‑1 (MQ=255)                                                                                                                            aTGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                             >  2:1217591/1‑60 (MQ=255)                                                                                                                            aTGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                             <  1:1217591/60‑1 (MQ=255)                                                                                                                               aGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGaca                                                                                                >  2:855136/1‑38 (MQ=255)                                                                                                                               aGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGaca                                                                                                <  1:855136/38‑1 (MQ=255)                                                                                                                                ggTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgccgc                                                                              >  2:590864/1‑55 (MQ=255)                                                                                                                                ggTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgccgc                                                                              <  1:590864/55‑1 (MQ=255)                                                                                                                                      |                                                                                                                             ATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTG  >  NC_000913/2591470‑2591727

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