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.6%	277.9 / ‑4.5	83	V222A (GTA→GCA)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (40/41);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCT  >  NC_000913/953706‑953894                                                                                                 |                                                                                             ctgaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCaa                                                                                           >  1:98500/7‑100 (MQ=255)  tgaccgTGATGAATTTCATCACGGTCAAGTGAGAAAAATTTTCCGATGCAGAAGCGCCTGAAGACAAAAATTCCGGAGATAAAAAATCGCGGATAACaa                                                                                           >  2:1332464/6‑99 (MQ=255)  tATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAt                                                                                          >  1:1749619/1‑100 (MQ=255)     ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                       >  1:353076/3‑100 (MQ=255)     cAGTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                       >  1:1511213/1‑100 (MQ=255)     cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                       <  2:2097784/100‑1 (MQ=255)      cgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                      >  2:2086534/2‑100 (MQ=255)      aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                      >  1:1565502/1‑100 (MQ=255)      aGTGATCAAATTCATCACGGTCAGGGGAGAAAAATTTTCGGGTGCAGAAGCGGCTGCAGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCaag                                                                                          <  2:2236909/96‑2 (MQ=255)        tgatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCAt                                                                                    >  2:900852/1‑100 (MQ=255)         gatgaAATTCATCACGGTCATGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCAt                                                                                    >  2:1371281/1‑99 (MQ=255)           tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATc                                                                                   >  1:21795/1‑98 (MQ=255)           tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                 >  1:1388621/1‑100 (MQ=255)            gaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCggg                                                                                >  1:1181365/1‑100 (MQ=255)             aaaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGa                                                                               >  2:547537/1‑100 (MQ=255)             aaaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAAACGACTGCGGGCCAAAATTGCGGGGAGGCGAAGTTGCGGATTGCAATACCCATCGGGa                                                                               <  2:555963/100‑1 (MQ=255)              tttttttctcGGTCAGGTGAGAAAAGTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                              <  1:579264/91‑1 (MQ=255)              aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                              <  1:216142/100‑1 (MQ=255)               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                             <  2:1115514/100‑1 (MQ=255)               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                             <  2:1060406/100‑1 (MQ=255)                ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGACCa                                                                            >  1:1467823/1‑100 (MQ=255)                ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTACGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                             >  2:536253/1‑99 (MQ=255)                 tcatccCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCTGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCAt                                                                           >  2:446310/1‑100 (MQ=255)                  catcaCGGTCAGGTGAGAAAAATTTTCCGCTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCCAAGTCGCGGATTGCAATACCCATCGGGATCATa                                                                          <  2:1345352/100‑1 (MQ=255)                   atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGCGATGCGAAGTCGCGTATTGCAATACCCATCGGGATCATaa                                                                         >  1:191800/1‑100 (MQ=255)                    tcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaaa                                                                        <  2:1003400/100‑1 (MQ=255)                    tcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaaa                                                                        <  1:492645/100‑1 (MQ=255)                       cGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACAt                                                                     >  1:2091794/1‑100 (MQ=255)                         gTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATg                                                                    <  2:248233/99‑1 (MQ=255)                           cAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGttt                                                                 <  1:1322847/100‑1 (MQ=255)                             ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                               >  1:1879226/1‑100 (MQ=255)                             ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                               >  1:2074499/1‑100 (MQ=255)                               tGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGa                                                             >  1:1283515/1‑100 (MQ=255)                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                          >  1:1743169/1‑100 (MQ=255)                                    aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                        <  1:157121/100‑1 (MQ=255)                                       aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                     >  1:1192127/1‑100 (MQ=255)                                        ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                                                    <  2:5916/100‑1 (MQ=255)                                        ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                                                    <  2:58215/100‑1 (MQ=255)                                        ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                                                    <  1:1737116/100‑1 (MQ=255)                                            ccGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                                                <  1:1554354/100‑1 (MQ=255)                                              ggTGCAGAACCGACTGCGGTACAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaa                                                >  2:557443/1‑98 (MQ=255)                                               gTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAcc                                             >  2:1514628/1‑100 (MQ=255)                                                 gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                           >  1:341418/1‑100 (MQ=255)                                                 gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                           <  1:905895/100‑1 (MQ=255)                                                  cAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                          <  1:2067279/100‑1 (MQ=255)                                                   ataaCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                         <  1:539057/98‑1 (MQ=255)                                                   aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                          >  1:980371/1‑99 (MQ=255)                                                   aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                         >  1:230320/1‑100 (MQ=255)                                                      aCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca                                      >  1:1667701/1‑100 (MQ=255)                                                        cGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac                                    <  1:1702276/100‑1 (MQ=255)                                                           cttctgTCCAAAATTCCGGTTATGCGAAGTCGCTGATTGCAAGACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaacaa                                  <  1:93844/94‑1 (MQ=255)                                                           cTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa                                 >  1:803354/1‑100 (MQ=255)                                                           cTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa                                 >  1:1310295/1‑100 (MQ=255)                                                            tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                <  2:1049166/100‑1 (MQ=255)                                                            tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                <  1:1866120/100‑1 (MQ=255)                                                             gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                               <  2:2074485/100‑1 (MQ=255)                                                             gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                               <  1:407719/100‑1 (MQ=255)                                                               ggTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATc                             <  2:648941/100‑1 (MQ=255)                                                               ggTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATc                             <  1:1924946/100‑1 (MQ=255)                                                                gTCCAAAATTCCGAGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcg                            >  1:1519992/1‑100 (MQ=255)                                                                 tCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAGACATcgc                           <  1:1171903/100‑1 (MQ=255)                                                                 tCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgc                           <  2:304295/100‑1 (MQ=255)                                                                   cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                         >  1:254568/1‑100 (MQ=255)                                                                   cAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                         >  1:1773976/1‑100 (MQ=255)                                                                       aTTCCGGGGATGCTAAGTCGCTGATTGCAATTCTCATCGGGATCATAAACATGTTTGCGATTCTGTGCTCAAATCTGTTTGCAACAAACATCGCGACCg                      <  2:392618/99‑1 (MQ=255)                                                                       aTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCg                      >  1:40433/1‑99 (MQ=255)                                                                        ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                    <  2:1081785/100‑1 (MQ=255)                                                                         tCCGGGGAGGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgca                   <  2:1940304/100‑1 (MQ=255)                                                                          ccGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcag                  <  1:779509/100‑1 (MQ=255)                                                                           cGGGGATGCGAAGTCGCGGATTGCAATACCCATCTGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                 <  2:1181361/100‑1 (MQ=255)                                                                           cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgca                   >  1:261598/1‑98 (MQ=255)                                                                             gggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCAc               <  1:339653/100‑1 (MQ=255)                                                                             gggATGCGAAGTCGCGGATTGCAATACCCATCGGGAACATAAACATGTTTGCGATACTTTGCTCAAAACCGCTGGCAACAAACATCGCGACACgcagc                 >  2:1564367/1‑98 (MQ=255)                                                                              ggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagca                >  2:1781190/1‑98 (MQ=255)                                                                                  gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAATACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa          <  2:1256545/100‑1 (MQ=255)                                                                                   cGAAGTCGCGGATTGCAATACCCATCTGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGTCAACAAACATCGCGACCGGCAGCACCATaa          <  2:537204/99‑1 (MQ=255)                                                                                   cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt         >  2:304191/1‑100 (MQ=255)                                                                                    gAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg        >  1:1109379/1‑100 (MQ=255)                                                                                    gAAGTCCCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg        <  1:2125535/100‑1 (MQ=255)                                                                                     aaGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGa       <  1:1694040/100‑1 (MQ=255)                                                                                        tCGCGGATTGCAATTCCCATCGGGATCATAAACATTTTTGCGATACTGTGCTCAAATCCGCTGGTAACAAACATCGCGTCCGGCAGCACCATAATGAACg    <  2:910565/100‑1 (MQ=255)                                                                                         cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGc   >  2:1653922/1‑100 (MQ=255)                                                                                           cgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACAATAATGaaagat  >  1:1303780/1‑95 (MQ=255)                                                                                                 |                                                                                             TTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCT  >  NC_000913/953706‑953894

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