Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 660,596 T→C 9.8% Y548H (TAC→CAC)  speA → biosynthetic arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947660,5960TC 9.8% 134.8 / 5.6 51Y548H (TAC→CAC) speAbiosynthetic arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (20/26);  new base C (3/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGTGCACCTG  >  NC_002947/660459‑660739
                                                                                                                                         |                                                                                                                                               
tATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATc                                                                                                                                      <  7:207709/149‑1 (MQ=255)
tATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGGGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATc                                                                                                                                      >  4:65693/1‑149 (MQ=255)
     cAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTgg                                                                                                                                  <  8:118891/148‑1 (MQ=255)
     cAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTgg                                                                                                                                  <  4:310417/148‑1 (MQ=255)
         aCCTGAATTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTACTGGTGGGTGCCTAGCAGGAAATCCTGGGTGa                                                                                                                              >  4:36105/1‑148 (MQ=255)
          ccTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACa                                                                                                                            <  2:216882/149‑1 (MQ=255)
           cTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACAt                                                                                                                           >  1:1703/1‑149 (MQ=255)
                ttGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGcac                                                                                                                       >  5:200647/1‑148 (MQ=255)
                   cGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAAc                                                                                                                    <  6:200647/148‑1 (MQ=255)
                     aTTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGATTTCCTggtag                                                                                                                               <  1:54845/135‑4 (MQ=255)
                     aTTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGATTTCCTggtag                                                                                                                               >  2:54845/1‑132 (MQ=255)
                      ttCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGt                                                                                                                >  2:206364/1‑149 (MQ=255)
                         cGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCg                                                                                                             <  3:65693/149‑1 (MQ=255)
                               cAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAGGaa                                                                                                                                         >  6:237305/1‑115 (MQ=255)
                               cAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAGGaa                                                                                                                                         <  5:237305/115‑1 (MQ=255)
                                aaGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAGGaac                                                                                                                                        >  4:218905/1‑114 (MQ=255)
                                aaGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAGGaac                                                                                                                                        <  3:218905/115‑2 (MQ=255)
                                  gATCAAGCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACAcc                                                                                                     >  2:133351/1‑148 (MQ=255)
                                       aCCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCg                                                                                               <  2:1703/149‑1 (MQ=255)
                                        ccAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCgg                                                                                              <  8:116437/149‑1 (MQ=255)
                                         cAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGt                                                                                             <  5:220625/149‑1 (MQ=255)
                                           gTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGtt                                                                                                               >  7:229112/1‑129 (MQ=255)
                                           gTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGtt                                                                                                               <  8:229112/129‑1 (MQ=255)
                                           gTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTACCTGGTGGGTGCCTACCAGGAAATCCAGGGTGACATGCACAACCTGTACGGTGACACCGATTCGGtc                                                                                            >  1:203555/1‑147 (MQ=255)
                                              cGTTGACGAGCAGTGCATCGAAACCAGCATGCCGGTGCATGCGTTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGTCATGCACAACCTGTTCGGTGACACCGATTCTGTGAAc                                                                                         <  3:36105/148‑1 (MQ=255)
                                                 tGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTg                                                                                            <  1:295480/142‑1 (MQ=255)
                                                 tGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTg                                                                                            >  2:295480/1‑142 (MQ=255)
                                                                          aTGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTa                                                            >  6:102012/1‑149 (MQ=255)
                                                                              cGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCAc                                                        >  3:298096/1‑149 (MQ=255)
                                                                                           gAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGaa                                            <  3:294178/148‑1 (MQ=255)
                                                                                           gAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGaa                                            >  4:294178/1‑148 (MQ=255)
                                                                                              ggaggGTGAGGCTTACCTGCTGGGCGTGTTCCTGGGGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAAccc                                        <  8:307731/149‑1 (MQ=255)
                                                                                                aggGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCAc                                      <  6:206849/149‑1 (MQ=255)
                                                                                                    tGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATc                                               >  3:96702/1‑136 (MQ=255)
                                                                                                    tGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATc                                               <  4:96702/136‑1 (MQ=255)
                                                                                                        gATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCAt                              <  1:206364/149‑1 (MQ=255)
                                                                                                               tgctgGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGaa                          >  3:9590/1‑146 (MQ=255)
                                                                                                               tgctgGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGaa                          <  4:9590/146‑1 (MQ=255)
                                                                                                                gctgGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCg                                                    <  6:121432/119‑1 (MQ=255)
                                                                                                                    caggatttCCTGGTAGGCACCCACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACAtgct                  >  1:128663/7‑149 (MQ=255)
                                                                                                                      cgggtTCCTGGTGGGTGCCTTCCAGGAAATCCTGGGGGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCGACGACACCATCGAAGACTtgctgc                <  2:203555/146‑1 (MQ=255)
                                                                                                                      cGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGACACCCACGACACCATCGAAGACAtgctgc                >  8:164649/1‑149 (MQ=255)
                                                                                                                      cGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAAccc                                        <  7:163101/125‑1 (MQ=255)
                                                                                                                      cGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAAccc                                        >  8:163101/1‑125 (MQ=255)
                                                                                                                      cGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACAtgctgc                <  4:109137/149‑1 (MQ=255)
                                                                                                                          ttCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGa                                    >  8:40834/1‑125 (MQ=255)
                                                                                                                          ttCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGa                                    <  7:40834/125‑1 (MQ=255)
                                                                                                                            ccTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTAc           <  1:133351/148‑1 (MQ=255)
                                                                                                                             cTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGt         >  6:140909/1‑149 (MQ=255)
                                                                                                                                gtgggtgCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGTGCa      >  4:110753/1‑149 (MQ=255)
                                                                                                                                    gtgCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGTGCACCTg  <  3:110753/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                               
TATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGTGCACCTG  >  NC_002947/660459‑660739

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