Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 660,941 A→T 9.6% intergenic (+73/‑94) speA → / → PP_0568 biosynthetic arginine decarboxylase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947660,9410AT 9.6% 138.9 / 6.1 52intergenic (+73/‑94)speA/PP_0568biosynthetic arginine decarboxylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base T (3/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGG  >  NC_002947/660795‑661087
                                                                                                                                                  |                                                                                                                                                   
cTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAg                                                                                                                                                   <  6:7122/149‑1 (MQ=255)
     cgTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTc                                                                                                                                              <  7:217209/149‑1 (MQ=255)
         aGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCAc                                                                                                                                            >  1:92369/1‑147 (MQ=255)
          gCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCatatc                                                                                                                                         <  8:270228/149‑5 (MQ=255)
            ggACCCACTATCTGGATGCATTGGGCCTGGGCTTGACCCGGTCTTCTTCCCTCTCGTCCTAATTGTCGGGGCCGCTCCTGCAAAAGCCTGTCCTGCCTTTTGTGGGACCGGGCTTGCGTCGCTATATGAGGGCGAACCCCTCACCAGCg                                                                                                                                       <  6:125071/149‑1 (MQ=255)
               cccAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTca                                                                                                                                    >  4:62028/1‑149 (MQ=255)
                  aGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACaa                                                                                                                                 <  3:146301/149‑1 (MQ=255)
                      tCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCg                                                                                                                              <  2:119301/148‑1 (MQ=255)
                            tGCCTTGCGCCTTGGGTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAATGCCTGTGCTGCCTTTTGTGGTAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAata                                                                                                                       <  5:74368/149‑1 (MQ=255)
                                  gcgcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCgg                                                                                                                  <  2:92369/148‑1 (MQ=255)
                                  gcgcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCg                                                                                                                   <  6:125337/147‑1 (MQ=255)
                                     ccTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCt                                                                                                              <  1:254788/149‑1 (MQ=255)
                                      cTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                                                             >  4:133442/1‑149 (MQ=255)
                                            ttGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCAt                                                                                                       >  1:67846/1‑149 (MQ=255)
                                              gACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCAt                                                                                                       <  8:7136/147‑1 (MQ=255)
                                               aCCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGc                                                                                                    <  3:133442/149‑1 (MQ=255)
                                                               tGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCt                                                                                    >  2:27549/1‑149 (MQ=255)
                                                                    tcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTg                                                                                 >  7:261144/1‑147 (MQ=255)
                                                                    tcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTg                                                                                 <  8:261144/147‑1 (MQ=255)
                                                                    tcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGaa                                                                               >  6:312793/1‑149 (MQ=255)
                                                                          ttGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAt                                                                         >  6:146006/1‑149 (MQ=255)
                                                                                  gtcgatCCTGCTATTGCCTGTGCTGGCTGTTGTGGGAGCGGCCTTGCGTAGCGATATGTGGGCGAAGGCCTCACCAGCGTCACAAGCGGTAAATATCGGGCCCTGGGCCATCGGACCTTTGCCGGTACCTTTGAATGCATT‑GTTCaaa                                                                 <  6:318890/143‑1 (MQ=255)
                                                                                  gCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAag                                                                <  8:125967/149‑1 (MQ=255)
                                                                                  gCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAag                                                                <  5:146006/149‑1 (MQ=255)
                                                                                   ccGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGAAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAagg                                                               <  8:119182/149‑1 (MQ=255)
                                                                                     gCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCaaa                                                                 >  5:318890/1‑145 (MQ=255)
                                                                                         cTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGAc                                                          >  5:49926/1‑148 (MQ=255)
                                                                                             aaaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGatat                                                     >  2:86671/1‑149 (MQ=255)
                                                                                                  ccTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGata                                                      <  6:97564/143‑1 (MQ=255)
                                                                                                  ccTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGata                                                      >  5:97564/1‑143 (MQ=255)
                                                                                                        ctgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAg                                          <  2:67846/149‑1 (MQ=255)
                                                                                                         tgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAgg                                         >  4:48441/1‑149 (MQ=255)
                                                                                                          gcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCaa                                                                  <  3:210080/123‑1 (MQ=255)
                                                                                                          gcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGAAAT‑GTTCaa                                                                  >  4:210080/1‑123 (MQ=255)
                                                                                                            cTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTc                                      <  5:312793/149‑1 (MQ=255)
                                                                                                             ttttGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTcc                                     >  6:119407/1‑149 (MQ=255)
                                                                                                                 gtgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATc                                 >  4:301169/1‑149 (MQ=255)
                                                                                                                   gggAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGt                               >  2:252734/1‑149 (MQ=255)
                                                                                                                    ggAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGtt                              >  4:77898/1‑149 (MQ=255)
                                                                                                                    ggAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGt                               <  1:276900/148‑1 (MQ=255)
                                                                                                                          gCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCATGGTCCATCGCACAAATGCCGCTACATTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCAAGGATCCCATCGTTTACGTg                        >  2:234096/1‑149 (MQ=255)
                                                                                                                                    gcgATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATATTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa                <  1:86671/148‑1 (MQ=255)
                                                                                                                                    gcgATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAgt              <  1:27549/149‑1 (MQ=255)
                                                                                                                                       tggtGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTc           >  3:268047/4‑149 (MQ=255)
                                                                                                                                       tggtGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAACTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTc           >  7:194005/4‑149 (MQ=255)
                                                                                                                                        ggtGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCt          >  5:19066/3‑149 (MQ=255)
                                                                                                                                         gtGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCt          <  6:19066/147‑1 (MQ=255)
                                                                                                                                             gggCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGt                                              >  3:69325/1‑108 (MQ=255)
                                                                                                                                             gggCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGt                                              <  4:69325/108‑1 (MQ=255)
                                                                                                                                             gggCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTc           >  6:293878/1‑143 (MQ=255)
                                                                                                                                             gggCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTc           <  5:293878/143‑1 (MQ=255)
                                                                                                                                                cGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCgg  <  5:119407/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAT‑GTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGG  >  NC_002947/660795‑661087

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