Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	108,835	A→G	7.1%	Y186C (TAC→TGC)	secA →	preprotein translocase subunit, ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	108,835	0	A	G	7.1%	164.6 / 3.3	56	Y186C (TAC→TGC)	secA	preprotein translocase subunit, ATPase
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

TGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGAAA  >  NC_000913/108692‑108981                                                                                                                                                |                                                                                                                                                   tGACGCCGAAAACACCCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTgc                                                                                                                                               >  7:207441/1‑149 (MQ=255)      ccGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGac                                                                                                                                          <  4:72285/149‑1 (MQ=255)          aaaCAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGTGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGacaaca                                                                                                                                      <  4:101238/149‑1 (MQ=255)              acccGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCTGGCATTCGAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGc                                                                                                                                  <  8:28372/147‑1 (MQ=255)                       cTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGccc                                                                                                                         >  3:205041/1‑149 (MQ=255)                           ttGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTg                                                                                                                       >  6:238398/1‑147 (MQ=255)                             gAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTGCCTGCGCGACAACATGGCGTTCAGCCCTgaag                                                                                                                    >  7:40564/1‑148 (MQ=255)                             gAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTGCCTGCGCGACAACATGGCGTTCAGCCCTgaag                                                                                                                    >  7:36538/1‑148 (MQ=255)                                 tCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgt                                                                                                               >  4:63868/1‑149 (MQ=255)                                  ccTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgt                                                                                                               <  7:143245/148‑1 (MQ=255)                                    ttGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTa                                                                                                            >  8:251957/1‑149 (MQ=255)                                    ttGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTa                                                                                                            >  4:229386/1‑149 (MQ=255)                                        ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              >  2:175752/1‑143 (MQ=255)                                        ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              >  6:131190/1‑143 (MQ=255)                                        ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              <  5:131190/143‑1 (MQ=255)                                        ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              <  1:175752/143‑1 (MQ=255)                                              tGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCCCTTACGGGACGAACGACGAATACGGCTTTGACTACCTGCGCGACGACCTGGCGTTCCGCCCTGAAGAACGGGTACAGCGTAAAc                                                                                                  >  6:155998/1‑149 (MQ=255)                                                   gTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCAc                                                                                             >  6:164059/1‑149 (MQ=255)                                                      tCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTAt                                                                                          <  1:141704/149‑1 (MQ=255)                                                           cTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                                                                                     >  7:137134/1‑149 (MQ=255)                                                              ccGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCtt                                                                                    <  8:36446/147‑2 (MQ=255)                                                                gggCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCtggtg                                                                                 >  1:64273/1‑148 (MQ=255)                                                                   cATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGAc                                                                              >  8:88179/1‑148 (MQ=255)                                                                        cAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg                                                                        >  3:178830/1‑149 (MQ=255)                                                                        cAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg                                                                        <  5:164059/149‑1 (MQ=255)                                                                          tgcccGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGa                                                                      <  4:104792/146‑1 (MQ=255)                                                                             ccGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTc                                                                   >  4:65062/1‑149 (MQ=255)                                                                              cGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTc                                                                   >  2:41422/1‑148 (MQ=255)                                                                                gCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCAt                                                                >  3:51166/1‑149 (MQ=255)                                                                                 cAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc                                                               <  4:51166/149‑1 (MQ=255)                                                                                  aaaGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATcc                                                              >  2:167519/1‑149 (MQ=255)                                                                                         gAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGa                                                       >  5:80382/1‑149 (MQ=255)                                                                                          aaGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGAt                                                      <  3:63868/149‑1 (MQ=255)                                                                                             cTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGaa                                                   <  2:64273/149‑1 (MQ=255)                                                                                                        gACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACAcc                                        >  1:6648/1‑149 (MQ=255)                                                                                                         aCATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCg                                       <  7:251957/149‑1 (MQ=255)                                                                                                          cATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGc                                      <  1:41422/149‑1 (MQ=255)                                                                                                                ttACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCa                                >  1:66413/1‑149 (MQ=255)                                                                                                                ttACGGTACGAACAACGAATACGGATTTGACTACCTGCGCGACAACATGGCGTTCAGACCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATTAAGCGAGTACACCGCTGATCa                                >  7:217120/1‑149 (MQ=255)                                                                                                                  aCGGTACGCACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGACCGTGTACAGCGTAAACTGCACTATGCGCGGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCa                                <  5:155998/147‑1 (MQ=255)                                                                                                                     gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATc                                 <  5:167420/143‑1 (MQ=255)                                                                                                                     gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATc                                 >  6:167420/1‑143 (MQ=255)                                                                                                                     gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTcc                           >  5:89392/1‑149 (MQ=255)                                                                                                                        cGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTccgg                         <  2:66413/148‑1 (MQ=255)                                                                                                                          aacaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGt                                             >  7:7212/1‑126 (MQ=255)                                                                                                                          aacaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGt                                             <  8:7212/126‑1 (MQ=255)                                                                                                                             aacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTccggcccgg                    >  3:82848/1‑148 (MQ=255)                                                                                                                                 aaTACGGCTTTTACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTTGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCagaa               <  8:217120/149‑1 (MQ=255)                                                                                                                                 aaTACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCagaa               <  3:65062/149‑1 (MQ=255)                                                                                                                                  aTACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCagaag              <  4:169976/149‑1 (MQ=255)                                                                                                                                  aTACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCag                 >  7:117267/1‑146 (MQ=255)                                                                                                                                     cGGCTTTGACTACCTGCGCGACAACCTGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGAc            <  7:88179/148‑1 (MQ=255)                                                                                                                                           tGACTGCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTc       <  8:36538/147‑1 (MQ=255)                                                                                                                                           tGACTGCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTc       <  8:40564/147‑1 (MQ=255)                                                                                                                                              cTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGaaa  >  7:251262/1‑149 (MQ=255)                                                                                                                                              cTACCTGCGCGACAACATGGCGGTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGaaa  <  5:130218/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   TGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGAAA  >  NC_000913/108692‑108981

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