Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	341,200	A→G	T26A (ACC→GCC)	yahJ →	putative deaminase with metallo‑dependent hydrolase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	341,200	0	A	G	100.0%	182.3 / NA	51	T26A (ACC→GCC)	yahJ	putative deaminase with metallo‑dependent hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/27);  total (24/27)

GCCTGCCTGCGGCGCTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGAC  >  NC_000913/341063‑341322                                                                                                                                          |                                                                                                                           gcctgcctGCGGCGCTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGcc                                                                                                                          >  1:291760/1‑140 (MQ=255)  cctgcctgCGGCGCTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCct                                                                                                                         >  1:709684/1‑140 (MQ=255)   ctgcctgcGGCGCTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCctc                                                                                                                        >  1:1007021/1‑140 (MQ=255)            gcgcTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCt                                                                                                               <  1:1076508/140‑1 (MQ=255)            gcgcTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCt                                                                                                               <  1:233478/140‑1 (MQ=255)            gcgcTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCt                                                                                                               <  2:88293/140‑1 (MQ=255)            gcgcTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCt                                                                                                               >  2:794825/1‑140 (MQ=255)                tGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGcc                                                                                                           >  2:1069002/1‑140 (MQ=255)                  ggcggcgAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCa                                                                                                         <  2:574503/139‑1 (MQ=255)                    gggcgAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATg                                                                                                       <  2:225221/139‑1 (MQ=255)                      gcgAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCC‑CTCGCCCATgcgg                                                                                                    <  2:709684/140‑1 (MQ=255)                       cgAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATTAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATgcgg                                                                                                    >  2:264631/1‑140 (MQ=255)                                cACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTg                                                                                           <  1:1069002/140‑1 (MQ=255)                                 aCTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTgg                                                                                          <  1:780062/140‑1 (MQ=255)                                 aCTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTg                                                                                           <  2:817706/139‑1 (MQ=255)                                 aCTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTg                                                                                           >  1:817706/1‑139 (MQ=255)                                  cTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGc                                                                                         <  1:907586/140‑1 (MQ=255)                                     atatTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCAcc                                                                                      >  1:652130/1‑140 (MQ=255)                                             aaaaCGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTg                                                                              >  1:700196/1‑140 (MQ=255)                                             aaaaCGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTg                                                                              >  1:851893/1‑140 (MQ=255)                                               aaCGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTgcg                                                                            <  2:700196/140‑1 (MQ=255)                                               aaCGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTgcg                                                                            >  2:741216/1‑140 (MQ=255)                                               aaCGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTgcg                                                                            >  2:764084/1‑140 (MQ=255)                                                 cGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTgcgaa                                                                          <  1:741216/140‑1 (MQ=255)                                                   tAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTgg                                                                                          <  1:125482/122‑1 (MQ=255)                                                   tAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTgg                                                                                          >  2:125482/1‑122 (MQ=255)                                                      ggACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTgcgaagcgaa                                                                     <  2:794615/140‑1 (MQ=255)                                                            aaGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAacaaca                                                                >  2:742141/1‑139 (MQ=255)                                                            aaGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAacaaca                                                                <  1:742141/139‑1 (MQ=255)                                                            aaGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAacac                                                               >  2:875851/1‑140 (MQ=255)                                                                  aaagaaagCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAacaaca                                                                >  2:1064359/1‑133 (MQ=255)                                                                  aaagaaagCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAacaaca                                                                <  1:1064359/133‑1 (MQ=255)                                                                  aaagaaagCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGaa                                                         >  1:56772/1‑140 (MQ=255)                                                                        agCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCAc                                                   <  2:551983/140‑1 (MQ=255)                                                                        agCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCAc                                                   >  2:389019/1‑140 (MQ=255)                                                                         gCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACt                                                  >  1:242347/1‑140 (MQ=255)                                                                         gCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACt                                                  >  1:828467/1‑140 (MQ=255)                                                                            aTAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGAc                                               <  1:373395/140‑1 (MQ=255)                                                                            aTAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGAc                                               <  1:389019/140‑1 (MQ=255)                                                                               gccgccGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCg                                            <  1:268334/140‑1 (MQ=255)                                                                               gccgccGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCg                                            <  2:652130/140‑1 (MQ=255)                                                                                   ccgTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCAtt                                        <  1:103570/140‑1 (MQ=255)                                                                                    cgTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTa                                       >  2:1012488/1‑140 (MQ=255)                                                                                     gTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTAc                                      <  2:56772/140‑1 (MQ=255)                                                                                     gTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTAc                                      <  1:1012488/140‑1 (MQ=255)                                                                                      tGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACt                                     <  2:1007021/140‑1 (MQ=255)                                                                                        aaTTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCGTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCg                                            >  2:970290/1‑131 (MQ=255)                                                                                        aaTTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCg                                            <  1:970290/131‑1 (MQ=255)                                                                                        aaTTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCAGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTAt                                   >  1:699274/1‑140 (MQ=255)                                                                                                     gagCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGtgc                      >  2:324570/1‑140 (MQ=255)                                                                                                              gATGGTCACCGCCGCCGCGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGaaa             <  2:828467/140‑1 (MQ=255)                                                                                                                         ccgccgcGCTGTTTGGTGCCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGAc  >  2:272930/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                           GCCTGCCTGCGGCGCTGCGCGGCGAAACGGGCACTCATATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGAC  >  NC_000913/341063‑341322

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