Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,340,588	T→C	L105S (TTA→TCA)	yrbG →	putative calcium/sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,340,588	0	T	C	100.0%	251.5 / NA	69	L105S (TTA→TCA)	yrbG	putative calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/34);  total (35/34)

CTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAG  >  NC_000913/3340458‑3340721                                                                                                                                   |                                                                                                                                      ctGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTaa                                                                                                                               <  1:431566/139‑1 (MQ=255) ctGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTaa                                                                                                                               <  1:170937/139‑1 (MQ=255)                       ccGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                         <  1:859480/139‑1 (MQ=255)                            ggTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCCTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                    <  1:2312/139‑1 (MQ=255)                                      tCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTa                                                                                          >  1:460297/1‑139 (MQ=255)                                      tCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTa                                                                                          >  1:223346/1‑139 (MQ=255)                                        ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                        >  2:486141/1‑139 (MQ=255)                                           tCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACg                                                                                     <  1:1011144/139‑1 (MQ=255)                                             aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                 >  2:418569/1‑125 (MQ=255)                                             aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                 <  1:418569/125‑1 (MQ=255)                                             aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGa                                                                                   >  1:1218620/1‑139 (MQ=255)                                                        atatatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCg                                                                        <  1:1306897/139‑1 (MQ=255)                                                         tatatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGc                                                                       <  2:153930/139‑1 (MQ=255)                                                          atatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCa                                                                      >  1:361420/1‑139 (MQ=255)                                                           tatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                     <  2:530381/139‑1 (MQ=255)                                                             ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCg                                                                   >  2:168820/1‑139 (MQ=255)                                                              tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                       <  1:1174907/102‑1 (MQ=255)                                                              tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                       >  2:1174907/1‑102 (MQ=255)                                                               gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGAt                                                                 >  1:735740/1‑139 (MQ=255)                                                                cTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                 <  2:533386/106‑1 (MQ=255)                                                                cTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                 >  1:533386/1‑106 (MQ=255)                                                                   aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                     >  1:579397/1‑131 (MQ=255)                                                                   aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                     <  2:579397/131‑1 (MQ=255)                                                                   aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                             >  2:25296/1‑139 (MQ=255)                                                                    tCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTa                                                                                                                                <  1:119098/71‑1 (MQ=255)                                                                    tCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTa                                                                                                                                >  2:119098/1‑71 (MQ=255)                                                                     ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                         <  2:52496/93‑1 (MQ=255)                                                                     ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                         >  1:52496/1‑93 (MQ=255)                                                                     ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACg                                                                                     <  1:375620/113‑1 (MQ=255)                                                                     ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACg                                                                                     >  2:375620/1‑113 (MQ=255)                                                                     ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATc                                                           >  2:125944/1‑139 (MQ=255)                                                                        cGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGcc                                                                                                        <  2:70500/91‑1 (MQ=255)                                                                        cGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGcc                                                                                                        >  1:70500/1‑91 (MQ=255)                                                                          gTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                         >  1:927124/1‑88 (MQ=255)                                                                          gTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTACTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                         <  2:927124/88‑1 (MQ=255)                                                                           tCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                     <  1:125944/139‑1 (MQ=255)                                                                               gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                         >  1:1311492/1‑83 (MQ=255)                                                                               gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                         <  2:1311492/83‑1 (MQ=255)                                                                               gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCctct                                                 >  2:921517/1‑139 (MQ=255)                                                                                     cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCt                                                                                                                                  >  1:1258425/1‑52 (MQ=255)                                                                                     cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCt                                                                                                                                  <  2:1258425/52‑1 (MQ=255)                                                                                     cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                     >  1:763425/1‑129 (MQ=255)                                                                                     cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                     <  2:763425/129‑1 (MQ=255)                                                                                     cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                           <  1:486141/139‑1 (MQ=255)                                                                                     cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                           <  2:1083234/139‑1 (MQ=255)                                                                                      tGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGc                                          <  1:921517/139‑1 (MQ=255)                                                                                             tCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTa                                   >  2:204975/1‑139 (MQ=255)                                                                                                   tACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                             <  2:460297/139‑1 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                    <  2:1253119/63‑1 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                    >  1:1253119/1‑63 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                        <  1:375673/75‑1 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                        >  2:375673/1‑75 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATgg                                                               >  1:1058371/1‑100 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATgg                                                               <  2:1058371/100‑1 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGc                                          <  1:1018640/121‑1 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGc                                          >  2:1018640/1‑121 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                             <  2:1205904/134‑1 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                             >  1:1205904/1‑134 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                        >  1:351609/1‑139 (MQ=255)                                                                                                        tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                        >  1:475261/1‑139 (MQ=255)                                                                                                          cATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCa                      >  1:1015099/1‑139 (MQ=255)                                                                                                                  tGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                        >  2:44723/1‑65 (MQ=255)                                                                                                                  tGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                        <  1:44723/65‑1 (MQ=255)                                                                                                                    ttCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACt            >  2:687988/1‑139 (MQ=255)                                                                                                                      cTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTg          <  2:223346/139‑1 (MQ=255)                                                                                                                      cTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTg          <  2:735740/139‑1 (MQ=255)                                                                                                                              tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCg                                                                   >  1:188354/1‑74 (MQ=255)                                                                                                                              tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCg                                                                   <  2:188354/74‑1 (MQ=255)                                                                                                                              tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAg  >  1:462302/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                      CTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAG  >  NC_000913/3340458‑3340721

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