Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	953,802	A→G	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	95.7%	226.8 / ‑2.4	70	V222A (GTA→GCA)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (33/34);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGC  >  NC_000913/953705‑953893                                                                                                  |                                                                                            gTTATCAGTGATGAAATTCATCACGGTCAGTTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCCAAATTCCGGTGATGCGAAGTCGCGGATTGCa                                                                                           <  1:801031/100‑1 (MQ=255) gTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATGTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCa                                                                                           <  1:1637651/100‑1 (MQ=255)  ttATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCaa                                                                                          <  1:802735/100‑1 (MQ=255)  ttATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACaa                                                                                          >  2:731678/1‑100 (MQ=255)    gaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATa                                                                                        >  2:1002796/5‑100 (MQ=255)    aTCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATa                                                                                        <  1:329725/100‑1 (MQ=255)       aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                     <  2:515895/100‑1 (MQ=255)       aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                     <  1:236068/100‑1 (MQ=255)       aGTGATGAAATTCATCACGGGCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATAccc                                                                                     >  1:1730222/1‑100 (MQ=255)          gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATc                                                                                  >  2:1132411/1‑100 (MQ=255)            tgaAATTTATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgn                                                                                <  1:975372/100‑2 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                <  2:902991/100‑1 (MQ=255)            tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                <  1:1129628/100‑1 (MQ=255)              aatttCATCACGGTCAGGTGAGAAAAATTTTCCAGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGa                                                                              >  2:263680/4‑100 (MQ=255)              aaaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGa                                                                              <  1:1752155/100‑1 (MQ=255)                 ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                           <  2:792047/100‑1 (MQ=255)                 ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                           <  2:1113714/100‑1 (MQ=255)                 ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                           <  1:70003/100‑1 (MQ=255)                 ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                           >  1:361074/1‑100 (MQ=255)                  tcatcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCAt                                                                          >  1:610845/1‑100 (MQ=255)                   catcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGAAATACACATCAGGATCATa                                                                         >  2:70699/1‑100 (MQ=255)                    atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTACAAAATTCCGGGGATGCGAAGTCGCAGATTGCAATACCCATCAGGATCATaa                                                                        >  1:1044510/1‑100 (MQ=255)                       aCGGTCAGGTGAGAAAACTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCTAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                        <  1:1565179/97‑1 (MQ=255)                        cGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACAt                                                                    >  2:1338369/1‑100 (MQ=255)                             aGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGttt                                                                >  1:144241/1‑99 (MQ=255)                              ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                              <  2:1594874/100‑1 (MQ=255)                              ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                              >  1:964407/1‑100 (MQ=255)                              ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                              >  2:1348759/1‑100 (MQ=255)                                 gagaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                           >  2:742499/1‑100 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATa                                                          >  2:1180998/1‑99 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                         >  2:772922/1‑100 (MQ=255)                                   gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                         >  1:729984/1‑100 (MQ=255)                                     aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                       <  1:1230217/100‑1 (MQ=255)                                       aaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtgtg                                                     <  1:742495/100‑1 (MQ=255)                                        aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                    <  1:1699425/100‑1 (MQ=255)                                        aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                    >  1:1187609/1‑100 (MQ=255)                                         ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCGGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                                                   >  1:1135467/1‑100 (MQ=255)                                         ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                                                   <  1:1690500/100‑1 (MQ=255)                                           ttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCa                                                 <  2:102482/100‑1 (MQ=255)                                           ttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCa                                                 >  2:1164808/1‑100 (MQ=255)                                               ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                             <  2:818779/100‑1 (MQ=255)                                                 tGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg                                           >  1:1588540/1‑100 (MQ=255)                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATACACATGTTTGCGATACTGTGCTCAAAACCGCt                                         >  1:1277796/1‑99 (MQ=255)                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                        >  1:498673/1‑100 (MQ=255)                                                         cGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaa                                    <  2:172895/99‑1 (MQ=255)                                                           actgttgTCCAAAATTCCGGGTATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                         <  2:1485720/86‑1 (MQ=255)                                                             tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                               >  2:1757419/1‑100 (MQ=255)                                                              gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                              >  1:676815/1‑100 (MQ=255)                                                              gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                              >  2:1389698/1‑100 (MQ=255)                                                              gCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca                              >  1:1179328/1‑100 (MQ=255)                                                                ggTCCAAAATTCCGTGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTGTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATc                            <  2:589924/100‑1 (MQ=255)                                                                 gTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcg                           >  1:1451190/1‑100 (MQ=255)                                                                   ccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                         >  2:1339406/1‑100 (MQ=255)                                                                   ccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                         >  1:140093/1‑100 (MQ=255)                                                                     aaaaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAATGTGCTCAAACCCGCTGGCAACAAACATcg                           >  1:1780989/1‑96 (MQ=255)                                                                       aaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCg                     <  1:104206/100‑1 (MQ=255)                                                                        aTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGAcc                      <  2:482869/98‑1 (MQ=255)                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                   <  1:551046/100‑1 (MQ=255)                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                   <  2:621940/100‑1 (MQ=255)                                                                          tCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgca                  <  2:1215381/100‑1 (MQ=255)                                                                            cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                >  1:548748/1‑100 (MQ=255)                                                                            cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                <  2:1280191/100‑1 (MQ=255)                                                                             ggggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCAGGAAACAAACATCACGACCGgcagca               >  2:1632727/1‑100 (MQ=255)                                                                                gATGCGAAGTCGCGGATTGCAATACCCATCGGCATCATAAACATGTTTGCGATACTGTGCTCAAATCCGCTGGCACCAAACATCGCGACCGGCAGCAcc             <  1:1716230/99‑1 (MQ=255)                                                                                   gCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATaa         >  1:891695/1‑100 (MQ=255)                                                                                   gCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCACCAAACATCGCGACCGGCCGCACCATaa         >  2:490809/1‑100 (MQ=255)                                                                                    cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt        <  1:772917/100‑1 (MQ=255)                                                                                    cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt        <  1:579274/100‑1 (MQ=255)                                                                                       aGTCGCGGATTGCAATATCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt        <  1:1717753/97‑1 (MQ=255)                                                                                          cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGc  >  1:1638056/1‑100 (MQ=255)                                                                                                  |                                                                                            GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGC  >  NC_000913/953705‑953893

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