Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,596,326 C→G 7.7% R227P (CGG→CCG)  nemA ← FMN‑linked N‑ethylmaleimide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,596,3260CG 7.7% 156.4 / 2.4 51R227P (CGG→CCG) nemAFMN‑linked N‑ethylmaleimide reductase
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base G (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

CCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTGCAGCAGGTAACCGTGGGCGGCGTG  >  NC_002947/3596179‑3596471
                                                                                                                                                   |                                                                                                                                                 
ccACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCg                                                                                                                                                  >  4:154721/1‑149 (MQ=255)
   ccGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTc                                                                                                                                               >  3:106786/1‑149 (MQ=255)
      gCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGgc                                                                                                                                            <  7:82621/149‑1 (MQ=255)
               ggCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGtt                                                                                                                                   <  3:154721/149‑1 (MQ=255)
               ggCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGtt                                                                                                                                   <  1:41971/149‑1 (MQ=255)
                        aGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGc                                                                                                                          <  3:128798/149‑1 (MQ=255)
                                             cgcTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATCCCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGca                                                                                                      <  4:106786/148‑1 (MQ=255)
                                              gcTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGcac                                                                                                     >  7:41695/1‑148 (MQ=255)
                                                ttGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGc                                                                                                                          >  2:237098/1‑125 (MQ=255)
                                                    ccAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATgcgc                                                                                              >  3:325183/1‑149 (MQ=255)
                                                    ccAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATgcgc                                                                                              >  3:106450/1‑149 (MQ=255)
                                                    ccAGTTCGCGGAGCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACACTgcgc                                                                                              <  1:177028/149‑1 (MQ=255)
                                                        ttCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGcc                                                                                                                         <  7:296695/118‑1 (MQ=255)
                                                        ttCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGcc                                                                                                                         >  8:296695/1‑118 (MQ=255)
                                                            cgGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGtt                                                                                      >  3:309706/1‑149 (MQ=255)
                                                                  aGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAAc                                                                                <  8:41695/149‑1 (MQ=255)
                                                                   ggTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACg                                                                               >  1:216588/1‑149 (MQ=255)
                                                                              gcggcCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTa                                                                    <  4:20840/149‑1 (MQ=255)
                                                                                 gcCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCg                                                                 <  5:19227/149‑1 (MQ=255)
                                                                                    tcctGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTc                                                              >  6:352956/1‑149 (MQ=255)
                                                                                      ctGGTCTTCGCCATTGTCCACACCATTGATACCACCCAAAGGGAACCCGCGAATACCGACCCGGTAGGCGCTCCAGTTTGCAATGGCCCCATACACCGCCACCAGCAAAATGCGCGGACGGTTTTAAACGCTTCCGCCGTAGCGGTCtt                                                            >  1:203080/1‑149 (MQ=255)
                                                                                         gTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTcgc                                                         >  5:68863/1‑149 (MQ=255)
                                                                                            ttCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCCGTCTTCGCGCt                                                      >  8:59216/1‑149 (MQ=255)
                                                                                            ttCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTGTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCt                                                      >  6:251185/1‑149 (MQ=255)
                                                                                                  aTTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTgg                                                <  4:106450/149‑1 (MQ=255)
                                                                                                  aTTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTgg                                                <  4:325183/149‑1 (MQ=255)
                                                                                                   ttGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTc                                                                                   >  6:16089/1‑113 (MQ=255)
                                                                                                   ttGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTc                                                                                   <  5:16089/113‑1 (MQ=255)
                                                                                                      tCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACt                                            >  7:5044/1‑149 (MQ=255)
                                                                                                       ccACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTg                                           >  1:165112/1‑149 (MQ=255)
                                                                                                        cacaCCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTgg                                          >  3:156750/1‑149 (MQ=255)
                                                                                                          cacCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTgggg                                        <  2:165112/149‑1 (MQ=255)
                                                                                                            ccATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGt                                      <  5:352956/149‑1 (MQ=255)
                                                                                                               ttGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAg                                   >  3:61657/1‑149 (MQ=255)
                                                                                                               ttGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAg                                   >  5:217399/1‑149 (MQ=255)
                                                                                                                     ccaccCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTg                             <  4:156750/149‑1 (MQ=255)
                                                                                                                     ccaccCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTg                             <  6:68863/149‑1 (MQ=255)
                                                                                                                         ccAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTg                                     >  5:212605/1‑137 (MQ=255)
                                                                                                                         ccAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTg                                     <  6:212605/137‑1 (MQ=255)
                                                                                                                         ccAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTg                          >  1:333941/1‑148 (MQ=255)
                                                                                                                          cAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTgca                        <  4:309706/149‑1 (MQ=255)
                                                                                                                          cAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTg                          >  4:153478/1‑147 (MQ=255)
                                                                                                                                    acGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTGCAGCAGGTAAcc              <  7:59216/149‑1 (MQ=255)
                                                                                                                                        gAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGGCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTGCAGCAGGTAACCGTgg          <  8:5044/149‑1 (MQ=255)
                                                                                                                                          gcgccGACCGGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGt                                      >  6:61635/4‑119 (MQ=255)
                                                                                                                                          gcgccGACCGGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGt                                      >  2:22732/4‑119 (MQ=255)
                                                                                                                                          gcgccGACCGGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGt                                      <  5:61635/116‑1 (MQ=255)
                                                                                                                                          gcgccGACCGGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGt                                      <  1:22732/116‑1 (MQ=255)
                                                                                                                                              cGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTgg                            >  3:307814/1‑125 (MQ=255)
                                                                                                                                              cGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTgg                            <  4:307814/125‑1 (MQ=255)
                                                                                                                                                aCCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTGCAGCAGGTAACCGTGGGCGGCGTg  <  2:333941/149‑1 (MQ=255)
                                                                                                                                                 cccGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTg                             <  1:309986/121‑1 (MQ=255)
                                                                                                                                                 cccGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTGCAGCAGGTAACCGTGGGCGGCGTg  >  3:56486/1‑148 (MQ=255)
                                                                                                                                                 cccGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCAAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTg                             >  2:309986/1‑121 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCACCGGCCCAGTCCGGCTCGGACAGGTGCAGGTAGGCGAGGTTGCGCTTGGCCAGTTCGCGGATCAGGTACAGGCCGGCGGCCTCCTGGTCTTCGCCATTGTCCACACCATTGAAACCACCCAACGGGAACACGCGAATACCGACCCGGTCGGCGCTCCAGTTGGCAATGGCCGCATCCACCGCCTCCAGCACAATGCGCGCACGGTTTTCAACGCTGCCGCCGTAGCGGTCTTCGCGCTGGTTGGCACTGGGGGTGAGGAACTGGTGCAGCAGGTAACCGTGGGCGGCGTG  >  NC_002947/3596179‑3596471

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