Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 718,457 C→G 9.8% intergenic (+27/+50) thiO → / ← PP_0613 FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947718,4570CG 9.8% 98.5 / 5.0 47intergenic (+27/+50)thiO/PP_0613FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base G (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-02

ATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGC  >  NC_002947/718312‑718602
                                                                                                                                                  |                                                                                                                                                  
aTTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCAGGTAAACCCGCTCCTAc                                                                                                                                                 >  3:167338/1‑149 (MQ=255)
   aTCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAACCGCGCCGAagcg                                                                                                                                              >  1:286545/1‑149 (MQ=255)
    tCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc                                                                                                                                               <  1:238548/148‑1 (MQ=255)
       cAATGGGCTGGTGCTGGCGCCGGCTTCGTGGCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCg                                                                                                                                          >  1:25131/1‑149 (MQ=255)
        aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                         <  5:90960/149‑1 (MQ=255)
        aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                         >  1:59971/1‑149 (MQ=255)
            ggCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGtt                                                                                                                                      <  2:59971/148‑1 (MQ=255)
            ggCTGGTGCTGGCGCCGGCTGCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGtt                                                                                                                                      <  5:266900/148‑1 (MQ=255)
              cTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTAc                                                                                                                                   <  7:236197/149‑1 (MQ=39)
                 gTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTAcc                                                                                                                                  >  5:129350/1‑147 (MQ=39)
                 gTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTAcc                                                                                                                                  <  6:129350/147‑1 (MQ=39)
                       gcgcCGGCTTCGGGCCTGCTGTTTACCGTTTTCCTTCCCGGGGCCGAGCGTATAAGCGCCCCGGCGCCTTATGCCCCAGACGGGGGTTTGGGCTGAGGACCTCTTCGCGGTTAAACCCGCTCCAACGCGAGCGGGTTTACCCGCGAAAg                                                                                                                         <  2:255381/149‑1 (MQ=255)
                         gcCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAAACGACACGACACGAGGGGGTTTTCCCGAGaaaacg                                                                                                                       >  8:289717/1‑146 (MQ=255)
                                  gTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCAGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGaa                                                                                                               <  4:167338/148‑1 (MQ=255)
                                   tGCCTGCTTTTTACCGATTTGCTGCCCGGGGCCGCGCGTATTAGCGACCCGGCGCGGTATGCGCCAGATGGGCGTTTGGGCTGAGCACCTCTTCGCGGTTAACCCCGCTCCTACGCGCGCGGTTTTACCCGCGAAAGGGCCGCAGAAGc                                                                                                             <  2:286545/149‑1 (MQ=255)
                                    gCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCTCCTAAGCGAGCGGGTTTAACCGCGAAAGGGCCGCAGAAGc                                                                                                             >  7:172468/1‑148 (MQ=255)
                                             ttACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGGTCCTACGCGAGCGGGTTTTCCCGCGCAACGGGCCCCGAAGCCAATACAggt                                                                                                   >  5:299278/1‑147 (MQ=255)
                                             ttACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCTACGCGCGCGGGTTTTACCGCGAAAGGGCCGCAGAAGCCAATACAGAt                                                                                                   >  7:5518/1‑149 (MQ=255)
                                              tACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACATCTTCGCGGGTAAAACCGCGCCTACGCGGGCCGGGTTTCCCGCGAAAGTGCCGCAGAACCCAATAAAGAtt                                                                                                  >  7:145098/1‑149 (MQ=255)
                                                 ggATTTGCTGACAGGGGCCGAGCCTATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGa                                                                                                                <  2:223579/131‑1 (MQ=255)
                                                   aTTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCAACAAACCGAGCGGGTTAAACCCCGAAACGGCCCTATAAGACAATAAAGAtttactg                                                                                             >  6:167922/1‑145 (MQ=255)
                                                        cTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAACCGCGACGACGCGAGCGGGGTTTACCGGGAAAGGGGCGCAGAAGCCAATAAAAATTTGAGGCAGt                                                                                         >  4:258153/1‑148 (MQ=255)
                                                            cAGGGTCCGAGCCGATTATCTACCGGGGGCCGTATGCGCAAGATGGTGGTTTGGGGTTAGGACCTCTTCTCGGGTAAACCCGCTCCTACGCTTGGGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCg                                                                                    <  8:172468/149‑1 (MQ=37)
                                                               gggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCGCCAAAGCGAGCGGGTTTTCCCGCGAAAGGGCCCCCGAAGCCAATAAAAATTTAAGCCAGTAATAGGTg                                                                                 >  2:208743/1‑149 (MQ=255)
                                                                         gaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTg                                                                       >  1:29273/1‑149 (MQ=255)
                                                                                  cccGGCGCGGTATGCGCCTGAAGGGGGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGTGCGGTTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCgg                                                              <  8:5518/149‑1 (MQ=39)
                                                                                        gcCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGTAAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                        <  3:258153/149‑1 (MQ=255)
                                                                                         ccGTATGCGCCAGACGGGCGTTTGGGCTGAGCACCTCATCGCGGGAAAACCCGCCCCAACGCGCGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCaa                                                                                                          <  5:9950/98‑1 (MQ=37)
                                                                                           gTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCGGCAAAACGAGCGGGGGTTACCGCGAAAAGGCCCCAGAAACCAAAAAAAATTTAAGGAAGGAATAGGTGCTTTCAGGTGGCCGCGGGGATAAtta                                                       >  8:290266/1‑145 (MQ=255)
                                                                                           gTAAGCGCCAGAACGGAGTTTCGGCTGAGGACCTCTTCGCGGGTAAAACCACACCAACACCACCCGGGTTACAAGCGAAACGGCCGCAGAAGCCAAAACAGAAATCAACCAGAAATCAATACTTTCATGTTCCGACTAGGATTAtgctg                                                     >  8:131978/1‑149 (MQ=255)
                                                                                                        aGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAACCGCTACAAAGCGAGCGGGTTTACCCGCCAAAGGGGCGGAGAAGCCAAAACAGa                                                                                                    >  4:175835/1‑90 (MQ=255)
                                                                                                        aGGGCGTTTGGGCTGAGGACCTCTTCGCGGGAAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGa                                                                                                    <  3:175835/90‑1 (MQ=255)
                                                                                                         gggCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg                                        <  8:242458/149‑1 (MQ=255)
                                                                                                         gggCGTTGGGGCGGCGCCCGTCTTCGCGGGTAAACCCGCTCCTACGCCAGGGGTTTTACCCGCGAAAGGGCCGTAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg                                        <  5:167922/149‑1 (MQ=255)
                                                                                                           gCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCACCAACGCGAGCGGGGTTTCCCCCGCAAAGGGCCCAGAAACCCAAAAAGAATTAAGCCACTAATACGTTCTTTCAAGTGGCGGCGGGGATAAAACAGGTCAAACAccgatg                                       >  4:279323/1‑144 (MQ=255)
                                                                                                                    cTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgcag                             <  2:29273/149‑1 (MQ=255)
                                                                                                                           ccTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTAACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc                      <  1:208743/149‑1 (MQ=255)
                                                                                                                           ccTCTTCGCGGGTAAACCCGCGCCTACGCGTGCGGTTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTg                       <  7:131978/148‑1 (MQ=255)
                                                                                                                              tttCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtca                   >  6:10614/2‑149 (MQ=38)
                                                                                                                               ttCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAAAACAGATTTCAGCCAGTAATCGGTGCTTGCACGTGCCGGCACGGATAATGCTGCTCATACACCTTGCACCCCCGCACCACTTGCtca                   >  4:260497/1‑148 (MQ=39)
                                                                                                                                  gcgGGTAAACCCTCTCTTACGCGTGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca                <  7:289717/148‑1 (MQ=255)
                                                                                                                                        aaaCCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAgc         >  2:274518/1‑149 (MQ=255)
                                                                                                                                          aCCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCTGCAAAgcgc       >  3:163783/1‑149 (MQ=255)
                                                                                                                                          aCCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCACAgcgc       >  4:184557/1‑149 (MQ=255)
                                                                                                                                            ccGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc     >  3:226565/1‑149 (MQ=255)
                                                                                                                                             cGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCg    <  1:274518/149‑1 (MQ=255)
                                                                                                                                               cTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGc  >  3:3704/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGC  >  NC_002947/718312‑718602

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