Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,258,991	A→G	E67E (GAA→GAG)	ychH →	DUF2583 family putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,258,991	0	A	G	100.0%	206.4 / NA	55	E67E (GAA→GAG)	ychH	DUF2583 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/27);  total (28/27)

GGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGT  >  NC_000913/1258856‑1259121                                                                                                                                        |                                                                                                                                   ggtCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAg                                                                                                                                 >  2:1234243/1‑139 (MQ=255)     ggCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGgtgt                                                                                                                             <  2:211691/139‑1 (MQ=255)     ggCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGgtgt                                                                                                                             <  2:61190/139‑1 (MQ=255)          ggggTATTCAATCCTCAACCAGTTACCATAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGAc                                                                                                                        <  2:789594/139‑1 (MQ=255)            ggTATTCAATCCTCAACCAGTTACCACAGTTTAATATACCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCg                                                                                                                      >  2:459829/1‑139 (MQ=255)                 tCAATCCTCAACCAGTTTCCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTAGTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACt                                                                                                                 >  2:822356/1‑139 (MQ=255)                     tCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACtggtg                                                                                                             >  1:623540/1‑139 (MQ=255)                            ccAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGc                                                                                                      <  2:1271512/139‑1 (MQ=255)                                      caGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACa                                                                                            <  1:459829/139‑1 (MQ=255)                                         tttAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAAc                                                                                         >  2:1179814/1‑139 (MQ=255)                                                 gCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCTCCACTATGACAAACGTTcccgc                                                                                 >  2:1135496/1‑139 (MQ=255)                                                 gCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTgccgc                                                                                 >  2:354563/1‑139 (MQ=255)                                                   cccAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGccgccg                                                                               >  1:475450/1‑139 (MQ=255)                                                     cAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACtggt                                                                                                              <  1:345424/106‑1 (MQ=255)                                                     cAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACtggt                                                                                                              >  2:345424/1‑106 (MQ=255)                                                         aTTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGa                                                                         <  1:1062882/139‑1 (MQ=255)                                                             cGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGTCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGccgccg                                                                               >  2:496778/1‑129 (MQ=255)                                                             cGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGTCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGccgccg                                                                               <  1:496778/129‑1 (MQ=255)                                                             cGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTgccgcc                                                                                >  2:350765/1‑128 (MQ=255)                                                             cGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTgccgcc                                                                                <  1:350765/128‑1 (MQ=255)                                                             cGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAAt                                                                     >  1:684014/1‑139 (MQ=255)                                                                      tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCAt                                                            <  2:684014/139‑1 (MQ=255)                                                                      tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCAt                                                            <  1:235287/139‑1 (MQ=255)                                                                      tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCACGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCAt                                                            <  2:296281/139‑1 (MQ=255)                                                                           tGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                       <  1:1268066/139‑1 (MQ=255)                                                                           tGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                       <  1:1234243/139‑1 (MQ=255)                                                                             cTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAAc                                                     >  1:1220313/1‑139 (MQ=255)                                                                              tAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACa                                                    >  2:61235/1‑139 (MQ=255)                                                                                aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTgc                                                                                    >  2:1194923/1‑105 (MQ=255)                                                                                aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTgc                                                                                    <  1:1194923/105‑1 (MQ=255)                                                                                aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACaaa                                                  >  1:718456/1‑139 (MQ=255)                                                                                aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACaaa                                                  >  1:51161/1‑139 (MQ=255)                                                                                      ttCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCTCCACTATGACAAACGTTCCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGc                                            <  1:1135496/139‑1 (MQ=255)                                                                                          tcgGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAg                                                          <  1:88547/121‑1 (MQ=255)                                                                                          tcgGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAg                                                          >  2:88547/1‑121 (MQ=255)                                                                                            ggTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGATCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTg                                      >  1:709686/1‑139 (MQ=255)                                                                                             gTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGc                                     <  2:652334/139‑1 (MQ=255)                                                                                              tGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                   <  2:109521/108‑1 (MQ=255)                                                                                              tGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                   >  1:109521/1‑108 (MQ=255)                                                                                                   ttCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTc                                                                                                        <  2:96021/66‑1 (MQ=255)                                                                                                   ttCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTc                                                                                                        >  1:96021/1‑66 (MQ=255)                                                                                                         ggctggcGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACcagtcagt                         <  1:324926/139‑1 (MQ=255)                                                                                                          gctggcGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGtt                        >  2:1281737/1‑139 (MQ=255)                                                                                                             ggcGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTggg                                                                                                           >  2:1134679/1‑53 (MQ=255)                                                                                                             ggcGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTggg                                                                                                           <  1:1134679/53‑1 (MQ=255)                                                                                                             ggcGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                   <  1:8679/93‑1 (MQ=255)                                                                                                             ggcGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                   >  2:8679/1‑93 (MQ=255)                                                                                                             ggcGGGTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTgc                     >  2:902910/1‑139 (MQ=255)                                                                                                                  gTGCCCGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                <  1:902910/139‑1 (MQ=255)                                                                                                                     cccGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGc                                                                                                      >  2:48017/1‑50 (MQ=255)                                                                                                                     cccGTGTTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGc                                                                                                      <  1:48017/50‑1 (MQ=255)                                                                                                                          gtTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATc                                                                    >  1:67860/1‑79 (MQ=255)                                                                                                                          gtTGGCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATc                                                                    <  2:67860/79‑1 (MQ=255)                                                                                                                             ggCGGGCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAAtct     >  1:302039/1‑139 (MQ=255)                                                                                                                                gggCATGAGCAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGt  <  1:1179814/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   GGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGT  >  NC_000913/1258856‑1259121

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