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

ACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTC  >  NC_000913/3340466‑3340719                                                                                                                           |                                                                                                                                    acaacGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGttgttg                                                                                                                     >  1:151404/1‑139 (MQ=255)    acGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTc                                                                                                                  >  2:701023/1‑139 (MQ=255)    acGCGATTTAGCCGTTGGTACAGCCCTCGGCTC‑AACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCa                                                                                                                 <  2:770202/139‑1 (MQ=255)             aGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtggtg                                                                                                         <  1:345968/139‑1 (MQ=255)                    ggTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGtt                                                                                                  >  1:431661/1‑139 (MQ=255)                       aCAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               <  1:157471/139‑1 (MQ=255)                                     aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAc                                                                                                                                   >  1:30258/1‑89 (MQ=255)                                     aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGa                                                                                 <  1:227714/139‑1 (MQ=255)                                            aTCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                     <  2:520986/112‑1 (MQ=255)                                            aTCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                     >  1:520986/1‑112 (MQ=255)                                                   tatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                   <  2:551514/139‑1 (MQ=255)                                                       gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGAt                                                               >  1:691139/1‑139 (MQ=255)                                                           aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                      <  1:231479/112‑1 (MQ=255)                                                           aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                      >  2:231479/1‑112 (MQ=255)                                                           aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           >  2:189137/1‑139 (MQ=255)                                                           aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           <  2:431661/139‑1 (MQ=255)                                                           aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           >  2:490193/1‑139 (MQ=255)                                                             ccccGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATc                                                         >  2:516045/4‑139 (MQ=255)                                                             ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               <  1:482158/101‑1 (MQ=255)                                                             ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               >  2:482158/1‑101 (MQ=255)                                                                cGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                           >  1:561026/1‑102 (MQ=255)                                                                cGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                           <  2:561026/102‑1 (MQ=255)                                                                        cTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCtt                                              <  1:516045/139‑1 (MQ=255)                                                                         tGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTACTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCttt                                             <  2:49761/139‑1 (MQ=255)                                                                          gcgcTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTc                                            <  2:151404/139‑1 (MQ=255)                                                                             cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACg                                                                                   >  2:297287/1‑97 (MQ=255)                                                                             cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACg                                                                                   <  1:297287/97‑1 (MQ=255)                                                                             cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                         <  2:148318/139‑1 (MQ=255)                                                                                           tACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                           <  2:167102/139‑1 (MQ=255)                                                                                           tACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                           <  1:189137/139‑1 (MQ=255)                                                                                           tACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                           <  1:490193/139‑1 (MQ=255)                                                                                               gTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               <  2:193927/67‑1 (MQ=255)                                                                                               gTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               >  1:193927/1‑67 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           <  2:83054/102‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           >  1:83054/1‑102 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCtt                                                       <  1:372000/106‑1 (MQ=255)                                                                                                tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCtt                                                       >  2:372000/1‑106 (MQ=255)                                                                                                         aTGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTaa             <  2:691139/139‑1 (MQ=255)                                                                                                                 cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               >  1:761195/1‑49 (MQ=255)                                                                                                                 cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               <  2:761195/49‑1 (MQ=255)                                                                                                                 cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCAc     >  1:463823/1‑139 (MQ=255)                                                                                                                   ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCgg                                                                                                    <  1:376326/42‑1 (MQ=255)                                                                                                                   ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCgg                                                                                                    >  2:376326/1‑42 (MQ=255)                                                                                                                    cGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTc  >  2:707276/1‑139 (MQ=255)                                                                                                                      tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGc                                        >  1:203130/1‑99 (MQ=255)                                                                                                                      tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGc                                        <  2:203130/99‑1 (MQ=255)                                                                                                                           |                                                                                                                                    ACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTC  >  NC_000913/3340466‑3340719

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