Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,537,058 A→G 11.8% intergenic (+35/+68) 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,0580AG11.8% 83.3 / 3.8 34intergenic (+35/+68)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (11/19);  new base G (2/2);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

CAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGC  >  NC_002947/5536913‑5537203
                                                                                                                                                 |                                                                                                                                                 
cAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa                                                                                                                                                <  8:150526/149‑1 (MQ=255)
 aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCgtc                                                                                                                                               >  3:174699/1‑147 (MQ=255)
  gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTTGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                              <  5:46342/149‑1 (MQ=255)
  gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                              <  4:33295/149‑1 (MQ=255)
   aTCGTTGGCTTCCTGCATTGGCTGCCGCTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACGCAGCACTAGCTCGAACTGCCGTTGGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCTCCCGCGAATACGAgg                                                                                                                                              <  8:302447/148‑1 (MQ=255)
       aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                         <  8:224503/149‑1 (MQ=255)
              ccTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAc                                                                                                                                    <  4:282604/147‑1 (MQ=255)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgcg                                                                                                                       <  6:74151/149‑1 (MQ=255)
                          gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                      <  4:174699/149‑1 (MQ=25)
                          gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                      <  6:129650/149‑1 (MQ=255)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                              <  1:152438/121‑1 (MQ=255)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                              >  2:152438/1‑121 (MQ=255)
                                 gACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                <  5:251839/148‑1 (MQ=21)
                                  aCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAcc                                                                                                               <  2:38990/148‑1 (MQ=255)
                                  aCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTAGCGGGTGAAccc                                                                                                              >  1:114183/1‑149 (MQ=255)
                                                             ccACCTACCCAACTACCCACCACACGCGCTACATCCCTTTCTCTCCTTTAAAAGTGGCCTCCAGCCGCGGGTTCACCCGCTAATACTAGGGTTAACCAACCCACGATTTCGCGGGTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCtt                                                                                   <  6:265069/149‑1 (MQ=255)
                                                                    ccccacgACCCCGCACACGCGCGAACTGCCTTTCGCTCTTTTATAAGGGGCCTCCCGGCGCGGGTTCACCCGCTAATAGTCGGGTCAACCCACCGACGCATTCGCGGTTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAc                                                                            <  2:114183/143‑1 (MQ=255)
                                                                     cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTACAGGGTTCGCCGTCa                                                                                        <  1:108040/136‑1 (MQ=255)
                                                                     cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATACGAGGGTGAAACCACCCACGCATTAGCGGGTGAACCCGCTCCTACAGGGGTCGcccgca                                                                                        >  2:108040/1‑132 (MQ=37)
                                                                      gcaacccccAGCACAAGCGCCAACTGCCGTTCGGCCCTTGATAAGTTGGCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                          <  4:226471/143‑1 (MQ=38)
                                                                      gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCCAAAACAAAGGTGAACCCACCCACGCATTAGCCGGTGAACCCCCACCAGCAGGGTTAGCCGTAAGCCATTGATTaaa                                                                           >  2:64084/1‑146 (MQ=255)
                                                                            cccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATACGAGGGTGAACCCCCCGACGCATTCGCGGGTGAACCCGCTACTGCAGGGTTAGCCGGCAGCCTTTGATTAATTGACGAt                                                                    >  6:88631/1‑149 (MQ=38)
                                                                             ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGGATACGAGGGTGGACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCCGGGGTCGCCGTCAGCCGTTGATTTCTTGACGCTg                                                                   >  1:13951/1‑149 (MQ=38)
                                                                                  acaAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCtt                                                              >  1:62734/1‑149 (MQ=255)
                                                                                         ccaaCTGCCGTTCGCTCATTTATAAGTGCCCTGCTGGAGCGGGTTCACCCGCTAATACGAGGGTCAACCCACCGACGCATTCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                       <  7:308759/147‑1 (MQ=255)
                                                                                          gAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCGCGAATAAGAGGGTGAACCCACCCACGCATTAGCGGGTGAACCCGCACCGGGAGGGTTCGCCGGCAGCCGTTGATTTCTTGGCGGTGACCGTGCCTTtc                                                       >  1:220156/1‑148 (MQ=255)
                                                                                          gAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                      <  7:273654/149‑1 (MQ=255)
                                                                                            aCTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcatc                                                    <  1:64084/149‑1 (MQ=255)
                                                                                                                                        cccGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc        >  8:99823/1‑149 (MQ=25)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaa                       >  4:265853/1‑133 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaa                       <  3:265853/133‑1 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca       >  5:150440/1‑149 (MQ=255)
                                                                                                                                           gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                    >  8:294469/1‑86 (MQ=255)
                                                                                                                                              aaTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGc  <  2:13951/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGC  >  NC_002947/5536913‑5537203

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: