Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	973,635	G→A	7.0%	V18V (GTG→GTA)	secF →	protein translocase subunit SecF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	973,635	0	G	A	7.0%	168.1 / 2.1	57	V18V (GTG→GTA)	secF	protein translocase subunit SecF
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

TTCACCGCCGTCATGGTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCTGTGCGTGCCGAGTTGGTC  >  NC_002947/973494‑973770                                                                                                                                               |                                                                                                                                        ttCACCGCCGTCATGGTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGctcctc                                                                                                                                   <  5:295834/148‑1 (MQ=255) ttCACCGCCGTCATGGTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCa                                                                                                                                  >  5:277568/1‑149 (MQ=255)        aggTCAGGTTGACCAATGCACGGGTCACCATGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGtt                                                                                                                           >  4:2323/3‑149 (MQ=255)         cgTCATGGTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTg                                                                                                                          <  6:277568/149‑1 (MQ=255)              tGGTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCt                                                                                                                     <  7:121560/149‑1 (MQ=255)               ggTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTgg                                                                                                                         <  1:18804/144‑1 (MQ=255)               ggTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTgg                                                                                                                         >  2:18804/1‑144 (MQ=255)                   agccGTGCATTTGTCAAGCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTc                                                                                                                 <  2:89414/146‑1 (MQ=255)                   aCCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAACATTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGtt                                                                                                                  <  1:110869/148‑1 (MQ=255)                          cATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCtggtg                                                                                                         <  3:240017/149‑1 (MQ=255)                          cATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCtggtg                                                                                                         >  4:240017/1‑149 (MQ=255)                                tCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCaaaa                                                                                                   <  1:302842/149‑1 (MQ=255)                                  aaCCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAg                                                                                                  <  7:344908/148‑1 (MQ=255)                                  aaCCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAg                                                                                                  >  8:344908/1‑148 (MQ=255)                                            gtgGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACtt                                                                                       >  3:67387/1‑149 (MQ=255)                                               gcgggcgTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTg                                                                                                                          >  8:249799/1‑111 (MQ=255)                                               gcgggcgTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTg                                                                                                                          <  7:249799/111‑1 (MQ=255)                                               gcgggcgTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCaaa                                                                                                    >  6:217554/1‑133 (MQ=255)                                               gcgggcgTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCaaa                                                                                                    <  5:217554/133‑1 (MQ=255)                                                 gggcgTGACATCAAGAAGCTGAGGGTTTGAGGGAGCTGAGATGAAAACCATCAA‑TTTCAAGGGTCTGCGCAATGTCGCGTTCGCCGTCAACGTGCTCCTCAGCGTGATGGAGCTGTACACCTGTTCGCAAAAGCGCCTGAACTTCggg                                                                                   >  8:316480/1‑148 (MQ=255)                                                     gTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGa                                                                              <  2:52106/149‑1 (MQ=255)                                                                 agCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCgggg                                                                   <  3:2323/148‑1 (MQ=255)                                                                  gCTGTGGGGTTGAGGGAGCTGCGATGAATACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCTCCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGt                                                                 <  7:316480/149‑1 (MQ=255)                                                                    tgtgGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCtggt                                                                                                          <  6:251972/106‑1 (MQ=255)                                                                    tgtgGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCtggt                                                                                                          >  5:251972/1‑106 (MQ=255)                                                                      tgtgTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCAGCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGTCTTCACCGGGGGTACGc                                                             <  4:67387/146‑1 (MQ=255)                                                                             gAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATc                                                        >  3:90234/1‑147 (MQ=255)                                                                               gggAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGc                                                                                               <  7:350030/106‑1 (MQ=255)                                                                               gggAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGc                                                                                               >  8:350030/1‑106 (MQ=255)                                                                                       cGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCATCTGGTGGCAAAATGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCt                                              <  3:287365/147‑1 (MQ=255)                                                                                             aaaCCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAgcg                                       <  4:90234/148‑1 (MQ=255)                                                                                               aCCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGccc                                    >  2:217666/1‑149 (MQ=255)                                                                                                ccATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACtt                                                                           <  1:247499/109‑1 (MQ=255)                                                                                                ccATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACtt                                                                           >  2:247499/1‑109 (MQ=255)                                                                                                  aTCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCgggg                                                                   >  1:133044/1‑115 (MQ=255)                                                                                                  aTCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCgggg                                                                   <  2:133044/115‑1 (MQ=255)                                                                                                   tCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGcc                                >  3:242228/1‑149 (MQ=255)                                                                                                   tCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGcc                                >  7:182149/1‑149 (MQ=255)                                                                                                          tCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGt                           >  5:156895/1‑148 (MQ=255)                                                                                                             tGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAg                       >  7:212289/1‑149 (MQ=255)                                                                                                              gggTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAgc                                        <  6:153487/131‑1 (MQ=255)                                                                                                              gggTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAgc                                        >  5:153487/1‑131 (MQ=255)                                                                                                               ggTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCCCCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAgg                      <  7:290699/148‑1 (MQ=255)                                                                                                                gtgtGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCt                                                            >  1:182490/1‑109 (MQ=255)                                                                                                                gtgtGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCt                                                            <  2:182490/109‑1 (MQ=255)                                                                                                                gtgtGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCt                    >  5:346548/1‑149 (MQ=255)                                                                                                                   ggcgcAATGGCGCGTTCGCCGTCACCTTGCTCCTCACCGTGTTTGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCtgtg                 <  4:237336/148‑1 (MQ=255)                                                                                                                      gcAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCTgtgcgt              >  6:311667/1‑149 (MQ=255)                                                                                                                        aaTGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCTgtgcgtgc            <  4:242228/149‑1 (MQ=255)                                                                                                                            tCGCGTTCGCCGTCACCGTACTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAg                       <  1:151282/134‑1 (MQ=255)                                                                                                                            tCGCGTTCGCCGTCACCGTACTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAg                       >  2:151282/1‑134 (MQ=255)                                                                                                                            tCGCGTTCGCCGTCACCGTACTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAg                       >  6:95840/1‑134 (MQ=255)                                                                                                                            tCGCGTTCGCCGTCACCGTACTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAg                       <  5:95840/134‑1 (MQ=255)                                                                                                                                 ttCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCTGTGCGTGCCGAGTTGGt   >  5:248920/1‑149 (MQ=255)                                                                                                                                   cGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCTGTGCGTGCCGAGTTGGTc  <  6:346548/148‑1 (MQ=255)                                                                                                                                     ccGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAgg                                                                                                 <  6:261310/51‑1 (MQ=255)                                                                                                                                     ccGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAgg                                                                                                 >  5:261310/1‑51 (MQ=255)                                                                                                                                               |                                                                                                                                        TTCACCGCCGTCATGGTGACCCGTGCATTGGTCAACCTGACCTGTGGCGGGCGTGACATCAAGAAGCTGTGGGTTTGAGGGAGCTGCGATGAAAACCATCAA‑TTTCATGGGTGTGCGCAATGTCGCGTTCGCCGTCACCGTGCTCCTCACCGTGTTGGCGCTGTTCAGCTGGTGGCAAAAGGGCCTGAACTTCGGGCTGGACTTCACCGGGGGTACGCTGATCGAGCTGACCTACGAGCGCCCGGCCGACGTCAAGGCTGTGCGTGCCGAGTTGGTC  >  NC_002947/973494‑973770

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