Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,281,575	(A)5→4	coding (1596/1650 nt)	yjcE →	putative cation/proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,281,571	0	A	.	100.0%	196.5 / NA	41	coding (1592/1650 nt)	yjcE	putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/22);  total (19/22)

GCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTC  >  NC_000913/4281422‑4281715                                                                                                                                                      |                                                                                                                                                 gCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCa                                                                                                                                                  <  1:667730/150‑1 (MQ=255)  cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCaaa                                                                                                                                                <  1:175127/151‑1 (MQ=255)     cgtcgCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCaaaat                                                                                                                                              <  1:50987/150‑2 (MQ=255)                            cgTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTg                                                                                                                     <  1:491904/151‑1 (MQ=255)                             gTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTg                                                                                                                     >  1:632638/1‑150 (MQ=255)                                  aaGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGt                                                                                                               >  1:60430/1‑151 (MQ=255)                                   aGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGtt                                                                                                              >  2:53676/1‑151 (MQ=255)                                    gTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTg                                                                                                             <  2:337355/151‑1 (MQ=255)                                       ccGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCt                                                                                                           >  2:281599/1‑150 (MQ=255)                                       ccGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTa                                                                                                          <  2:376127/151‑1 (MQ=255)                                                 agagAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGaa                                                                                                 >  2:332389/1‑150 (MQ=255)                                                     aaCCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCa                                                                                            <  1:89754/151‑1 (MQ=255)                                                      aCCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAg                                                                                           <  1:423757/151‑1 (MQ=255)                                                           gAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAg                                                                                       <  1:364285/150‑1 (MQ=255)                                                             gCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCg                                                                                    <  2:91446/151‑1 (MQ=255)                                                               gTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGaa                                                                                  <  1:233012/151‑1 (MQ=255)                                                                  gCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGc                                                                               >  1:213347/1‑151 (MQ=255)                                                                    ttCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCg                                                                              >  2:341739/1‑150 (MQ=255)                                                                    ttCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCg                                                                              >  1:70515/1‑150 (MQ=255)                                                                       cGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACg                                                                           >  2:397479/1‑150 (MQ=255)                                                                                  aTTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCg                                                               >  1:61104/1‑151 (MQ=255)                                                                                    tGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGa                                                              <  2:647912/150‑1 (MQ=255)                                                                                     gCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAg                                                             <  2:67801/150‑1 (MQ=255)                                                                                      cGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGt                                                           >  2:33468/1‑151 (MQ=255)                                                                                                  gCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGaaaaa                                                <  1:426490/150‑1 (MQ=255)                                                                                                  gCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGaaaaa                                                <  2:255266/150‑1 (MQ=255)                                                                                                    tGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGt                                             >  1:273188/1‑151 (MQ=255)                                                                                                       aCTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCa                                          <  2:271686/151‑1 (MQ=255)                                                                                                          ttACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCaa                                         <  2:311385/149‑1 (MQ=255)                                                                                                              caccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                   >  2:613934/1‑151 (MQ=255)                                                                                                                   gcgcgcCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAcc                              >  1:177009/1‑151 (MQ=255)                                                                                                                     gcgcCACACGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                            <  2:177004/151‑1 (MQ=255)                                                                                                                          acacGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACa                       <  2:215275/151‑1 (MQ=255)                                                                                                                             cGGGAGATCAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGt                    >  1:156235/1‑151 (MQ=255)                                                                                                                                     cAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                           >  2:352741/1‑136 (MQ=255)                                                                                                                                     cAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                           <  1:352751/136‑1 (MQ=255)                                                                                                                                     cAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCgg            <  2:641223/151‑1 (MQ=255)                                                                                                                                     cAGCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCg             <  2:493284/150‑1 (MQ=255)                                                                                                                                      tgCAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCggg           >  2:176666/2‑151 (MQ=255)                                                                                                                                        cAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCt         <  1:133192/151‑1 (MQ=255)                                                                                                                                        cAACGAAACGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCt         >  1:565946/1‑151 (MQ=255)                                                                                                                                              aaCGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACAtt   <  2:306304/151‑1 (MQ=255)                                                                                                                                              aaCGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACAtt   >  1:158894/1‑151 (MQ=255)                                                                                                                                               aCGCTGC‑AAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTc  <  2:213342/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTC  >  NC_000913/4281422‑4281715

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