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%	196.7 / NA	55	L105S (TTA→TCA)	yrbG	putative calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/29);  total (26/29)

ACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGT  >  NC_000913/3340466‑3340718                                                                                                                           |                                                                                                                                   acaacGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGttgttg                                                                                                                    >  2:725056/1‑139 (MQ=255)               ccGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                      <  1:249289/139‑1 (MQ=255)                 gTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                    <  2:880206/139‑1 (MQ=255)                 gTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                    <  1:926063/139‑1 (MQ=255)                 gTTGGTACAGCCCTCGGCTCAAACATCATCAATATATAGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                    <  2:592380/139‑1 (MQ=255)                    ggTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                 <  1:903253/139‑1 (MQ=255)                    ggTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                 >  1:459162/1‑139 (MQ=255)                         aGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTa                                                                                            <  1:365825/139‑1 (MQ=255)                           cccTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                          <  2:889441/139‑1 (MQ=255)                               cGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                 <  1:399958/128‑1 (MQ=255)                               cGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                 >  2:399958/1‑128 (MQ=255)                                     aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGa                                                                                >  2:280091/1‑139 (MQ=255)                                                 tatatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGc                                                                    >  1:403186/1‑139 (MQ=255)                                                     ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCg                                                                >  2:682706/1‑139 (MQ=255)                                                      tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                          >  2:643159/1‑112 (MQ=255)                                                      tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                          <  1:643159/112‑1 (MQ=255)                                                      tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                               <  1:280091/139‑1 (MQ=255)                                                       gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCGCCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                              <  1:840636/107‑1 (MQ=255)                                                       gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTATTGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                              >  2:840636/1‑107 (MQ=255)                                                        cTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATg                                                             <  1:682706/139‑1 (MQ=255)                                                           aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                          >  2:843722/1‑139 (MQ=255)                                                                   tCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                  <  2:751168/139‑1 (MQ=255)                                                                     tGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCg                                                                                                             <  2:350380/78‑1 (MQ=255)                                                                     tGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCg                                                                                                             >  1:350380/1‑78 (MQ=255)                                                                                     tCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                     <  2:668673/86‑1 (MQ=255)                                                                                     tCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                     >  1:668673/1‑86 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCg                                                                                                             >  2:277152/1‑51 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCg                                                                                                             <  1:277152/51‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                      <  1:12296/58‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                      >  2:12296/1‑58 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGt                                                                       <  1:449715/89‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGt                                                                       >  2:449715/1‑89 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGc                                                                    >  1:483471/1‑92 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGc                                                                    <  2:483471/92‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                               <  2:220856/97‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                               >  1:220856/1‑97 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCtg                                     <  1:438537/123‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCtg                                     >  2:438537/1‑123 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                     >  1:932451/1‑139 (MQ=255)                                                                                                      cTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGt               <  1:42852/139‑1 (MQ=255)                                                                                                        gATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTa             <  1:725056/139‑1 (MQ=255)                                                                                                          tGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTaaa           >  1:227296/1‑139 (MQ=255)                                                                                                            ttCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACt         >  1:224737/1‑139 (MQ=255)                                                                                                            ttCTAAGCCGTAAATCACCCTTAATATTATTGGTCAGCGTCGTCGTCGGTTACGTACTCTATGACGGACACCTTAGTCGCATCGATGGTATCTTTCTCCTCTTTATGGCTTTCCTATGGCTGCTGTTCAtgctcaaacc         >  2:68861/1‑130 (MQ=255)                                                                                                               tACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                        <  2:432647/105‑1 (MQ=255)                                                                                                               tACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                        >  1:432647/1‑105 (MQ=255)                                                                                                                 cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTAt                               <  2:481674/112‑1 (MQ=255)                                                                                                                 cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTAt                               >  1:481674/1‑112 (MQ=255)                                                                                                                   ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                              <  2:275396/47‑1 (MQ=255)                                                                                                                   ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                              >  1:275396/1‑47 (MQ=255)                                                                                                                   ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                          <  2:405568/115‑1 (MQ=255)                                                                                                                   ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                          >  1:405568/1‑115 (MQ=255)                                                                                                                   ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGt  <  2:227296/139‑1 (MQ=255)                                                                                                                      tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGAc                                                                               >  2:286131/1‑59 (MQ=255)                                                                                                                      tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGAc                                                                               <  1:286131/59‑1 (MQ=255)                                                                                                                           |                                                                                                                                   ACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGT  >  NC_000913/3340466‑3340718

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