Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,537,075 G→C 14.6% intergenic (+52/+51) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,537,0750GC14.6% 110.3 / 9.5 47intergenic (+52/+51)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base C (3/4);  total (22/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 = 9.09e-01

GGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGT  >  NC_002947/5536934‑5537220
                                                                                                                                              |                                                                                                                                                 
gggCAGCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCGGAATACGCGGGGGCACCCACCGACGCCtt                                                                                                                                            >  8:252032/1‑149 (MQ=255)
   cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGACGGTGAACCCACCGACGCATTcg                                                                                                                                          >  1:257753/1‑149 (MQ=255)
   cAGCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                         <  2:313188/149‑1 (MQ=255)
   cAGCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                         <  5:153625/149‑1 (MQ=255)
   cAGCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                         <  4:18873/149‑1 (MQ=255)
     gCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGTAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                       <  3:9038/149‑1 (MQ=25)
     gCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                       <  5:211390/149‑1 (MQ=25)
      ccGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCggg                                                                                                                                      <  8:84124/149‑1 (MQ=255)
       cGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAACCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                                     <  6:37553/149‑1 (MQ=255)
         gTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCGCGCATACGAGGGTTAACCCACCCACCCATTAGCGGGTg                                                                                                                                    >  3:337611/1‑148 (MQ=255)
            gACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                                 <  5:156190/148‑1 (MQ=21)
               cAGGCAGTGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGc                                                                                                                             <  2:7000/149‑1 (MQ=255)
                                 atcG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCt                                                                                                                         <  3:24169/135‑1 (MQ=255)
                                      tcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGCATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCGGGGTTCGCCGTCAGc                                                                                                       >  5:347581/1‑149 (MQ=255)
                                          cACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATATGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTGCGGGCGGGTTCGGCGTCAGCCgtt                                                                                                   >  6:83460/1‑149 (MQ=37)
                                                  gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAATAAGAGGGTGAACCCCCCGAACCAATACCCGGTGAACCCCCGCCAGCCCGGTTCGCCGGGTGCCGTTGATTaaat                                                                                           >  8:187790/1‑146 (MQ=255)
                                                    acgccccACCCCAAGCGCCACCTCCCGTTCGCTCCTTTATAAGTGCCCTGCTGGAGCGGGTTCACCCGCTAATACGAGGGTCAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGa                                                                                         <  3:237711/145‑1 (MQ=25)
                                                    aaGCCCCAGCCCAAGCGCGAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGa                                                                                         <  7:252032/149‑1 (MQ=255)
                                                    aaGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGATTACGAGGGTGAACCCACCGACGCATTCCCGGGTGCTCCCGCTCCTGCACGGTTCGCCGGCAGCCATTGATTACTTGa                                                                                         >  7:258953/1‑149 (MQ=255)
                                                        cccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGCATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCGCCAGCAGGGTTCGCCGTCAGCCTTTGATTACATGACGAt                                                                                     >  6:397433/1‑149 (MQ=255)
                                                         ccaccaCACGCGCCAACTGCCGTACGCTCATTGATATGGGGCCTCCAGGTGCGGGTACACCCGCTAATACTAGGGTTAACCCACCGACGCATGCGGGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATg                                                                                    <  7:148356/149‑1 (MQ=255)
                                                              acacGCGCGAACTGCCGTTCGCTCCTTGATAAGTGCCGTGCAGCAGCGGGTTCACCCGCGTATACGAGGGGTAACCCACCGACGCATTGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCtt                                                                               <  5:390270/149‑1 (MQ=255)
                                                                aaGCGCGAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc                                                                             <  5:397433/149‑1 (MQ=255)
                                                                 aGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCTGCCTTTGATTACTTGACGATGACCTTGcc                                                                            <  7:187790/149‑1 (MQ=255)
                                                                           gcggtgcGCTCCTTTATACGGGGCCGCCCGCCGCGGGTCCCCCCGCTAATACGCGGGCCACCCCACCCACGCATTCGCGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                  <  2:206459/143‑1 (MQ=255)
                                                                            ccGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                  <  1:131243/148‑1 (MQ=255)
                                                                             cGTGCGCTCCTTGATAAGGGGCCTGCCGGAGCGGGTTCCCCCGCTAATACGCGGGTGAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCGGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCga                                                                <  6:347581/149‑1 (MQ=255)
                                                                                                               gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg                              >  8:229950/1‑149 (MQ=255)
                                                                                                                   aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcac                          >  8:355349/1‑149 (MQ=25)
                                                                                                                    cccGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcac                          >  6:39464/1‑148 (MQ=25)
                                                                                                                    cccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                         >  1:269272/1‑149 (MQ=255)
                                                                                                                     ccGCGAATACGAGGGTGACCCCACCGACGCCTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGCGAAGGTGCAGAAGAAGATGTAGTcacca                        >  5:112244/1‑149 (MQ=255)
                                                                                                                     ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgt                                                           <  6:134124/114‑1 (MQ=255)
                                                                                                                     ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgt                                                           >  5:134124/1‑114 (MQ=255)
                                                                                                                     ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                        >  7:317449/1‑149 (MQ=255)
                                                                                                                     ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                        >  2:151459/1‑149 (MQ=255)
                                                                                                                       gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                               >  1:189692/1‑140 (MQ=255)
                                                                                                                       gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                               <  2:189692/140‑1 (MQ=255)
                                                                                                                       gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                      >  5:120084/1‑149 (MQ=255)
                                                                                                                          aaTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGc                   <  8:197497/149‑1 (MQ=255)
                                                                                                                                 gggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGc            >  6:168440/1‑149 (MQ=255)
                                                                                                                                  ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCt           >  6:166212/1‑149 (MQ=255)
                                                                                                                                   gTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt          >  8:177024/1‑149 (MQ=255)
                                                                                                                                       aCCCACCGCCGGATTCGCGGGGGAACCCGCTCCTGCTGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa      <  8:258953/149‑1 (MQ=255)
                                                                                                                                       aCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                    <  4:59798/119‑1 (MQ=255)
                                                                                                                                       aCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                    >  3:59798/1‑119 (MQ=255)
                                                                                                                                           aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt  <  4:337611/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCG‑TCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGT  >  NC_002947/5536934‑5537220

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