Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,802	A→G	100%	V222A (GTA→GCA)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	953,802	0	A	G	98.8%	272.2 / ‑4.5	81	V222A (GTA→GCA)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (41/39);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTT  >  NC_000913/953705‑953895                                                                                                  |                                                                                              gTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGAGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCa                                                                                             >  1:686587/1‑100 (MQ=255)   tATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAt                                                                                           >  2:2299246/1‑100 (MQ=255)     tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                         >  1:817288/1‑100 (MQ=255)     tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                         <  2:202170/100‑1 (MQ=255)      ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                        >  2:87099/3‑100 (MQ=255)        gTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCa                                                                                      <  2:1735180/100‑1 (MQ=255)         tgatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTACAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCAt                                                                                     >  1:725724/1‑100 (MQ=255)          gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGACCCCACTGCGGACCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCAt                                                                                     >  2:1244713/1‑99 (MQ=255)          gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCAt                                                                                     >  1:2243030/1‑99 (MQ=255)          gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATc                                                                                    >  2:1334535/1‑100 (MQ=255)           atgaAGTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACAGCGGTCCAGAACTCCGGAGATGCGAAGTCGCGGATTACAATACCCATCg                                                                                   <  2:2588193/100‑1 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACTGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                  <  2:719739/100‑1 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                  >  1:977706/1‑100 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                  >  2:2529745/1‑100 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                  <  1:1521672/100‑1 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                  >  2:805009/1‑100 (MQ=255)             gaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCggg                                                                                 >  1:476879/1‑100 (MQ=255)               aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                               >  1:685236/1‑100 (MQ=255)                 ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                             <  1:896499/100‑1 (MQ=255)                  tcgtcaCGGTCAGGTGAGAAAAATTTTCCCGTGCAGAACCGACTGCGGTCCAAATTTCCTGGGATGCGATGTCGCTCATTGCAATACCCATGAGGATc                                                                              <  2:1794050/95‑1 (MQ=255)                   catcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATa                                                                           <  1:670371/100‑1 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                          >  1:2332794/1‑100 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                          <  1:1662154/100‑1 (MQ=255)                     tcaTGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaaa                                                                         <  1:1334875/100‑1 (MQ=255)                     tcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaaa                                                                         >  1:560945/1‑100 (MQ=255)                     tcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaaa                                                                         >  1:578111/1‑100 (MQ=255)                      caCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                        <  1:1654496/100‑1 (MQ=255)                      caCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                        <  2:284615/100‑1 (MQ=255)                       aCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACa                                                                       >  2:435352/1‑100 (MQ=255)                           tCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGtt                                                                   >  1:463776/1‑100 (MQ=255)                              ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                >  1:1679640/1‑100 (MQ=255)                              ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                <  2:977696/100‑1 (MQ=255)                               gTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCg                                                               <  2:2553791/100‑1 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                           >  1:2070757/1‑100 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                           <  1:1993333/100‑1 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATAATAAACATGTTTGCGATAc                                                           >  1:1146115/1‑100 (MQ=255)                                     aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                         <  2:2121837/100‑1 (MQ=255)                                      naaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtgt                                                        >  1:119030/2‑100 (MQ=255)                                        aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGTATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                      <  1:2109620/100‑1 (MQ=255)                                        aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                      <  2:2243003/100‑1 (MQ=255)                                        aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                      <  1:239946/100‑1 (MQ=255)                                        aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                      <  2:1904524/100‑1 (MQ=255)                                          tttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTc                                                    <  1:668096/100‑1 (MQ=255)                                             ccGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                                                 >  1:1739738/1‑100 (MQ=255)                                              cGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaaa                                                >  2:1968549/1‑100 (MQ=255)                                              cGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATATTTGCGATACTGTGCTCaaaa                                                >  1:1688034/1‑100 (MQ=255)                                              cGGTGCAGAACCGACTGCGGTCCAAAAATCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                                                 >  2:815822/1‑99 (MQ=255)                                               ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                               <  1:153979/100‑1 (MQ=255)                                               ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                               <  2:682309/100‑1 (MQ=255)                                                  gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                            >  2:2465932/1‑100 (MQ=255)                                                  gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                            <  2:2040349/100‑1 (MQ=255)                                                   cAGAACCGACTGCGGTCCCAAATTCCGTGGATCCGAAGTCTCGGATTGCAATACCCATCGAGATCATAAACATGTTAGCGATACTGTGCTCAAAACCGCt                                           >  1:884680/1‑100 (MQ=255)                                                   cAGAACCGACTGCGGTCCAAAATTCGGGGGATCCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                               <  2:904306/96‑1 (MQ=255)                                                    aGAACCGACTGCAGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                          <  1:228563/100‑1 (MQ=255)                                                       aCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca                                       <  1:2299274/100‑1 (MQ=255)                                                       aCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGCca                                       >  2:766166/1‑100 (MQ=255)                                                        ccGCCTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                        <  2:2145941/98‑1 (MQ=255)                                                         cGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGTATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac                                     <  2:317827/100‑1 (MQ=255)                                                          gACTGCGGTCCAAAATTCCGGGGATGCAACGTCGCGGATTGCAATACGTATCGGGATCATAAACATGTTTGTGATACTGTGCTCAAAACCGCTGGcaaca                                    <  2:1188824/100‑1 (MQ=255)                                                             tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                 <  2:476875/100‑1 (MQ=255)                                                               cGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACAt                               >  2:1234476/1‑100 (MQ=255)                                                                ggTCCAAAATTCCGGGGATGCGAAGAAGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATc                              >  2:213854/1‑100 (MQ=255)                                                                 gTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcg                             >  2:1293212/1‑100 (MQ=255)                                                                 gTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcg                             >  1:2067189/1‑100 (MQ=255)                                                                  gccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgc                            <  1:1244725/99‑1 (MQ=255)                                                                   ccAAAATTCCGGGGATGCAAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                           >  2:1156177/1‑100 (MQ=255)                                                                    cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                          >  2:1807992/1‑100 (MQ=255)                                                                    cAAAATTCCGGGGATGCGAAGTCGCAGATTGCAAAACCCATCGGGATCATAAACATGTTTGCGATACTTTGCTCAAAACCGCTGGCAACAAACATCGCGa                          >  2:693208/1‑100 (MQ=255)                                                                      aaaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGAcc                        >  2:2525195/1‑100 (MQ=255)                                                                             ggggATGCGAAGTCGCGGATTGCAATACCCAGCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagca                 <  2:722947/100‑1 (MQ=255)                                                                                   gCGAAGTCTCGTATTGTAATACTCATCTGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa           <  1:749989/100‑1 (MQ=255)                                                                                   gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa           <  2:2012881/100‑1 (MQ=255)                                                                                   gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa           >  2:1028458/1‑100 (MQ=255)                                                                                    cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt          <  1:1293225/100‑1 (MQ=255)                                                                                       aGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGan       <  2:973280/100‑2 (MQ=255)                                                                                       aGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGaa       <  1:2501589/100‑1 (MQ=255)                                                                                         tCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACg     <  2:2304015/100‑1 (MQ=255)                                                                                          cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGc    >  1:1881515/1‑100 (MQ=255)                                                                                          cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGc    >  2:1265878/1‑100 (MQ=255)                                                                                           gcgGATTGCAATACACAACGGAATAATAAACATGTTTGCGATACTTTGCGCAAAACCGCGGTCAACAAACATCGCGACTTGCAGCACCATAATGAAGGCt   >  2:486373/1‑100 (MQ=255)                                                                                            cgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCtt  <  1:600887/100‑1 (MQ=255)                                                                                                  |                                                                                              GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTT  >  NC_000913/953705‑953895

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑