Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,707,320	A→T	9.4%	intergenic (‑196/+2)	lepA ← / ← rseC	back‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,707,320	0	A	T	9.4%	170.9 / 7.2	53	intergenic (‑196/+2)	lepA/rseC	back‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base T (2/3);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAG  >  NC_000913/2707175‑2707468                                                                                                                                                  |                                                                                                                                                     cGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATc                                                                                                                                                   <  2:207249/149‑1 (MQ=255)   cTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCa                                                                                                                                                  <  3:102698/148‑1 (MQ=255)     gTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTg                                                                                                                                               <  2:151108/149‑1 (MQ=255)      ttCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTgg                                                                                                                                              <  1:397155/149‑1 (MQ=255)         tGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtt                                                                                                                                           <  1:267714/149‑2 (MQ=255)             cTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgt                                                                                                                                        <  3:78198/148‑5 (MQ=255)             cTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGt                                                                                                                                       >  2:315891/1‑149 (MQ=255)              ttAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTc                                                                                                                                      >  1:300045/1‑149 (MQ=255)                 aCGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGTGAGGAAATGTTAATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATCCAAAATGCAACAAATCCAGTGAAATCACTGGCTCGAGTTTTc                                                                                                                                   >  2:380971/1‑149 (MQ=255)                     cTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGa                                                                                                                                <  4:4581/148‑1 (MQ=255)                       gTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGa                                                                                                                             <  2:354950/149‑1 (MQ=255)                            aaTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGtt                                                                                                                         <  3:287945/148‑1 (MQ=255)                               aGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACATCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCa                                                                                                                      >  3:243287/1‑148 (MQ=255)                                  aaCGCCGCATTCTACAGTACAACATTGAGGCGAGGAAATGGTAATACCATATGGATTGTGGTATCTGGGAACGTCCTCGCATTTGTTATGCGAAATGCAACAAAGCCAGTGAACTCAGTGGCTCGCGTCTTCCGAAGATGTTTCAAAt                                                                                                                   >  4:250149/1‑148 (MQ=255)                                      ccgcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGT‑GTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGcaca                                                                                                             >  1:35111/1‑149 (MQ=255)                                         cATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAg                                                                                                           <  2:300045/149‑1 (MQ=255)                                         cATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAg                                                                                                           <  2:5173/149‑1 (MQ=255)                                                   tacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGaa                                                                                                                               <  2:380377/119‑1 (MQ=255)                                                   tacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGaa                                                                                                                               >  1:380377/1‑119 (MQ=255)                                                     caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCa                                                                                                <  3:242936/148‑1 (MQ=255)                                                          ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTc                                                                                                                                   <  4:63689/108‑1 (MQ=255)                                                          ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTc                                                                                                                                   >  3:63689/1‑108 (MQ=255)                                                              ggCGAGGAAAGGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGc                                                                                       >  3:266447/1‑148 (MQ=255)                                                                cGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCtt                                                                                     >  4:80152/1‑148 (MQ=255)                                                                 gAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCTCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTa                                                                                    <  4:165336/148‑1 (MQ=255)                                                                   ggAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAg                                                                                  >  3:82614/1‑148 (MQ=255)                                                                    gAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                                                                >  1:350818/1‑149 (MQ=255)                                                                    gAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                                                                <  1:315891/149‑1 (MQ=255)                                                                    gAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAga                                                                                 <  4:52782/148‑1 (MQ=255)                                                                    gAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCACACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                                                                >  1:79713/1‑149 (MQ=255)                                                                           tCATACCGTATGGATTGT‑GTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGtt                                                                        <  2:35111/149‑1 (MQ=255)                                                                              tACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGc                                                                      >  2:135069/1‑149 (MQ=255)                                                                                     tGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTg                                                               >  1:123043/1‑149 (MQ=255)                                                                                            tgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCggg                                                         <  4:243287/148‑1 (MQ=255)                                                                                              gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACaa                                                                                                            <  2:172710/95‑1 (MQ=255)                                                                                              gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACaa                                                                                                            >  1:172710/1‑95 (MQ=255)                                                                                                aTCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcagca                                                    >  1:314208/1‑149 (MQ=255)                                                                                                    ggAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaa                                                                                                                                                      >  4:354975/1‑47 (MQ=255)                                                                                                    ggAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaa                                                                                                                                                      <  3:354975/47‑1 (MQ=255)                                                                                                      aaaCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGt                                                                         <  4:3163/122‑1 (MQ=255)                                                                                                      aaaCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGt                                                                         >  3:3163/1‑122 (MQ=255)                                                                                                      aaaCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTa                                              <  1:159028/149‑1 (MQ=255)                                                                                                         cGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcgc                                           <  1:135069/149‑1 (MQ=255)                                                                                                           tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcgcg                                          >  3:194744/1‑148 (MQ=255)                                                                                                           tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGa                                         <  1:322475/149‑1 (MQ=255)                                                                                                             cTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAgg                                                                                                          <  2:358319/82‑1 (MQ=255)                                                                                                             cTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAgg                                                                                                          >  1:358319/1‑82 (MQ=255)                                                                                                               cGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTa                                      >  4:114941/1‑148 (MQ=255)                                                                                                                            cAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCa                        <  2:356607/149‑1 (MQ=255)                                                                                                                               aaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAgg                      <  3:62984/148‑1 (MQ=255)                                                                                                                                        gaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCaa            >  2:262731/2‑149 (MQ=255)                                                                                                                                         aGCTAGTGATTTCACTGGCTCGCGTCTTCTGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCaa            >  4:340570/1‑148 (MQ=255)                                                                                                                                          gCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAAt           >  3:398809/1‑148 (MQ=255)                                                                                                                                            cAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCAtt                                                                  <  3:98332/91‑1 (MQ=255)                                                                                                                                            cAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCAtt                                                                  >  4:98332/1‑91 (MQ=255)                                                                                                                                                 aaaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGa   <  1:380971/149‑1 (MQ=255)                                                                                                                                                  aaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGa   <  3:142064/148‑1 (MQ=255)                                                                                                                                                  aaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAg  >  1:261154/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAG  >  NC_000913/2707175‑2707468

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