Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,835,836	A→T	7.5%	L152M (TTG→ATG) ‡	yohC ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,835,836	0	A	T	7.5%	152.3 / 3.1	53	L152M (TTG→ATG) ‡	yohC	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base T (2/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

TTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGAC  >  NC_002947/1835688‑1835984                                                                                                                                                     |                                                                                                                                                     ttCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCaa                                                                                                                                                      >  8:80455/1‑149 (MQ=255)   cTAGTTGGTGTAGACGGGCCCCACACCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCAGCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGc                                                                                                                                                    >  5:294464/1‑149 (MQ=255)      gTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGc                                                                                                                                                    <  1:136490/146‑1 (MQ=255)           tgtAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCTACGGTATGTTCATGAACGTCGGCAAGCCGACATa                                                                                                                                             >  7:27320/1‑148 (MQ=255)            gtAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATa                                                                                                                                             >  1:302830/1‑147 (MQ=255)             ttgACGTGCCCCAAGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCTTGCCGAcgttca                                                                                                                                           <  6:260240/146‑6 (MQ=255)               gACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATAcagca                                                                                                                                        <  4:147404/149‑1 (MQ=255)                  gggCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAgg                                                                                                                                      >  8:131702/1‑148 (MQ=255)                        cACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCg                                                                                                                               >  4:82659/1‑149 (MQ=255)                           gCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCTACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCgtgt                                                                                                                            <  8:27320/149‑1 (MQ=255)                           gCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCgtgt                                                                                                                            <  1:357613/149‑1 (MQ=255)                           gCCCAGACCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCgtgt                                                                                                                            <  2:275746/149‑1 (MQ=255)                            cccAGTCCCAAGATAATCACGGTAGCGGCCATGATTGACACCAGTACTACCAAGCCCACTGCCAGCACCGAGATGGAGAACAGGAAGCCTTCATCCGACGGTATGTTAATGAACGTCGGCAACCCGACATACACCAGCTACACCGTCTa                                                                                                                           >  1:231440/1‑149 (MQ=255)                               aGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCa                                                                                                                        <  2:245731/149‑1 (MQ=255)                                    ccAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATgg                                                                                                                   >  4:66897/1‑149 (MQ=255)                                    ccAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCAACTGCCACCACCGACATGCAGAACAGCAAGCGTTCATCACACGTTATGTTCATGAACCTCGGCAAGCCGACATAAAGCAGGTACACCGTGTTGCAGATgg                                                                                                                   >  6:6435/1‑149 (MQ=255)                                     cAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGAACCTAAACCGTGTATCAGATgac                                                                                                                  >  6:259521/1‑147 (MQ=255)                                           aTCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAg                                                                                                             <  7:131702/148‑1 (MQ=255)                                           aTCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAg                                                                                                             <  3:82659/148‑1 (MQ=255)                                                            gATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGcc                                                                                             >  4:4384/1‑147 (MQ=255)                                                                  gaacaGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCca                                                                                       >  2:245888/4‑147 (MQ=255)                                                                   accaGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAg                                                                                    >  8:304064/1‑149 (MQ=255)                                                                          ctacCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCggg                                                                             <  7:24432/149‑1 (MQ=255)                                                                          ctacCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCgg                                                                              <  7:80455/148‑1 (MQ=255)                                                                           tacCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCggg                                                                             >  8:110173/1‑148 (MQ=255)                                                                                         cctgCACGGAGCTGGAGAACTGGAAGCCTTGATCCGTGGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGTTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCTTGCGCAGCGACAGGTGCGGGTAAAGCGCCGCCatg                                                              <  2:231440/146‑3 (MQ=255)                                                                                          cAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACAcc                                                                                                                                <  6:280105/82‑1 (MQ=255)                                                                                          cAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACAcc                                                                                                                                >  5:280105/1‑82 (MQ=255)                                                                                          cAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGgc                                                             >  2:39760/1‑149 (MQ=255)                                                                                             cACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGgccgc                                                          >  7:40607/1‑149 (MQ=255)                                                                                                      ggAGAACAGGAAGCCTTCATTCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAAc                                                  >  1:353327/1‑148 (MQ=255)                                                                                                           aCAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAgggg                                            >  8:286194/1‑149 (MQ=255)                                                                                                               gAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTgg                                        >  5:160924/1‑149 (MQ=255)                                                                                                                  gCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAggggg                                           >  1:360859/1‑143 (MQ=255)                                                                                                                  gCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAggggg                                           <  2:360859/143‑1 (MQ=255)                                                                                                                   ccTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGgt                                    <  2:302830/149‑1 (MQ=255)                                                                                                                    cTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCGGGCCGCCGATGAACAAGGGGGTGGCGgtg                                   >  4:218073/1‑149 (MQ=255)                                                                                                                    cTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGgtg                                   >  1:137624/1‑149 (MQ=255)                                                                                                                         tCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGt                              <  1:39760/149‑1 (MQ=255)                                                                                                                           cGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTgg                            <  5:6435/149‑1 (MQ=255)                                                                                                                                 tATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGGGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAg                      >  3:277239/1‑149 (MQ=255)                                                                                                                                     ttCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGa                   <  3:66897/148‑1 (MQ=255)                                                                                                                                      tCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATg                 >  3:61659/1‑149 (MQ=255)                                                                                                                                       gctgtatgTCGGCTTGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCCGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATg                 >  2:67873/8‑148 (MQ=255)                                                                                                                                        ctgtatgTCGGCTTGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCa               <  1:67873/143‑1 (MQ=255)                                                                                                                                        aTGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCa               <  7:304064/149‑1 (MQ=255)                                                                                                                                        aTGAACGTCGGCAAGCCGACATACAGCAGGTAAACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCa               >  6:115868/1‑149 (MQ=255)                                                                                                                                          gAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCaaag                       <  6:185178/139‑2 (MQ=255)                                                                                                                                          gAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCaaag                       >  5:185178/1‑138 (MQ=255)                                                                                                                                           tatgTCGGCTTGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTg            >  4:83472/4‑149 (MQ=255)                                                                                                                                           aaCGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCATCCACAGGTGCGGGTAAAGCGGAGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATTCATTg            <  5:259521/149‑1 (MQ=255)                                                                                                                                            aCGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTGAGTGGCAAAAGCGATGCATTgc           >  3:43057/1‑149 (MQ=255)                                                                                                                                                 ggCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCAt      >  1:336762/1‑149 (MQ=255)                                                                                                                                                  gCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATg     <  6:294464/149‑1 (MQ=255)                                                                                                                                                     aGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGAc  <  4:61659/149‑1 (MQ=255)                                                                                                                                                     aGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGAc  <  1:245888/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGAC  >  NC_002947/1835688‑1835984

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