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

TCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAG  >  NC_000913/3340455‑3340721                                                                                                                                      |                                                                                                                                        tctcTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCt                                                                                                                                    >  2:837825/1‑139 (MQ=255)           aacaacGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATgttgtt                                                                                                                          <  1:396775/139‑1 (MQ=255)                     tttAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                                <  1:738693/139‑1 (MQ=255)                            gTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           <  1:23536/137‑1 (MQ=255)                            gTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           >  2:23536/1‑137 (MQ=255)                                  aCAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   <  1:186348/139‑1 (MQ=255)                                        cTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                         >  1:772601/1‑127 (MQ=255)                                        cTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGACTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                         <  2:772601/127‑1 (MQ=255)                                          cGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCTGGTTCCGTA‑‑CTCTAt                                                                                         >  1:428091/1‑139 (MQ=255)                                           ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                         <  2:519083/124‑1 (MQ=255)                                           ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                         >  1:519083/1‑124 (MQ=255)                                                       aTCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   <  1:215758/118‑1 (MQ=255)                                                       aTCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   >  2:215758/1‑118 (MQ=255)                                                            tatatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           >  1:799963/1‑105 (MQ=255)                                                            tatatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           <  2:799963/105‑1 (MQ=255)                                                            tatatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGc                                                                       >  2:1044346/1‑139 (MQ=255)                                                              tatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGttgttg                                                                                                                         <  1:820777/89‑1 (MQ=255)                                                              tatTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGttgttg                                                                                                                         >  2:820777/1‑89 (MQ=255)                                                                ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                                <  2:766772/96‑1 (MQ=255)                                                                ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                                >  1:766772/1‑96 (MQ=255)                                                                ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTCTATGACGGACAACTTAGTCGCAg                                                                     >  1:817514/1‑139 (MQ=255)                                                                ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGc                                                                    <  2:30970/138‑1 (MQ=255)                                                                ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGc                                                                    >  1:30970/1‑138 (MQ=255)                                                                  gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAc                                                                                                                                       >  1:230317/1‑71 (MQ=255)                                                                      caccTCGATCTGGCTGCGCTGGTTCGATATTTTACCGCCCATTCTGATGATCTGTGCCGTGAATCACCCTTAAATTTGTTGGTCAGCGTGGTTGACGCTTCCGTA‑‑CTCTATGACAGGCAACTTAGTCGCAGCGTTGGTa                                                             >  1:111659/3‑139 (MQ=255)                                                                      aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTcc                                                                                                    >  2:921641/1‑102 (MQ=255)                                                                      aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTcc                                                                                                    <  1:921641/102‑1 (MQ=255)                                                                      aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                             >  1:212765/1‑139 (MQ=255)                                                                                tGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGa                                                                  >  1:663081/1‑124 (MQ=255)                                                                                tGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGa                                                                  <  2:663081/124‑1 (MQ=255)                                                                                  gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   >  1:804126/1‑91 (MQ=255)                                                                                  gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   <  2:804126/91‑1 (MQ=255)                                                                                      cgcTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAg                                                                     <  1:895543/115‑1 (MQ=255)                                                                                      cgcTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAg                                                                     >  2:895543/1‑115 (MQ=255)                                                                                        cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                           >  2:1020766/1‑139 (MQ=255)                                                                                        cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                           <  1:837825/139‑1 (MQ=255)                                                                                         tGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGc                                          <  1:368223/139‑1 (MQ=255)                                                                                                 ccTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTAt                                  <  1:390645/139‑1 (MQ=255)                                                                                                           tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                        >  1:1013536/1‑139 (MQ=255)                                                                                                           tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                        <  1:1020766/139‑1 (MQ=255)                                                                                                           tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGACAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGTCGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                        >  2:766493/1‑139 (MQ=255)                                                                                                             cATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGc                                     <  2:752410/124‑1 (MQ=255)                                                                                                             cATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGc                                     >  1:752410/1‑124 (MQ=255)                                                                                                                          tACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCTGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGc         <  2:428091/139‑1 (MQ=255)                                                                                                                              ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                        <  1:285262/120‑1 (MQ=255)                                                                                                                              ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                        >  2:285262/1‑120 (MQ=255)                                                                                                                                 tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAg                                                                     <  2:528003/72‑1 (MQ=255)                                                                                                                                 tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAg                                                                     >  1:528003/1‑72 (MQ=255)                                                                                                                                 tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                             <  1:1195508/112‑1 (MQ=255)                                                                                                                                 tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgctg                             >  2:1195508/1‑112 (MQ=255)                                                                                                                                 tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAg  >  1:638251/1‑139 (MQ=255)                                                                                                                                 tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAg  >  2:850934/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        TCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTA‑‑CTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAG  >  NC_000913/3340455‑3340721

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