Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	473,618	A→G	H218R (CAC→CGC)	amtB →	ammonia/ammonium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	473,618	0	A	G	100.0%	87.2 / NA	26	H218R (CAC→CGC)	amtB	ammonia/ammonium transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/16);  total (10/16)

TGGTTTGCTGGCTTCTCACGGTGCGCTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCACAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACG  >  NC_000913/473481‑473745                                                                                                                                          |                                                                                                                                tggtTTGCTGGCTTCTCACGGTGCGCTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGcgca                                                                                                                               >  2:631142/1‑140 (MQ=255) tggtTTGCTGGCTTCTCACGGTGCGCTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGcgca                                                                                                                               >  1:681661/1‑140 (MQ=255)          ggCTTCTCACGGTGCGCTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGa                                                                                                                      <  1:22464/140‑1 (MQ=255)                   cGGTGCGCTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCa                                                                                                             <  1:603675/140‑1 (MQ=255)                         gcTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGa                                                                                                       <  1:595155/140‑1 (MQ=255)                              aTTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGcc                                                                                                  <  1:631142/140‑1 (MQ=255)                               tttCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCa                                                                                                 <  2:188103/140‑1 (MQ=255)                                        tGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCtata                                                                                        >  1:597527/1‑140 (MQ=255)                                                 ggtgCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTctct                                                                                           <  2:698268/128‑1 (MQ=255)                                                 ggtgCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTctct                                                                                           >  1:698268/1‑128 (MQ=255)                                                        aTTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCtt                                                                        >  1:757523/1‑140 (MQ=255)                                                            aCGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAAc                                                                    >  2:736152/1‑140 (MQ=255)                                                             cgccgcAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACg                                                                   >  1:586922/1‑140 (MQ=255)                                                                                   ggCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTCTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGaa                                             >  2:57296/1‑140 (MQ=255)                                                                                     cgcgTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACG‑CGGGTCAGCGGGCACGGCGAATGa                                          >  1:800450/1‑140 (MQ=255)                                                                                         gttCTGATGGGAAAACAAGTGGGCTTCGGTATAGAGGCGTTTAAACCGTGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAAt                                       <  2:404650/138‑1 (MQ=255)                                                                                           tCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATcg                                     <  1:57296/140‑1 (MQ=255)                                                                                           tCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATcg                                     <  2:83221/140‑1 (MQ=255)                                                                                                   gAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTg                             >  2:158813/1‑140 (MQ=255)                                                                                                    aaaaCGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCGGGGTCAGGGGGAACGGCGAATGAAATCGCGGCACTgg                            >  1:317691/1‑140 (MQ=255)                                                                                                      aaCGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGTGTCAGCGGGCAAGGAGAATGAGATCGCGGCACTGGCa                          >  1:242087/1‑140 (MQ=255)                                                                                                      aaCGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCa                          <  1:158813/140‑1 (MQ=255)                                                                                                      aaCGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCa                          <  1:736152/140‑1 (MQ=255)                                                                                                      aaCGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCa                          <  2:597527/140‑1 (MQ=255)                                                                                                      aaCGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACAGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCa                          <  2:485527/140‑1 (MQ=255)                                                                                                      aaCGCGTGGGCTCCGGTTAAGAAGCGTTTAACCCGCGCACCCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCa                          <  2:317691/140‑1 (MQ=255)                                                                                                        cgcgTGGGCTTCGGTAAAGAGGCGTTTAAACCGCGCAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCAtt                        <  2:757523/140‑1 (MQ=255)                                                                                                                              cGTTTAAACCGCGCAACTTGCCCATGGTCTTCACCTGGAATGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAAcg  <  2:242087/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                TGGTTTGCTGGCTTCTCACGGTGCGCTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCCGGTCTGGTGGGCGCGTATCTGATAGGAAAACGCGTGGGCTTCGGTAAAGAGGCGTTTAAACCGCACAACCTGCCGATGGTCTTCACCGGGACTGCCATTCTCTATATCGGTTGGTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACG  >  NC_000913/473481‑473745

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