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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,835,835	0	A	T	7.5%	152.4 / 2.8	53	L152* (TTG→TAG) ‡	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 = 8.60e-01

TTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATG  >  NC_002947/1835688‑1835981                                                                                                                                                    |                                                                                                                                                   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                                                                                                                             >  5:280105/1‑82 (MQ=255)                                                                                          cAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACAcc                                                                                                                             <  6:280105/82‑1 (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)                                                                                                                                                    |                                                                                                                                                   TTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAACAGGAAGCCTTCATCCGACGGTATGTTCATGAACGTCGGCAAGCCGACATACAGCAGGTACACCGTGTAGCAGATGGCGGCAGTGCCGATCAGCATGCCCAGCCACAGGTGCGGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATG  >  NC_002947/1835688‑1835981

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