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

CCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGA  >  NC_000913/2591476‑2591737                                                                                                                                |                                                                                                                                       ccTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                             <  2:1006990/139‑1 (MQ=255)             cATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGccc                                                                                                                          <  2:325292/130‑1 (MQ=255)             cATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGccc                                                                                                                          >  1:325292/1‑130 (MQ=255)                  aTGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGa                                                                                                            <  2:165626/139‑1 (MQ=255)                       tGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTa                                                                                                                     <  2:393041/125‑1 (MQ=255)                       tGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTa                                                                                                                     >  1:393041/1‑125 (MQ=255)                         aaaTGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACacaaca                                                                                                     >  1:354441/1‑139 (MQ=255)                           aTGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGaca                                                                                                          >  2:815571/1‑132 (MQ=255)                           aTGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGaca                                                                                                          <  1:815571/132‑1 (MQ=255)                                       aTCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                           >  2:251711/1‑103 (MQ=255)                                       aTCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                           <  1:251711/103‑1 (MQ=255)                                        tCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                           <  2:689254/102‑1 (MQ=255)                                        tCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                           >  1:689254/1‑102 (MQ=255)                                              gTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCt                                                                                <  1:628122/139‑1 (MQ=255)                                                         tGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                           <  2:1212657/85‑1 (MQ=255)                                                         tGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                           >  1:1212657/1‑85 (MQ=255)                                                         tGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtgatg                                                                     >  1:819505/1‑139 (MQ=255)                                                          gCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtgatga                                                                    <  2:68579/139‑1 (MQ=255)                                                                     aaGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATg                                                           >  1:641356/1‑137 (MQ=255)                                                                     aaGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATg                                                           <  2:641356/137‑1 (MQ=255)                                                                     aaGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGcc                                                         <  1:321632/139‑1 (MQ=255)                                                                     aaGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGcc                                                         <  2:354441/139‑1 (MQ=255)                                                                         cTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGt                                                                            >  2:810555/1‑116 (MQ=255)                                                                         cTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGt                                                                            <  1:810555/116‑1 (MQ=255)                                                                           aTGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                            <  1:1143515/98‑1 (MQ=255)                                                                           aTGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                            >  2:1143515/1‑98 (MQ=255)                                                                            tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCttt                                                  >  1:34100/1‑139 (MQ=255)                                                                              ctgctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgcc                                                                                          >  2:707972/1‑97 (MQ=255)                                                                              ctgctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgcc                                                                                          <  1:707972/97‑1 (MQ=255)                                                                              ctgctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAg                                                       <  1:748568/132‑1 (MQ=255)                                                                              atgctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGTGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGTAGGATGCCAg                                                       >  2:748568/2‑132 (MQ=255)                                                                                 ctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGttga                                             >  2:1052684/1‑139 (MQ=255)                                                                                  tgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCt                                                                                >  2:1084872/1‑103 (MQ=255)                                                                                  tgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCt                                                                                <  1:1084872/103‑1 (MQ=255)                                                                                      tttCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGaa                                        >  2:314009/1‑139 (MQ=255)                                                                                           gTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTcc                                                                                  >  1:674951/1‑92 (MQ=255)                                                                                           gTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTcc                                                                                  <  2:674951/92‑1 (MQ=255)                                                                                                tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCa                                                                >  1:718785/1‑105 (MQ=255)                                                                                                tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCa                                                                <  2:718785/105‑1 (MQ=255)                                                                                                     ttATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATc                         <  2:1218217/139‑1 (MQ=255)                                                                                                          aGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgccg                                                                                         >  2:885479/1‑70 (MQ=255)                                                                                                          aGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgccg                                                                                         <  1:885479/70‑1 (MQ=255)                                                                                                            cAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTa                  <  1:1052684/139‑1 (MQ=255)                                                                                                             aaTTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                            <  2:1082991/64‑1 (MQ=255)                                                                                                             aaTTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                            >  1:1082991/1‑64 (MQ=255)                                                                                                               tttttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCg               >  2:834286/1‑139 (MQ=255)                                                                                                                  ttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTTGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTg            <  2:768728/139‑1 (MQ=255)                                                                                                                  ttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTg            <  1:314009/139‑1 (MQ=255)                                                                                                                         aGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCAt     >  2:67743/1‑139 (MQ=255)                                                                                                                            tgtAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGa  >  1:99598/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       CCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGA  >  NC_000913/2591476‑2591737

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