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

CCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAA  >  NC_000913/2591467‑2591729                                                                                                                                         |                                                                                                                               ccGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATc                                                                                                                              >  2:1014310/1‑139 (MQ=255)          ccTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                     <  1:125954/139‑1 (MQ=255)             ctGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGccc                                                                                                                  <  1:1014310/139‑1 (MQ=255)              tGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCg                                                                                                                 >  1:411634/1‑139 (MQ=255)                    ggcATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCg                                                                                                                       >  2:20978/1‑127 (MQ=255)                    ggcATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCg                                                                                                                       <  1:20978/127‑1 (MQ=255)                      cATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGa                                                                                                         <  2:411634/139‑1 (MQ=255)                        ttAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGc                                                                                                       >  2:1761045/1‑139 (MQ=255)                               cTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGAcaca                                                                                                <  1:1761045/139‑1 (MQ=255)                               cTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGAcaca                                                                                                >  1:1233866/1‑139 (MQ=255)                                tGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACacaa                                                                                               >  1:491017/1‑139 (MQ=255)                                            aaTTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGa                                                                                                         >  1:1690312/1‑117 (MQ=255)                                            aaTTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGa                                                                                                         <  2:1690312/117‑1 (MQ=255)                                              ttATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgccg                                                                                 <  2:1233866/139‑1 (MQ=255)                                                aTCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                               <  2:1454382/139‑1 (MQ=255)                                                aTCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                               >  1:1454382/1‑139 (MQ=255)                                                 tCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGccc                                                                              >  1:460651/1‑139 (MQ=255)                                                 tCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGccc                                                                              <  2:1825295/139‑1 (MQ=255)                                                   aaaCGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGa                                                                                                    >  2:745271/1‑115 (MQ=255)                                                   aaaCGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGa                                                                                                    <  1:745271/115‑1 (MQ=255)                                                      cGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGa                                                                                                    >  2:1584766/1‑112 (MQ=255)                                                      cGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGa                                                                                                    <  1:1584766/112‑1 (MQ=255)                                                              gCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCt                                                                 <  2:842884/139‑1 (MQ=255)                                                               ctctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                   >  1:129181/1‑88 (MQ=255)                                                               ctctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                   <  2:129181/88‑1 (MQ=255)                                                                tctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtga                                                               <  2:858614/139‑1 (MQ=255)                                                                tctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtga                                                               <  2:53864/139‑1 (MQ=255)                                                                    cGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtga                                                               <  2:1105801/135‑1 (MQ=255)                                                                    cGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtga                                                               >  1:1105801/1‑135 (MQ=255)                                                                         ccGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAgg                                                      >  1:846484/1‑139 (MQ=255)                                                                                     tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgcc                                                                                  <  2:460114/99‑1 (MQ=255)                                                                                     tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcgcc                                                                                  >  1:460114/1‑99 (MQ=255)                                                                                           tgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                     >  2:1130481/1‑58 (MQ=255)                                                                                           tgGGTTCCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTg                                                                                                                     <  1:1130481/58‑1 (MQ=255)                                                                                                    gTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTa                                                                                                             >  1:24150/1‑57 (MQ=255)                                                                                                    gTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTa                                                                                                             <  2:24150/57‑1 (MQ=255)                                                                                                       aTTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTcc                                                                          >  1:709500/1‑89 (MQ=255)                                                                                                       aTTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTcc                                                                          <  2:709500/89‑1 (MQ=255)                                                                                                         tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCt                                                                             >  1:1506599/1‑84 (MQ=255)                                                                                                         tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCt                                                                             <  2:1506599/84‑1 (MQ=255)                                                                                                         tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtgatga                                                            <  1:1483767/101‑1 (MQ=255)                                                                                                         tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtgatga                                                            >  2:1483767/1‑101 (MQ=255)                                                                                                             gTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGAt                  <  1:444218/139‑1 (MQ=255)                                                                                                                   aGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGttgattga                                 <  1:635745/118‑1 (MQ=255)                                                                                                                   aGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGttgattga                                 >  2:635745/1‑118 (MQ=255)                                                                                                                        tttttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGAcaca                                                                                                <  1:71902/50‑1 (MQ=255)                                                                                                                        tttttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGAcaca                                                                                                >  2:71902/1‑50 (MQ=255)                                                                                                                           ttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTg    <  2:1644720/139‑1 (MQ=255)                                                                                                                             ccATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGaa  <  2:491017/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                               CCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAA  >  NC_000913/2591467‑2591729

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: ‑