Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,249,628	G→A	71.8%	G85D (GGC→GAC)	yqjD →	membrane‑anchored ribosome‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,249,628	0	G	A	71.8%	96.5 / 52.5	76	G85D (GGC→GAC)	yqjD	membrane‑anchored ribosome‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (11/10);  new base A (31/23);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

ATCGCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCT  >  NC_000913/3249529‑3249726                                                                                                    |                                                                                                   aTCGCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTggg                                                                                                    <  2:411358/100‑1 (MQ=255) aTCGCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTggg                                                                                                    <  2:2030187/100‑1 (MQ=255)    gCCTGGGTGAAAGCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACAt                                                                                                 >  1:1469987/1‑100 (MQ=255)    gCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCAt                                                                                                 >  2:1775693/1‑100 (MQ=255)       tGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTgg                                                                                              >  2:1922965/1‑100 (MQ=255)        gggTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGgc                                                                                             >  2:1677409/1‑100 (MQ=255)        gggTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTAGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGgc                                                                                             >  1:1496898/1‑100 (MQ=255)         ggTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGgcg                                                                                            <  2:749926/100‑1 (MQ=255)          gTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGgcgc                                                                                           >  1:249148/1‑100 (MQ=255)          gTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGgcgc                                                                                           <  2:379523/100‑1 (MQ=255)           agAAACCGGTGATGCCATTGCCAAACAAACCCGTCTCGCGGCGACGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTgg                                                                                              >  2:550008/2‑96 (MQ=255)            gAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTTGACGGGCGTGGACATTGGCGCTg                                                                                         <  1:2089523/100‑1 (MQ=255)             aaaCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGc                                                                                        <  1:1381893/100‑1 (MQ=255)             aaaCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTg                                                                                         >  1:1672589/1‑99 (MQ=255)             aaaCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGc                                                                                        <  1:662390/100‑1 (MQ=255)             aaaCCGGGGATGCCATTGCCAAACAAATCCGTGGCGCGGCGGCTCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGc                                                                                        <  2:1211663/100‑1 (MQ=255)              aaCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCa                                                                                       >  1:813785/1‑100 (MQ=255)                        gCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAg                                                                             <  2:956339/100‑1 (MQ=255)                         ccATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGt                                                                            >  1:1740212/1‑100 (MQ=255)                              gCCAAACAAACCCGTGTCGCGGCGTCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCg                                                                       >  1:1571896/1‑100 (MQ=255)                              gCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGCTGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTc                                                                        >  2:784475/1‑99 (MQ=255)                              gCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCg                                                                       >  2:1211605/1‑100 (MQ=255)                              gCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGATGTAGTGCTCg                                                                       >  1:320675/1‑100 (MQ=255)                              gCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTc                                                                        >  1:337148/1‑99 (MQ=255)                                caaacaaaCCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGc                                                                     <  1:891097/100‑1 (MQ=255)                                  aacaaacaCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGt                                                                   >  1:1169871/1‑100 (MQ=255)                                     aaacCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTct                                                                <  2:58213/100‑1 (MQ=255)                                       acCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                              >  2:2140423/1‑100 (MQ=255)                                        cccGTTTCGTTGCGGCGCGTGCGGTTGATTATGTGCGCGAAAATCCGTGGACGGGCTTTGTCATTGGCGCTGCAATCTGTTTAGTGCTCGGCGTTctgc                                                              <  2:296509/99‑1 (MQ=255)                                        cccGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                              <  1:1659632/99‑1 (MQ=255)                                        cccGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGACGTTctgct                                                             >  1:807187/1‑100 (MQ=255)                                          cGTGTTGCGGCGGCGCGTGTCGATGAGTATGTGCGCGAAGATCCGTGGACGGGCGTGGACATTGGCGTTGCAATCGGTGTAGTGCTCGGCGTTct                                                                <  2:1012071/95‑1 (MQ=255)                                               cgcgGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGcgt                                                      >  1:1065675/1‑100 (MQ=255)                                                gcggcggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGcgtc                                                     <  1:321623/100‑1 (MQ=255)                                                gcggcggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTcgcggc                                                     >  2:1895239/1‑98 (MQ=255)                                                gcggcggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGcgtc                                                     >  2:1365411/1‑100 (MQ=255)                                                  ggcggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGATGTAGTGCTCGGCGTTCTGCTGTCGcgtcgt                                                   <  2:320676/100‑1 (MQ=255)                                                  ggcggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCCATCGGTGTAGTGCTCGGCGTTCTGCTGTCGcgtcgt                                                   >  1:1260180/1‑100 (MQ=255)                                                    cggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTa                                                 <  1:1777083/100‑1 (MQ=255)                                                    cggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTa                                                 <  2:1486186/100‑1 (MQ=255)                                                      gcgcGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACTGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAAt                                               <  2:1245691/100‑1 (MQ=255)                                                       cgcgTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAAtt                                              >  1:1543942/1‑100 (MQ=255)                                                       cgcgTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAAtt                                              <  2:1483305/100‑1 (MQ=255)                                                          gTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATg                                           <  1:1354638/100‑1 (MQ=255)                                                           tGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATgg                                          >  1:2169062/1‑100 (MQ=255)                                                            gCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGc                                         <  1:444532/100‑1 (MQ=255)                                                             ccGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCg                                        >  1:1172698/1‑100 (MQ=255)                                                             ccGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCg                                        <  2:2169047/100‑1 (MQ=255)                                                               gATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGa                                      >  1:775171/1‑100 (MQ=255)                                                               gATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGAACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTGGTTAATTATGGCGGa                                      <  2:1830979/100‑1 (MQ=255)                                                                    gTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCCATCGGTGTATTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACt                                  >  1:75502/1‑99 (MQ=255)                                                                    gTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtc                                 >  1:2067207/1‑100 (MQ=255)                                                                     tATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtca                                <  1:150859/100‑1 (MQ=255)                                                                       tgtgCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtcatc                              >  1:738266/1‑100 (MQ=255)                                                                         tgCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCAc                            >  2:2123625/1‑100 (MQ=255)                                                                           cgcgAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGc                          <  2:660340/100‑1 (MQ=255)                                                                            gcgAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGca                         <  2:346785/100‑1 (MQ=255)                                                                             cgAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGcac                        >  1:1989968/1‑100 (MQ=255)                                                                             cgAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGATCTCGGCGTTCTCCGGTCGCTTCAATAAAAAAGGCGGACACTCATCACGcac                        >  2:1604901/1‑100 (MQ=255)                                                                              cgaaaTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGATCTGCTGTCGCGTCGTTAATTATGGCGGAAACTCATCACGcac                        >  1:2102586/3‑99 (MQ=255)                                                                                aaaTCCGTGGACGGGCGTGGACATTGGCACTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAg                     <  2:1249497/100‑1 (MQ=255)                                                                                 aaTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAgg                    <  2:2067193/100‑1 (MQ=255)                                                                                 aaTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAgg                    <  2:242191/100‑1 (MQ=255)                                                                                  aTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAggg                   <  2:1421062/100‑1 (MQ=255)                                                                                     cGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCAAGCACAAGGGccc                <  1:1471410/100‑1 (MQ=255)                                                                                     cGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGccc                >  2:1936354/1‑100 (MQ=255)                                                                                     cGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGccc                >  1:1350405/1‑100 (MQ=255)                                                                                      gTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCg               >  1:1170918/1‑100 (MQ=255)                                                                                        ggACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCgg              <  2:738266/99‑1 (MQ=255)                                                                                        ggACGGGCGTAGACATTGACGCTGCAATCGGAGTAGTGCTCGGCGATATGCTGTCGCGTAGTTAATAATGGCGGACACTCATCACGCACAAGGGCCCg               >  2:493053/1‑98 (MQ=255)                                                                                           cGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTTCTCGGCGTTCTGCTGTCCCTTCGTTAATTATGGAGGACACTCATCACGCACAAGGGCCCGGTaaa          >  2:1551483/1‑100 (MQ=255)                                                                                           cGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTaaa          >  2:1987786/1‑100 (MQ=255)                                                                                           cGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTaaa          >  1:1370151/1‑100 (MQ=255)                                                                                             ggCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAg        <  1:2140438/100‑1 (MQ=255)                                                                                              gCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACCCATCACGCACAAGGGCCCGGTAAAAGc       >  1:1663747/1‑100 (MQ=255)                                                                                              gCGTGGACATTGGCGCTGCAATCGGTGTAGTGTTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAg        <  1:136738/99‑1 (MQ=255)                                                                                              gCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAg        <  1:1834671/99‑1 (MQ=255)                                                                                                 tGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGtt    >  1:597781/1‑100 (MQ=255)                                                                                                   ggCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAATAGCGTTCt  >  1:1554497/1‑100 (MQ=255)                                                                                                    |                                                                                                   ATCGCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCT  >  NC_000913/3249529‑3249726

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