Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,928,420 C→T 8.3% intergenic (+96/‑44) peaE → / → PP_3464 phenylacetaldehyde dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,928,4200CT 8.3% 165.1 / 4.2 60intergenic (+96/‑44)peaE/PP_3464phenylacetaldehyde dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (3/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

AATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCACCACCCTCAACGCCGACCT  >  NC_002947/3928302‑3928567
                                                                                                                      |                                                                                                                                                   
aaTCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAAc                                                                                                                        <  1:281126/148‑1 (MQ=255)
     gtgtTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGTAGCGGGTTAACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGa                                                                                                                  <  6:186212/149‑1 (MQ=255)
      tgtTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATgg                                                                                                                   >  1:262185/1‑147 (MQ=255)
      tgtTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATgg                                                                                                                   <  2:262185/147‑1 (MQ=255)
        ttGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTAGCGGGTAAAACCGCTACAACACCACGTATAACATGCCCTGTAAGAGCGGGTTAACCCGCGAAGAGGCCCGCACAGGCAAAAAAAAAAAAACAAAAAAAAACATTGatc                                                                                                                >  2:291450/1‑146 (MQ=255)
              aaGCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCCCGCCGAAACCAGGTATAACATGCCCTGTAGGAGCGGGGTAACCCGCGAAGAGGCCGGCCCAGGCAACAAAAAAACAACAAAAAAAACCCTGGAGCCAcccc                                                                                                          >  2:286160/1‑148 (MQ=37)
                gCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCCCATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGa                                                                                                                  <  8:167215/138‑1 (MQ=255)
                gCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATAACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGa                                                                                                                  >  7:167215/1‑138 (MQ=255)
                 ctctGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGa                                                                                                      <  1:286160/149‑1 (MQ=255)
                      aTCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGcac                                                                                                 >  3:250633/1‑149 (MQ=255)
                         cTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGAGATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAAC‑AAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCacacac                                                                                             >  1:85019/1‑149 (MQ=255)
                          tttttCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTTTCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAaacaa                                                                                                                                  >  3:211386/1‑112 (MQ=255)
                               ccTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGc                                                                                         <  5:283054/148‑1 (MQ=255)
                               ccTGTACCGGCCTCTTCGCGGGTAAACCCCCGGCTACACCAGGGATGACATGCCCGGTAGGAGCGGGGTGACCCGCGAAGAGGCCGGCCCAGGCAACAAAAAAACAACCAAAAAAACCAAGGAGGAACCCCAAGAGGcacacccacccc                                                                                        >  6:308621/1‑142 (MQ=25)
                                   tACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTAGCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGc                                                                                     <  3:52864/148‑1 (MQ=255)
                                    aCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCTGGTATCACATGCCCTGTAGGCGCGGTTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGca                                                                                    <  5:308621/148‑1 (MQ=255)
                                      CGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACAt                                                                                 <  4:218568‑M1/149‑1 (MQ=255)
                                       GGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACAtt                                                                                >  4:250633‑M1/1‑149 (MQ=255)
                                           TCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGcc                                                                            <  1:319249‑M1/149‑1 (MQ=255)
                                           TCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGcc                                                                            <  6:310297‑M1/149‑1 (MQ=255)
                                             TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGcac                                                                                                 <  5:14303‑M1/126‑1 (MQ=255)
                                             TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGcac                                                                                                 >  6:14303‑M1/1‑126 (MQ=255)
                                             TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCct                                                                          <  3:309260‑M1/149‑1 (MQ=255)
                                              gaGCGGGTTTACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCctg                                                                         <  6:38587‑M1/147‑1 (MQ=255)
                                                  gggtgaacccgctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGt                                                                     >  7:164825‑M2/69‑149 (MQ=255)
                                                         cccgctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTg                                                                               >  6:231082‑M2/62‑132 (MQ=255)
                                                         cccgctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTg                                                                               <  5:231082‑M2/71‑1 (MQ=255)
                                                          ccgctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggTACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAAcc                                                                <  6:1715‑M2/86‑1 (MQ=255)
                                                          ccgctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggTACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAAcc                                                                >  5:1715‑M2/61‑146 (MQ=255)
                                                          ccgctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggTACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGcc                                                             >  7:165465‑M2/61‑149 (MQ=255)
                                                           cgctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggTACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGccc                                                            >  3:351316‑M2/60‑149 (MQ=255)
                                                            gctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGtt                                                                    >  6:250952‑M2/59‑140 (MQ=255)
                                                            gctcctacagggcatgtgatacctggtgtaggagcgggtttacccgcgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGtt                                                                    <  5:250952‑M2/82‑1 (MQ=255)
                                                                 tACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCt                                                                                       <  6:77617/116‑1 (MQ=255)
                                                                 tACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCt                                                                                       >  5:77617/1‑116 (MQ=255)
                                                                      cAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCaa                                                 <  5:310297/149‑1 (MQ=255)
                                                                         gTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGcc                                              >  7:85843/1‑149 (MQ=255)
                                                                              acaTGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTAc                                           >  8:23363/1‑147 (MQ=255)
                                                                                aTGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAAC‑AAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCt                                      <  2:85019/149‑1 (MQ=255)
                                                                                    ccTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGt                                    >  1:7922/1‑148 (MQ=255)
                                                                                      tGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTAc                                  >  2:341510/1‑148 (MQ=255)
                                                                                        tagtagCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGATGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACaa                                <  1:291450/148‑1 (MQ=255)
                                                                                           gAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGccccc                                                          <  1:125767/119‑1 (MQ=255)
                                                                                           gAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGccccc                                                          >  2:125767/1‑119 (MQ=255)
                                                                                               gggTTCACCCGCTAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGAt                         <  3:283520/148‑1 (MQ=255)
                                                                                               gggTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGATGTACAACCACGATg                        >  3:297505/1‑149 (MQ=255)
                                                                                                 gTTCACCCGTTAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGc                      <  1:122648/149‑1 (MQ=255)
                                                                                                 gTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTa                                            <  4:157939/127‑1 (MQ=255)
                                                                                                 gTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTa                                            >  3:157939/1‑127 (MQ=255)
                                                                                                       ccgcgaagcggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCAGACACTCGCCTGCACATTTGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCaccacc                >  7:77312‑M2/16‑149 (MQ=255)
                                                                                                       ccgcgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCaccacc                >  8:26656‑M2/16‑149 (MQ=255)
                                                                                                       ccgcgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCaccacc                >  2:98178‑M2/16‑149 (MQ=255)
                                                                                                          cgaagaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCACCACCCTc             <  4:297505‑M2/137‑1 (MQ=255)
                                                                                                             agaggccggTACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCACCACCCTCAAc          <  4:351316‑M2/140‑1 (MQ=255)
                                                                                                             agaggccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCACCACCCTCAAc          <  2:319249‑M2/140‑1 (MQ=255)
                                                                                                                 gccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTaa                                                                  <  3:236509‑M2/84‑1 (MQ=255)
                                                                                                                 gccggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTaa                                                                  >  4:236509‑M2/6‑89 (MQ=255)
                                                                                                                   cggCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCACCACCCTCAACGCCGAc    <  2:7922‑M2/146‑1 (MQ=255)
                                                                                                                     gCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCACCACCCTCAACGCCGACCt  >  6:290466‑M2/2‑149 (MQ=255)
                                                                                                                     gCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCACCGCCAGCCACGCCTACGAGCTGTACAACCACGAACCCACCACCCTcaacgccccact  >  5:30509‑M2/2‑144 (MQ=255)
                                                                                                                      |                                                                                                                                                   
AATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTACGAGCTGTACAACCACGATGGCACCACCCTCAACGCCGACCT  >  NC_002947/3928302‑3928567

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