Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,537,067 A→G 14.7% intergenic (+44/+59) 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,0670AG14.7% 56.0 / 3.0 27intergenic (+44/+59)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base G (2/2);  total (12/15)
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.80e-01

ATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTG  >  NC_002947/5536920‑5537213
                                                                                                                                                   |                                                                                                                                                  
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)
                          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)
                                                                                                                                                 tGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACAGGCTGCCAGCTTg  >  8:239690/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTG  >  NC_002947/5536920‑5537213

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: