Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,537,061 A→T 24.1% intergenic (+38/+65) 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,0610AT24.1% 46.1 / 12.5 35intergenic (+38/+65)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base T (2/6);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

ATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGT  >  NC_002947/5536920‑5537193
                                                                                                                                             |                                                                                                                                    
aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCTCGTATTCGCGGGTGaa                                                                                                                               <  5:75489/149‑1 (MQ=255)
  gCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTgggtt                                                                                                                             >  5:186697/1‑145 (MQ=255)
     ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTgggttcac                                                                                                                          >  3:152909/1‑142 (MQ=38)
     ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATACCAGGGTGAACCCAc                                                                                                                          >  7:282682/1‑149 (MQ=255)
        cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGa                                                                                                                       <  8:222912/149‑1 (MQ=255)
         tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                      >  3:85555/1‑149 (MQ=255)
                 cAGCCGGTTGACCAGTCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGTACCCACCGACGCATTcgc                                                                                                              <  6:147508/149‑1 (MQ=255)
                 cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                              <  4:85555/149‑1 (MQ=255)
                 cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                              <  2:36592/149‑1 (MQ=255)
                   gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                            <  3:19668/149‑1 (MQ=25)
                   gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                            <  4:152909/149‑1 (MQ=25)
                   gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                            <  6:186697/149‑1 (MQ=25)
                   gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                            <  5:192460/149‑1 (MQ=25)
                   gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                            <  3:177887/149‑1 (MQ=255)
                             cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGc                                                                                                  <  6:115251/149‑1 (MQ=255)
                               ggCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTc                                                                                                <  1:20053/149‑1 (MQ=255)
                                           aatcttcGTCGCCACCTACCCCAAGACCCAGCCAATGCGCTAACTGCGTTTCGTTCTTTTATAAGTGGCGTGCTGTCGGGGGTTCACCCGCTAATACGTGGGTTAACCCACCGACTCATTCGGGGGTTAACCCGCTCCTGCAGGGTTc                                                                                     <  8:282682/147‑1 (MQ=255)
                                                 cgtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                    <  1:56260/95‑1 (MQ=255)
                                                 cgtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                    >  2:56260/1‑95 (MQ=255)
                                                       cACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCttt                                                                        >  6:170831/1‑149 (MQ=255)
                                                          ctcccGCAACACCCACCACAAGCGCCACCTCCCGTCCGCGCCTTGATAAGTGCCCTGCAGGAGCGGGTACACCCGCGAATACGAGGGTGAACCCACCGACCCATCCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAt                                                                     <  5:151430/146‑1 (MQ=37)
                                                          cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCGCGAATAAGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCGCCTGGAGGGTTAGCCGGCAGCCGTTGAt                                                                     >  6:273961/1‑149 (MQ=38)
                                                            aCCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTAGCGGGTGAACCCGCGACTGCAGGGTTAGCCGTCAGCCATTGAtt                                                                    >  4:228459/1‑148 (MQ=39)
                                                               gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                <  5:170831/149‑1 (MQ=255)
                                                                cAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                  <  5:184182/82‑1 (MQ=255)
                                                                cAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                  >  6:184182/1‑82 (MQ=255)
                                                                        aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGAc                                                       <  5:273961/149‑1 (MQ=255)
                                                                              ggagcgAACGGCAGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCGCGAATACCAGGGTGAACCCACCGACGCATTAGCGGGTGAACCCGCTCCTGGAGGGTTCGCCGTGAGCCTGTGATTTCGTGACGCTGACCGTggc                                                 >  4:74779/4‑147 (MQ=255)
                                                                                         ccGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                              <  5:197488/61‑1 (MQ=255)
                                                                                         ccGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                              >  6:197488/1‑61 (MQ=255)
                                                                                         ccGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                       <  3:228459/148‑1 (MQ=255)
                                                                                                   cTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCCCGAATAACAGGGTTAACCCACCCACGCCATAGCGGGTGGACCCCCGCCGGCACGGGGGGCCGACCGCCGTTTGTTAAATTACGATGACCATTGCGTTCATAATCCAGAAGTTGcc                            >  8:4353/1‑149 (MQ=255)
                                                                                                               ccTCCAGGAGGGGTTTCACCCGCTAAAAGGTGGGTCAACCCACGTACGCTTTGGGGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                <  3:298893/149‑1 (MQ=255)
                                                                                                                cTGCAGGAGCGGGTTCACCGGCGAATACTCGGGTCAACCCACCGACGCATTCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGACAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCag               <  7:4353/149‑1 (MQ=38)
                                                                                                                             ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt  >  7:222771/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                    
ATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGT  >  NC_002947/5536920‑5537193

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: