Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 718,460 C→T 8.4% intergenic (+30/+47) 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,4600CT 8.4% 127.2 / 2.4 49intergenic (+30/+47)thiO/PP_0613FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein
Reads supporting (aligned to +/- strand):  ref base C (14/29);  new base T (2/2);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCAC  >  NC_002947/718315‑718605
                                                                                                                                                 |                                                                                                                                                   
aTCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc                                                                                                                                                   >  2:296376/1‑148 (MQ=255)
 tCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc                                                                                                                                                   >  7:65360/1‑147 (MQ=255)
  cGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAg                                                                                                                                               >  1:82666/1‑149 (MQ=255)
     aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCgg                                                                                                                                            <  3:135674/149‑1 (MQ=255)
         ggCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAACCCGCGCCGAAGCGGGGCGGGtt                                                                                                                                         >  7:142363/1‑149 (MQ=255)
          gCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTa                                                                                                                                       <  1:302548/149‑1 (MQ=255)
           cTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑c                                                                                                                                     <  1:296376/149‑1 (MQ=255)
              gTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑ccc                                                                                                                                   <  5:59977/148‑1 (MQ=255)
                           cTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGttt                                                                                                                                        <  4:106776/131‑1 (MQ=255)
                           cTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGttt                                                                                                                                        >  3:106776/1‑131 (MQ=255)
                                         tttACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCGGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGa                                                                                                       <  7:315460/149‑1 (MQ=255)
                                         tttACCGATTTGCTGACAGGCGCCGACCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTCCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTAA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGa                                                                                                       <  1:123977/149‑1 (MQ=39)
                                              cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCTTATGCGCCAGAAGGGCGTTTGGGCTGAGTACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGTTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  7:28754/149‑1 (MQ=255)
                                              cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCTTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTTCGCGTTTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  1:95891/149‑1 (MQ=255)
                                              cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  8:65360/149‑1 (MQ=255)
                                              cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  5:186243/149‑1 (MQ=255)
                                              cGATTTGCTCACAGGGGCCCAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTTCGCGGGTAACCCCGCTCCAACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  7:195812/149‑1 (MQ=255)
                                               gATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGGTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGAtttacg                                                                                                 >  6:183808/1‑146 (MQ=255)
                                                          aGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCGCCTACGC‑GAGCGGGTTTT‑CCCGCGAAAGGGCCGCAGTAGCCAATAAAGATTTAATCCAGTAATCgg                                                                                      >  8:75121/1‑149 (MQ=255)
                                                           ggggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCACCTACGC‑GAGCGGGGTTA‑CCCCCGAAAAGGCCGCACAACCCAAAAAAGATTGAAGCCAGTAAACGGt                                                                                     >  3:75685/1‑149 (MQ=255)
                                                             ggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCTACAACGC‑GAGCGGGTTTA‑ACCGCGAAAGGGCCGCAGAAGCCAAAACAGATTTAAGCCAGTAAAAGGTGc                                                                                   >  6:144671/1‑149 (MQ=255)
                                                                  ggccgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑c                                                                                                                                     <  2:89735/93‑1 (MQ=25)
                                                                  ggccgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑c                                                                                                                                     >  1:89735/2‑94 (MQ=25)
                                                                     cgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGt                                                                           <  5:183808/149‑1 (MQ=255)
                                                                       aTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTCCGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGc                                                                         <  2:260830/149‑1 (MQ=255)
                                                                             gACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCAACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATAAAGATTTGAGCCAGTAATAGGTGGTTGCAGGTGCCGGcgc                                                                   >  6:119289/1‑147 (MQ=255)
                                                                             gACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCGCCGAAGC‑GAGCGGGGTTT‑ACCGCGAAAGGGGCGCCGAAGGCCATAAAGAGTTCAGCCAGTAATAGGTGGTTGCAGGTggcggc                                                                     >  4:79136/1‑147 (MQ=255)
                                                                                cgggggcCTTATGCGCCAGAAGGGCTTTTGGGCGGTGGACCTCTTCGCGGGTAAACCCGCTCCTCCGC‑TAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGa                                                                <  3:228175/144‑1 (MQ=38)
                                                                                 gggCGCCTTATCCCCCAGACGCGGGTTTGCGCGCAGCACGGCTTCGCGGTTAAACCCGCCCCACCGC‑GAGCGGTTTAA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGAt                                                               <  3:29266/148‑1 (MQ=255)
                                                                                     gcCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                           <  7:247735/149‑1 (MQ=255)
                                                                                         tATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCTACGC‑GAGCGGGTTTA‑ACCGCGAAAGGGCCGCAGAAGCCAAAACAGATTTAAGCCAGTAATAGGTGCTTGCAGGTGCCGGCTCGGCTAAtgctgc                                                       >  2:228366/1‑149 (MQ=255)
                                                                                                atactggggGTTTGGGCTGCGAACGCCTGCG‑GGGTAAAACCGCCACTACCC‑CCGGGGGTTTACCCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATaca                                                <  1:367396/144‑1 (MQ=255)
                                                                                                     gggCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg                                           <  2:82666/149‑1 (MQ=255)
                                                                                                              ggCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgc                                  <  1:228366/149‑1 (MQ=255)
                                                                                                                cTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgcag                                <  5:144671/149‑1 (MQ=255)
                                                                                                                       ccgctgcgcgGGTAAACCCGCCCCACCCC‑CAGCCGGATTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc                         <  8:142363/143‑1 (MQ=255)
                                                                                                                       ccTCTTCGCGGGTAACCCCGCTCCAACGC‑TCGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc                         <  3:79136/149‑1 (MQ=255)
                                                                                                                         tctTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGTTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  3:125841/74‑1 (MQ=255)
                                                                                                                              gcgGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAg                  >  3:211197/1‑149 (MQ=255)
                                                                                                                              gagtgaaaaCACGCTCATACGC‑CAGCGGTTTTA‑CCCGCGAAAGTGCCGCAGAATCCAATACATATTTCAGCCATTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAg                  <  6:117778/143‑1 (MQ=255)
                                                                                                                                 ggTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca                   >  4:208397/1‑145 (MQ=39)
                                                                                                                                 ggTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca                   <  3:208397/145‑1 (MQ=39)
                                                                                                                                    aaaCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                           <  4:141095/102‑1 (MQ=255)
                                                                                                                                    aaaCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                           >  3:141095/1‑102 (MQ=255)
                                                                                                                                      aCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGc                                          <  6:68743/117‑1 (MQ=255)
                                                                                                                                      aCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGc                                          >  5:68743/1‑117 (MQ=255)
                                                                                                                                        ccGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc        >  6:227171/1‑149 (MQ=255)
                                                                                                                                        ccGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc        >  2:272339/1‑149 (MQ=255)
                                                                                                                                        ccGCTCCTACCC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc        <  7:75121/149‑1 (MQ=255)
                                                                                                                                           cTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGc     <  5:119289/149‑1 (MQ=255)
                                                                                                                                              cTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGccac  >  8:80549/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCAC  >  NC_002947/718315‑718605

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: ‑