Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,467,909 A→G 7.8% S267G (AGC→GGC)  queF → NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,467,9090AG 7.8% 143.9 / 3.0 51S267G (AGC→GGC) queFNADPH‑dependent 7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base G (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

CCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCT  >  NC_002947/2467761‑2468052
                                                                                                                                                     |                                                                                                                                                
ccAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCa                                                                                                                                                  <  1:292541/149‑1 (MQ=255)
      tGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCg                                                                                                                                             >  6:334887/1‑148 (MQ=255)
         cGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTGATCGGCCCggtgc                                                                                                                                         <  7:15561/149‑4 (MQ=255)
          gACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAAGCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTAGGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATaa                                                                                                                                        <  3:95726/149‑1 (MQ=255)
                cATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACg                                                                                                                                  >  2:293497/1‑149 (MQ=255)
                    aGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCt                                                                                                                               <  5:234405/148‑1 (MQ=255)
                        gtgtgtGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCgg                                                                                                                                            >  3:17559/1‑131 (MQ=255)
                        gtgtgtGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCgg                                                                                                                                            <  4:17559/131‑1 (MQ=255)
                            gtgGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCAGATCAGCCCGGATAACAAACGGCTGGTCCGGCa                                                                                                                      >  7:224595/1‑149 (MQ=255)
                                aGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTaa                                                                                                                  >  3:168891/1‑149 (MQ=255)
                                 gCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCACCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTaa                                                                                                                  >  6:168823/1‑148 (MQ=255)
                                        tACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAggg                                                                                                           <  8:42489/148‑1 (MQ=255)
                                         aacTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAgggg                                                                                                          <  4:214003/146‑1 (MQ=255)
                                                 ttGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCAATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCACCCCGGATAACAAACGGCTGGTCCGGAAGTAATACAGAGCGCCCCATTg                                                                                                 >  4:280253/1‑149 (MQ=255)
                                                    aagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCg                                                                                               <  3:280253/148‑1 (MQ=255)
                                                         tCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGc                                                                                                                                >  1:119479/1‑110 (MQ=255)
                                                         tCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGc                                                                                                                                <  2:119479/110‑1 (MQ=255)
                                                           tgctgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGc                                                                                       >  7:147431/1‑149 (MQ=255)
                                                             ctgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCAAAGCGGccccc                                                                                         <  5:334887/145‑1 (MQ=255)
                                                              tgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAt                                                                                    >  8:130272/1‑149 (MQ=255)
                                                                 aGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCTCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTc                                                                                 <  1:41387/149‑1 (MQ=255)
                                                                 aGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCAAAGCGGccccc                                                                                         <  4:326373/141‑1 (MQ=255)
                                                                 aGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCAAAGCGGccccc                                                                                         >  3:326373/1‑141 (MQ=255)
                                                                      gAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATc                                                                            >  4:116532/1‑149 (MQ=255)
                                                                           ccTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACgg                                                                                                                                 <  3:160898/91‑1 (MQ=255)
                                                                           ccTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACgg                                                                                                                                 >  4:160898/1‑91 (MQ=255)
                                                                             tGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGc                                                                     >  1:228071/1‑149 (MQ=255)
                                                                                  gTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                 >  8:68900/1‑148 (MQ=255)
                                                                                   tCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                 >  3:150745/1‑147 (MQ=255)
                                                                                   tCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                 >  2:310016/1‑147 (MQ=255)
                                                                                      aTGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGc                                                            >  3:289216/1‑148 (MQ=255)
                                                                                       tGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTCGTTCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATCCAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCg                                                           >  1:102735/1‑149 (MQ=255)
                                                                                       tGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCg                                                           <  3:116532/148‑1 (MQ=255)
                                                                                         cgcgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAttt                                                                                  >  2:170924/1‑124 (MQ=255)
                                                                                         cgcgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAttt                                                                                  <  1:170924/124‑1 (MQ=255)
                                                                                           cgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTg                                                      <  7:130272/149‑1 (MQ=255)
                                                                                                        gcgGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCa                                         <  8:224595/149‑1 (MQ=255)
                                                                                                           gTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAggg                                      <  2:228071/149‑1 (MQ=255)
                                                                                                              ggCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTgg                                   <  8:147431/149‑1 (MQ=255)
                                                                                                              ggCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCACCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGTGTgg                                   >  2:365363/1‑149 (MQ=255)
                                                                                                              ggCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGGCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTgg                                   <  1:130993/149‑1 (MQ=255)
                                                                                                                  ggACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCa                                                               >  5:358902/1‑118 (MQ=255)
                                                                                                                  ggACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCa                                                               <  6:358902/118‑1 (MQ=255)
                                                                                                                    aCATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGt                               >  3:78689/1‑147 (MQ=255)
                                                                                                                        cAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGt                                     >  5:63220/1‑137 (MQ=255)
                                                                                                                         aaCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCg                        >  7:155944/1‑149 (MQ=255)
                                                                                                                         aaCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCg                        <  8:204135/149‑1 (MQ=255)
                                                                                                                                    cagcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGc             <  1:310016/149‑1 (MQ=255)
                                                                                                                                        tccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTc          >  1:189253/2‑148 (MQ=255)
                                                                                                                                         ccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCa                                         <  1:78301/116‑1 (MQ=255)
                                                                                                                                         ccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCa                                         >  2:78301/1‑116 (MQ=255)
                                                                                                                                                 gATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCATCCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCt  <  1:293497/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCT  >  NC_002947/2467761‑2468052

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: