Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,719,331 A→G 7.3% Q153R (CAG→CGG)  nfuA → Fe/S biogenesis protein NfuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,719,3310AG 7.3% 156.3 / 2.5 55Q153R (CAG→CGG) nfuAFe/S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

AACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCCTACTACTGAGTTGAAGG  >  NC_002947/2719186‑2719466
                                                                                                                                                 |                                                                                                                                       
aaCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAgg                                                                                                                                       >  8:274485/1‑148 (MQ=255)
aaCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAgg                                                                                                                                       <  7:274485/148‑1 (MQ=255)
         aGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGc                                                                                                                               >  6:43252/1‑147 (MQ=255)
                tCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                                                                                      <  6:168012/149‑1 (MQ=255)
                   acgaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGt                                                                                                                    >  4:157510/1‑148 (MQ=255)
                     gaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAc                                                                                                                 >  7:169009/1‑149 (MQ=255)
                      aacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAcc                                                                                                                <  1:183053/149‑1 (MQ=255)
                           aTCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATCGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACCGCATTGCCGTACTCCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                                                                                           >  5:2934/1‑149 (MQ=255)
                            tCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                                                                                           <  8:295337/148‑1 (MQ=255)
                            tCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                                                                                           <  8:169009/148‑1 (MQ=255)
                                ctactaCCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaag                                                                                                      <  7:160047/149‑1 (MQ=255)
                                 tactacCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCATGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGGTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACACTGaaggaagg                                                                                                     >  5:273960/1‑149 (MQ=255)
                                      cctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATc                                                                                                 <  5:43252/148‑1 (MQ=255)
                                          cAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCTGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                                                                                             <  4:345685/148‑1 (MQ=255)
                                          aaGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCTGGTGAGCCTGGTGTACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                                                                                             <  6:2934/147‑1 (MQ=255)
                                             aCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGa                                                                                               >  2:73697/1‑143 (MQ=255)
                                             aCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGa                                                                                               <  1:73697/143‑1 (MQ=255)
                                             aCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAc                                                                                         >  1:311143/1‑149 (MQ=255)
                                             aCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAc                                                                                         >  2:129850/1‑149 (MQ=255)
                                               cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCGGGGTTGCgcc                                                                                                                               >  5:303659/1‑107 (MQ=255)
                                               cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCGGGGTTGCGGc                                                                                                                               <  6:303659/109‑1 (MQ=255)
                                               cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAccc                                                                                                               <  3:310921/125‑1 (MQ=255)
                                               cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAccc                                                                                                               >  4:310921/1‑125 (MQ=255)
                                                        ccccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGCGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGc                                                                              >  6:161425/1‑149 (MQ=255)
                                                         cccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaag                                                                                                          <  8:191273/120‑1 (MQ=255)
                                                         cccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaag                                                                                                          >  7:191273/1‑120 (MQ=255)
                                                         cccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCg                                                                             >  3:56936/1‑149 (MQ=255)
                                                          ccGGCCTGGCCAGCCTTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCg                                                                             >  3:18194/1‑148 (MQ=255)
                                                                tGGCCAGCCATGGCGGCCAGGTGAGACTGGTGGAAGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCg                                                                      >  6:157421/1‑149 (MQ=255)
                                                                  gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATc                                                                                                 >  2:160016/1‑120 (MQ=255)
                                                                  gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATc                                                                                                 <  1:160016/120‑1 (MQ=255)
                                                                  gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa                                                                    <  1:129850/149‑1 (MQ=255)
                                                                   ccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAg                                                                   >  8:279154/1‑149 (MQ=255)
                                                                                 cAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGCGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTgcg                                                     <  5:161425/149‑1 (MQ=255)
                                                                                       aGCCTGGTGGACGTGGTCGATGACGGCAGTGCCGTACTGCAGTTGGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGt                                               <  5:196733/149‑1 (MQ=255)
                                                                                            ggtggACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGtgact                                           >  6:70306/1‑148 (MQ=255)
                                                                                               ggACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGAcc                                       <  4:18194/149‑1 (MQ=255)
                                                                                                gACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGAcc                                       <  3:157510/148‑1 (MQ=255)
                                                                                                 aCGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACca                                      >  5:155586/1‑148 (MQ=255)
                                                                                                   gTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAg                                   >  5:298061/1‑149 (MQ=255)
                                                                                                      gTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGc                                  >  2:194143/1‑147 (MQ=255)
                                                                                                       tCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGc                                  >  6:306126/1‑146 (MQ=255)
                                                                                                        cGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCag                               >  7:132342/1‑148 (MQ=255)
                                                                                                          aTGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCagaaa                            <  7:279154/149‑1 (MQ=255)
                                                                                                               ggCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGaaaa                       <  4:56936/149‑1 (MQ=255)
                                                                                                                      ccGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCagaa                             >  2:160068/1‑136 (MQ=255)
                                                                                                                      ccGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCagaa                             <  1:160068/136‑1 (MQ=255)
                                                                                                                       cGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCctac                >  6:310903/1‑148 (MQ=255)
                                                                                                                              cAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCcta                 <  1:285215/140‑1 (MQ=255)
                                                                                                                              cAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCcta                 >  2:285215/1‑140 (MQ=255)
                                                                                                                               aGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCCTACTACTGAGtt       <  5:70306/149‑1 (MQ=255)
                                                                                                                                 ttCGGTGGTGGCTGCCGGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCCtactac             <  8:261074/141‑1 (MQ=255)
                                                                                                                                 ttCGGTGGTGGCTGCCGGGGTGGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCCtactac             >  7:261074/1‑141 (MQ=255)
                                                                                                                                   cGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCCTACTACTGAGTTGaa    <  8:132342/148‑1 (MQ=255)
                                                                                                                                    ggtggtGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCCTACTACTGAGTTGAAgg  <  3:233900/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
AACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGCGTGCGCGACGTGACTGACCACAGCCAGAAAGAAAACGCCTACTACTGAGTTGAAGG  >  NC_002947/2719186‑2719466

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