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%	108.1 / NA	34	*119* (TAG→TAA)	yffB	putative ArsC family reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/18);  total (16/18)

TCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTG  >  NC_000913/2591478‑2591727                                                                                                                              |                                                                                                                             tctGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                 <  1:68487/139‑1 (MQ=255)         ggcATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTAtt                                                                                                         <  2:857910/139‑1 (MQ=255)            aTTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAg                                                                                                      >  1:563846/1‑139 (MQ=255)              tAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCt                                                                                                    <  1:325634/139‑1 (MQ=255)                     tGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCt                                                                                                                           >  1:393033/1‑109 (MQ=255)                     tGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCt                                                                                                                           <  2:393033/109‑1 (MQ=255)                              tGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                 <  1:536886/110‑1 (MQ=255)                              tGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                 >  2:536886/1‑110 (MQ=255)                                cAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                  >  2:321777/1‑139 (MQ=255)                                    tATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgccgc                                                                              <  1:321777/139‑1 (MQ=255)                                     aTCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCt                                                                                                    <  1:132820/116‑1 (MQ=255)                                     aTCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCt                                                                                                    >  2:132820/1‑116 (MQ=255)                                           cGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTccc                                                                       <  2:765908/139‑1 (MQ=255)                                             tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTAt                                                                                                          <  2:762905/102‑1 (MQ=255)                                             tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTAt                                                                                                          >  1:762905/1‑102 (MQ=255)                                             tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgcc                                                                                <  2:909931/128‑1 (MQ=255)                                             tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgcc                                                                                >  1:909931/1‑128 (MQ=255)                                                aTTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGt                                                                  >  1:28563/1‑139 (MQ=255)                                                 ttGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTc                                                                 >  2:242140/1‑139 (MQ=255)                                                              ccGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGa                                                                                                       <  2:17404/88‑1 (MQ=255)                                                              ccGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGa                                                                                                       >  1:17404/1‑88 (MQ=255)                                                                 gTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATg                                                 >  2:731357/1‑139 (MQ=255)                                                                          tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                  >  2:179637/1‑97 (MQ=255)                                                                          tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                  <  1:179637/97‑1 (MQ=255)                                                                            ctgctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTttgt                                      <  1:480373/139‑1 (MQ=255)                                                                               ctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtgat                                                            <  2:265065/114‑1 (MQ=255)                                                                               ctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtgat                                                            >  1:265065/1‑114 (MQ=255)                                                                                              tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGttgat                                  >  1:617114/1‑125 (MQ=255)                                                                                              tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGttgat                                  <  2:617114/125‑1 (MQ=255)                                                                                                   ttATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATc               <  2:477490/139‑1 (MQ=255)                                                                                                       cAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGtt           >  2:68483/1‑139 (MQ=255)                                                                                                                ttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTg  <  2:563846/139‑1 (MQ=255)                                                                                                                    aTGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATc               >  1:790737/1‑122 (MQ=255)                                                                                                                    aTGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATc               <  2:790737/122‑1 (MQ=255)                                                                                                                              |                                                                                                                             TCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTG  >  NC_000913/2591478‑2591727

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