Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,173 A→T 10.4% noncoding (168/184 nt) glmY ← sRNA activator of glmS mRNA, glmZ processing antagonist

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,691,1730AT10.4% 190.9 / 7.8 59noncoding (168/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base T (4/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATA  >  NC_000913/2691029‑2691321
                                                                                                                                                 |                                                                                                                                                    
gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCg                                                                                                                                                  >  4:557651/1‑149 (MQ=255)
    aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCtt                                                                                                                                              <  2:337275/149‑1 (MQ=255)
    aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCAGGCtt                                                                                                                                              >  4:557942/1‑149 (MQ=255)
         gCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTa                                                                                                                                         >  3:349659/1‑149 (MQ=255)
          cAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCTTTGttgt                                                                                                                                        <  4:411570/149‑3 (MQ=255)
          cAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                                        <  4:207968/149‑1 (MQ=255)
                  ccAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAg                                                                                                                                <  1:215638/149‑1 (MQ=255)
                       aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                           <  3:557651/149‑1 (MQ=255)
                       aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                           <  3:126264/149‑1 (MQ=255)
                              gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAcc                                                                                                                                                    >  3:453869/1‑117 (MQ=255)
                              gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAcc                                                                                                                                                    <  4:453869/117‑1 (MQ=255)
                                tGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTc                                                                                                                  >  3:71623/1‑149 (MQ=255)
                                    gTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACg                                                                                                              >  4:430921/1‑149 (MQ=255)
                                     tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGt                                                                                                             <  3:302901/149‑1 (MQ=255)
                                           gCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGcc                                                                                                       <  4:264407/149‑1 (MQ=255)
                                             ggggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGGCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCg                                                                                                     <  3:557942/149‑1 (MQ=255)
                                                   aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                               >  1:94514/1‑149 (MQ=255)
                                                    aaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTg                                                                                              >  3:292066/1‑149 (MQ=255)
                                                        ggCCAGCGTTTCAAGGTGTTACTCTCGACAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGatat                                                                                          >  4:311095/1‑149 (MQ=255)
                                                         gCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAccc                                                                                                                                                   >  3:423707/1‑91 (MQ=255)
                                                         gCCAGCGTTTCAAGGTGTTACTCTCGTAAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAccc                                                                                                                                                   <  4:423707/91‑1 (MQ=255)
                                                           cAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGc                                                                                       >  3:135662/1‑149 (MQ=255)
                                                             gCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata                                                                                           >  1:362776/1‑143 (MQ=255)
                                                             gCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata                                                                                           <  2:362776/143‑1 (MQ=255)
                                                                   cAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTg                                                                                                                                            >  2:405039/1‑88 (MQ=38)
                                                                   cAAGGTGTTACTCTCGTCAGACGCGAATAGCCTAATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTg                                                                                                                                            <  1:405039/88‑1 (MQ=38)
                                                                      ggTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt                                                                                                                                              <  3:325236/83‑1 (MQ=39)
                                                                      ggTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt                                                                                                                                              >  4:325236/1‑83 (MQ=38)
                                                                           tACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa                                                                       <  3:311095/149‑1 (MQ=255)
                                                                           tACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa                                                                       <  1:177478/149‑1 (MQ=255)
                                                                              tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                       <  3:294539/82‑1 (MQ=255)
                                                                              tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                       >  4:294539/1‑82 (MQ=255)
                                                                                  gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                                >  4:499278/1‑149 (MQ=255)
                                                                                              tAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGAt                                                                                                   <  3:192892/102‑1 (MQ=255)
                                                                                              tAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGAt                                                                                                   >  4:192892/1‑102 (MQ=255)
                                                                                                   tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                               >  2:151615/1‑149 (MQ=255)
                                                                                                    gATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCgg                                                                                                                                                 >  2:470521/1‑50 (MQ=25)
                                                                                                     aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                             <  1:391441/149‑1 (MQ=255)
                                                                                                          aaCCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt                                                                     <  2:543153/120‑1 (MQ=255)
                                                                                                          aaCCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt                                                                     >  1:543153/1‑120 (MQ=255)
                                                                                                                  ggAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAttt                                <  3:499278/149‑1 (MQ=255)
                                                                                                                      gTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATg                            <  3:430921/149‑1 (MQ=255)
                                                                                                                           gATACAACCAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                               <  3:456919/77‑1 (MQ=255)
                                                                                                                           gATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                               >  4:456919/1‑77 (MQ=255)
                                                                                                                            aTACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCAc                      >  4:61345/1‑149 (MQ=255)
                                                                                                                                aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata                                                                                         >  3:485480/1‑78 (MQ=255)
                                                                                                                                aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata                                                                                         <  4:485480/78‑1 (MQ=255)
                                                                                                                                aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                      <  4:487871/129‑1 (MQ=255)
                                                                                                                                aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                      >  3:487871/1‑129 (MQ=255)
                                                                                                                                aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                  <  4:71623/149‑1 (MQ=255)
                                                                                                                                 acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                            >  3:144327/1‑106 (MQ=255)
                                                                                                                                 acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                            <  4:144327/106‑1 (MQ=255)
                                                                                                                                 acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTccc                 >  4:44297/1‑149 (MQ=255)
                                                                                                                                    aaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa              >  3:105739/1‑149 (MQ=255)
                                                                                                                                     aGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg             >  2:566578/1‑149 (MQ=255)
                                                                                                                                       ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAggg           >  3:203305/1‑149 (MQ=255)
                                                                                                                                       ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAggg           >  1:346353/1‑149 (MQ=255)
                                                                                                                                         ggg‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGc         >  2:76443/1‑149 (MQ=255)
                                                                                                                                         ggg‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGc         <  4:135662/149‑1 (MQ=255)
                                                                                                                                                 aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa                                                                               >  3:621524/1‑72 (MQ=255)
                                                                                                                                                 aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa                                                                               <  4:621524/72‑1 (MQ=255)
                                                                                                                                                 aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATa  <  2:94514/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATA  >  NC_000913/2691029‑2691321

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