Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,445,152	G→A	R137C (CGC→TGC)	yfcA ←	TauE/TSUP family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,445,152	0	G	A	100.0%	177.1 / NA	48	R137C (CGC→TGC)	yfcA	TauE/TSUP family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/26);  total (22/26)

GCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTCATGG  >  NC_000913/2445004‑2445298                                                                                                                                                     |                                                                                                                                                   gCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACt                                                                                                                                                  <  2:16195/151‑1 (MQ=255) gCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGCAAACCGTACATCCGGCACt                                                                                                                                                  >  1:7164/1‑151 (MQ=255)                ttGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCa                                                                                                                                    >  1:154956/1‑150 (MQ=255)                    cGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCAcc                                                                                                                                  >  2:393927/1‑148 (MQ=255)                       ggTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCtt                                                                                                                            >  1:354192/1‑151 (MQ=255)                        gTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCttt                                                                                                                           <  1:336368/151‑1 (MQ=255)                          tAAACCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCGCCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTg                                                                                                                          <  1:660491/150‑1 (MQ=255)                               ccGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCa                                                                                                                    <  2:256770/151‑1 (MQ=255)                                cGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGAAGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGc                                                                                                                        >  1:524552/1‑146 (MQ=255)                                      gCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGAACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGa                                                                                                              <  2:752940/150‑1 (MQ=255)                                        gTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAAt                                                                                                           <  1:683635/151‑1 (MQ=255)                                        gTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAAt                                                                                                           <  1:678419/151‑1 (MQ=255)                                          aaCGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATag                                                                                                         <  1:378181/151‑1 (MQ=255)                                           aCGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAAt                                                                                                           >  2:733195/1‑148 (MQ=255)                                                      cgcgTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACa                                                                                              >  2:743415/1‑150 (MQ=255)                                                       gcgTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTTGCATCAGCAAGAAATAGAGGGCGATACaa                                                                                             >  1:533662/1‑150 (MQ=255)                                                           aaaaCGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCCTCAGCAAGAAATAGAGGCCGATACaaaat                                                                                          >  1:601328/1‑147 (MQ=255)                                                             aaCGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTAcc                                                                                       >  1:364502/1‑150 (MQ=255)                                                              aCGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAg                                                                                     <  1:139923/151‑1 (MQ=255)                                                               cGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGa                                                                                    <  1:644241/151‑1 (MQ=255)                                                                    cGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATggg                                                                               <  2:759728/151‑1 (MQ=255)                                                                      gCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCa                                                                             <  2:67486/151‑1 (MQ=255)                                                                      gCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCa                                                                             <  2:592435/151‑1 (MQ=255)                                                                        aGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTTACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAAAGGGCaa                                                                            >  2:141755/1‑150 (MQ=255)                                                                         gCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCaaa                                                                           <  1:107123/150‑1 (MQ=255)                                                                            gggCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCt                                                                       >  2:102187/1‑151 (MQ=255)                                                                              gCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTTTTAACCTAGCTTTTGCATCAGCAAGAAATAGAGGGCGAAACAAAATACCAGAAATGGGaaaa                                                                          >  1:663160/1‑146 (MQ=255)                                                                                 aaaGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGt                                                                  <  2:580919/151‑1 (MQ=255)                                                                                         catcatAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACa                                                            >  1:673148/1‑149 (MQ=255)                                                                                         catcatAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATc                                                          <  2:14662/151‑1 (MQ=255)                                                                                             atAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGcc                                                      <  2:154989/151‑1 (MQ=255)                                                                                                aaaaCCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGa                                                   <  2:85978/151‑1 (MQ=255)                                                                                                 aaaCCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGaa                                                  <  1:418500/151‑1 (MQ=255)                                                                                                     cGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGt                                               >  2:512982/1‑150 (MQ=255)                                                                                                              cACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCAcc                                       >  2:468611/1‑149 (MQ=255)                                                                                                                  ggCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACCAATTACCAGAATGGGCAAAATCTGCCGGAAGACATCAGCCTGAAAGTATTGCACcagca                                   >  1:419936/1‑149 (MQ=255)                                                                                                                   gCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACcagca                                   >  2:554160/1‑148 (MQ=255)                                                                                                                     tATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATAAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTgcg                              >  1:278072/1‑151 (MQ=255)                                                                                                                      aTTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTATCTTTGGCATCAGCAAGAAATAGAGGGCGATACAAATTTCCTGAATGGGCAAAAACTGCCGTAAGACATCAGCCTGAACGTATTGCACcagcact                                 >  1:357611/1‑145 (MQ=255)                                                                                                                         aaCGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTTAACGCATTGCACCAGCAGTGcgcc                            >  1:196952/1‑149 (MQ=255)                                                                                                                          aCGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCt                         <  2:74122/151‑1 (MQ=255)                                                                                                                           cgcgAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCt                         <  1:601931/150‑1 (MQ=255)                                                                                                                                       aCCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACa            <  2:364549/151‑1 (MQ=255)                                                                                                                                         cGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACaaa          <  1:171/151‑1 (MQ=255)                                                                                                                                           tACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAgg        >  2:137284/1‑151 (MQ=255)                                                                                                                                            aCATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGt       <  1:699549/151‑1 (MQ=255)                                                                                                                                             cATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTc      <  2:630699/151‑1 (MQ=255)                                                                                                                                                 cGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTCATgg  <  2:628930/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTCATGG  >  NC_000913/2445004‑2445298

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