Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,859,124	A→G	G331G (GGT→GGC)	yidK ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,859,124	0	A	G	100.0%	151.7 / NA	42	G331G (GGT→GGC)	yidK	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/14);  total (28/14)

TTTGCGCCCGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAAC  >  NC_000913/3858977‑3859268                                                                                                                                                    |                                                                                                                                                 tttGCGCCCGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGccc                                                                                                                                               >  1:599608/1‑151 (MQ=255)  ctgagcCCGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGCCGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCCCGCCGCAGAAGccc                                                                                                                                               >  2:81598/5‑150 (MQ=255)    gcgcCCGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCccacc                                                                                                                                            <  1:275090/151‑1 (MQ=255)       cccGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCAtt                                                                                                                                         <  2:202226/151‑1 (MQ=255)        ccGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTg                                                                                                                                        <  2:242452/151‑1 (MQ=255)        ccGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTg                                                                                                                                        >  1:588142/1‑151 (MQ=255)                              gCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGAc                                                                                                                  >  1:450734/1‑151 (MQ=255)                                   gCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCg                                                                                                             <  1:648150/151‑1 (MQ=255)                                     aTTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTc                                                                                                           >  1:552529/1‑151 (MQ=255)                                        cTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCggg                                                                                                        >  2:370797/1‑151 (MQ=255)                                            ttAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAgg                                                                                                    >  2:404796/1‑151 (MQ=255)                                                  atgCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGt                                                                                              <  2:551436/151‑1 (MQ=255)                                                     cGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCgg                                                                                           >  2:632598/1‑151 (MQ=255)                                                        cGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCtt                                                                                        >  2:667291/1‑151 (MQ=255)                                                            aaaTACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCgg                                                                                    >  2:384771/1‑151 (MQ=255)                                                                aCCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTaa                                                                                 >  2:286057/1‑150 (MQ=255)                                                                aCCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTaa                                                                                 <  1:203062/150‑1 (MQ=255)                                                                          taTAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCAACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCAATGTCGGCTTTCGGTAAATACTGATAc                                                                       >  2:655506/2‑150 (MQ=255)                                                                               cGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGa                                                                 <  1:234242/151‑1 (MQ=255)                                                                                             tCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAAccc                                                   >  2:51425/1‑151 (MQ=255)                                                                                              cAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCt                                                  >  2:443113/1‑151 (MQ=255)                                                                                                gAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACGCCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAAGCCTGATACAGATGAAACGCAATCAACCCTgg                                                <  2:225779/151‑1 (MQ=255)                                                                                                   aGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGgcag                                             <  2:397546/151‑1 (MQ=255)                                                                                                    gCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGgcagc                                            <  1:197988/151‑1 (MQ=255)                                                                                                     ccGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGgcagc                                            >  1:286188/1‑150 (MQ=255)                                                                                                      cGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAc                                          <  2:462396/151‑1 (MQ=255)                                                                                                        ttGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAcc                                         >  2:457900/1‑150 (MQ=255)                                                                                                          gAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGt                                      >  1:344441/1‑151 (MQ=255)                                                                                                              gTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTAcc                                   >  1:464617/1‑150 (MQ=255)                                                                                                                aCTGATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGc                                >  1:140500/1‑151 (MQ=255)                                                                                                                 cTGATCACCGCACCAAATACCTCCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCTTTAAATCCTTATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAg                                 <  1:41682/149‑1 (MQ=255)                                                                                                                   gATCACCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGGAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCggg                             >  2:340445/1‑151 (MQ=255)                                                                                                                       aCCGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCg                          >  1:637746/1‑150 (MQ=255)                                                                                                                        ccGCACCAAATAACGCCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCg                          >  1:232620/1‑149 (MQ=255)                                                                                                                        ccGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCg                          >  1:148091/1‑149 (MQ=255)                                                                                                                         cGCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGc                       >  1:301774/1‑151 (MQ=255)                                                                                                                          gCACCAAATAACACCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCa                      <  2:95473/151‑1 (MQ=255)                                                                                                                                    acacCGCGCCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCAc            <  1:174534/151‑1 (MQ=255)                                                                                                                                         gcgcCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCAccgccg       >  2:366244/1‑151 (MQ=255)                                                                                                                                          cgcCGAAGAAGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGt      >  1:443815/1‑151 (MQ=255)                                                                                                                                              gaagaaGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAAc  >  2:249898/1‑151 (MQ=255)                                                                                                                                              gaagaaGCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAAc  >  1:595921/1‑151 (MQ=255)                                                                                                                                               aagaagCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAAc  >  1:264796/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                 TTTGCGCCCGACGGTGACCAGCTGCTGCGGCTCGGCATTCTGGTTAATGATGCGACGGTAAATACCCATACTGAATAACGTACTGGCGCTATTCAGAAAGCCGTTGAAGGTACTGATCACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAAC  >  NC_000913/3858977‑3859268

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