Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,496,279 G→A 8.4% R331H (CGC→CAC)  ydcJ → DUF1338 domain‑containing protein YdcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,496,2790GA 8.4% 137.7 / 2.8 48R331H (CGC→CAC) ydcJDUF1338 domain‑containing protein YdcJ
Reads supporting (aligned to +/- strand):  ref base G (14/29);  new base A (2/2);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTGGTTG  >  NC_000913/1496130‑1496417
                                                                                                                                                     |                                                                                                                                          
gggCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCg                                                                                                                                            >  1:1039443/1‑150 (MQ=255)
   cACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCAc                                                                                                                                         >  1:305693/1‑150 (MQ=255)
      gCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACttt                                                                                                                                      <  1:1347759/150‑1 (MQ=255)
      gCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACttt                                                                                                                                      <  1:472753/150‑1 (MQ=255)
           ccgcccgcGTTGCTGAAATTGAGCAGCGTGGCGGGGCATTACAGCCGAAAGGGCGACAAATGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCtt                                                                                                                                 <  2:375631/150‑2 (MQ=255)
                 gcTTTGGTGAAATTGATCAGCGTGGCGTGGCATTAACGCCGAAAGTGCTACAACTGTATTATGATCTTCTGCGTAACGCTGGAACCGTGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTg                                                                                                                           <  2:1096370/150‑1 (MQ=255)
                   tttGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAg                                                                                                                         <  2:1039443/150‑1 (MQ=255)
                      ggTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGttt                                                                                                                      <  2:684217/150‑1 (MQ=255)
                            aTTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATg                                                                                                                <  2:1193726/150‑1 (MQ=255)
                                      gtggcgtggcATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCTCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCaat                                                                                                      <  2:146865/150‑2 (MQ=255)
                                             ggcATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGc                                                                                               <  1:669217/150‑1 (MQ=255)
                                                         gAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTa                                                                                   >  2:59203/1‑150 (MQ=255)
                                                              ttctataaCTTTATTATTATCTTCTGCGTAACGCTTGAACCGGGCAGGATAATCTCACTCACCAAATTCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTc                                                                              <  1:628605/144‑1 (MQ=255)
                                                              ggCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTc                                                                              <  1:1237173/150‑1 (MQ=255)
                                                                     aCTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGcc                                                                       >  2:793969/1‑150 (MQ=255)
                                                                         tATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCg                                                                   <  1:59203/150‑1 (MQ=255)
                                                                          atgatgATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCgg                                                                  >  1:630888/1‑150 (MQ=255)
                                                                          atgatgATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCgg                                                                  >  2:641514/1‑150 (MQ=255)
                                                                          atgatgATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCgg                                                                  <  1:46924/150‑1 (MQ=255)
                                                                           tgatgaTCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCggg                                                                 <  2:354417/150‑1 (MQ=255)
                                                                           tgatgaTCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCggg                                                                 <  2:251288/150‑1 (MQ=255)
                                                                               gatCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGCAGCCTTCGGGTGAg                                                             >  2:269296/1‑150 (MQ=255)
                                                                                   ttctGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACg                                                                         <  2:339798/134‑1 (MQ=255)
                                                                                   ttctGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACg                                                                         >  1:339798/1‑134 (MQ=255)
                                                                                    tctGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCa                                                        <  1:1375695/150‑1 (MQ=255)
                                                                                      tGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATc                                                      <  2:257557/150‑1 (MQ=255)
                                                                                         gTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGCAAACCTTACGCACTTTTCCTGACAGTGACTTTTTAATGCGTCAGCAAGGGTTGGCAGGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCTTCGTc                                                   <  2:980888/150‑1 (MQ=255)
                                                                                                       gggCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCt                                     >  2:1061062/1‑150 (MQ=255)
                                                                                                        ggCAGGATAATCTCACTCACCAAATGCATTTACAGGACACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGGCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTg                                    <  2:1063863/150‑1 (MQ=255)
                                                                                                         gCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTgg                                   <  2:798090/150‑1 (MQ=255)
                                                                                                            ggATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGgaga                                >  1:37616/1‑150 (MQ=255)
                                                                                                               tAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGa                             >  2:1211739/1‑150 (MQ=255)
                                                                                                                aaTCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGAt                            <  2:305693/150‑1 (MQ=255)
                                                                                                                  tctcACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATcc                          <  1:385981/150‑1 (MQ=255)
                                                                                                                      actcacCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTCAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCTGTTACTGTCTGAGGCCTTCGGGTGAGGCTCATGGTCAGTCGATTCATCCTTGAGAGTATCCACAg                      <  2:830324/150‑1 (MQ=255)
                                                                                                                        tcacCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGcc                    >  2:667746/1‑150 (MQ=255)
                                                                                                                         cacCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGGATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGccc                   <  1:751135/150‑1 (MQ=255)
                                                                                                                            cAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTa                <  1:793969/150‑1 (MQ=255)
                                                                                                                              aaTGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAAt              >  2:106248/1‑150 (MQ=255)
                                                                                                                                 gCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGa           >  2:395845/1‑150 (MQ=255)
                                                                                                                                       aCAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTgg     <  2:281039/150‑1 (MQ=255)
                                                                                                                                         aGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGGATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTGGtt   <  1:526738/150‑1 (MQ=255)
                                                                                                                                         aGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTGGtt   <  1:521709/150‑1 (MQ=255)
                                                                                                                                         aGGAAACCTTCCACACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTgg     >  1:973636/1‑148 (MQ=255)
                                                                                                                                         aGGAAACCTTCCACACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTgg     >  1:978628/1‑148 (MQ=255)
                                                                                                                                         aGGAAACCTTCCACACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTgg     <  2:973636/148‑1 (MQ=255)
                                                                                                                                         aGGAAACCTTCCACACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTgg     <  2:978628/148‑1 (MQ=255)
                                                                                                                                          ggAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTGGTTg  <  2:630888/150‑1 (MQ=255)
                                                                                                                                          ggAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTGGTTg  >  1:233149/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
GGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTACCGTCTGACGCCTTCGGGTGAGGCGCATCGTCAGGCGATTCATCCTGGAGACGATCCACAGCCCTTAATTGAACGTGGTTG  >  NC_000913/1496130‑1496417

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