Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,173 A→T 14.6% 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,1730AT14.6% 118.2 / 10.9 41noncoding (168/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base A (15/20);  new base T (3/3);  total (18/23)
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.00e+00

GCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACAT  >  NC_000913/2691026‑2691320
                                                                                                                                                   |                                                                                                                                                   
gcagcaGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTAc                                                                                                                                                    >  1:456851/1‑149 (MQ=255)
      gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCt                                                                                                                                              <  4:79332/149‑1 (MQ=255)
       aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCtt                                                                                                                                             <  4:115929/149‑1 (MQ=255)
        aTCAGCAAAGATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAGAAGCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGGCGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACTACAAAGCCGGGTAATTCCCGGCttt                                                                                                                                            <  3:203230/149‑1 (MQ=255)
        aTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCttt                                                                                                                                            <  3:454169/149‑1 (MQ=255)
              aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATg                                                                                                                                      <  2:305232/149‑1 (MQ=255)
                 aaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCCTTGTTATGGaa                                                                                                                                   <  3:307850/149‑1 (MQ=255)
                          aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                          <  2:139937/149‑1 (MQ=255)
                          aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                          <  3:270278/149‑1 (MQ=255)
                             aCCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGcc                                                                                                                       <  4:196550/149‑1 (MQ=255)
                              ccAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                      >  1:62008/1‑149 (MQ=255)
                               cAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTa                                                                                                                     >  1:6401/1‑149 (MQ=255)
                                   tGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTc                                                                                                                 <  4:273963/149‑1 (MQ=255)
                                      cGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGAc                                                                                                              <  2:456851/149‑1 (MQ=255)
                                           tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                         >  4:78735/1‑149 (MQ=255)
                                           tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                         <  3:78735/149‑1 (MQ=255)
                                                      aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                              <  1:37022/149‑1 (MQ=255)
                                                          gggCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata                                                                                          >  3:395061/1‑149 (MQ=255)
                                                                   tttCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGc                                                                                 <  2:62008/149‑1 (MQ=255)
                                                                          gtgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATc                                                                          <  4:395061/149‑1 (MQ=255)
                                                                                         gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                           <  3:344130/149‑1 (MQ=255)
                                                                                         gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                           <  1:52325/149‑1 (MQ=255)
                                                                                            gcgAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                      >  3:442727/1‑87 (MQ=255)
                                                                                            gcgAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                      <  4:442727/87‑1 (MQ=255)
                                                                                                    ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                                <  2:148293/149‑1 (MQ=255)
                                                                                                      tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                              >  3:296272/1‑149 (MQ=255)
                                                                                                      tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                              >  3:186139/1‑149 (MQ=255)
                                                                                                                     ggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAttt                               <  1:437164/149‑1 (MQ=255)
                                                                                                                             aGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGc                       <  4:186139/149‑1 (MQ=255)
                                                                                                                              gATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCa                      >  2:73561/1‑149 (MQ=255)
                                                                                                                               aTACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCAc                     >  1:456588/1‑149 (MQ=255)
                                                                                                                                 acaacaAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGt                   >  1:98362/1‑149 (MQ=255)
                                                                                                                                      aaaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                 >  1:422686/1‑146 (MQ=255)
                                                                                                                                      aaaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                 <  2:422686/146‑1 (MQ=255)
                                                                                                                                        aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg            >  3:102426/1‑149 (MQ=255)
                                                                                                                                        aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg            >  4:454277/1‑149 (MQ=255)
                                                                                                                                        aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg            >  3:125258/1‑149 (MQ=255)
                                                                                                                                        aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg            >  4:375266/1‑149 (MQ=255)
                                                                                                                                               aaTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGa     <  2:456588/149‑1 (MQ=255)
                                                                                                                                                 ttACCCGGCTTTGTTATGGAATAAGGAGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATAACAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACa   >  4:408154/1‑149 (MQ=255)
                                                                                                                                                  tACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACAt  >  2:17244/1‑149 (MQ=255)
                                                                                                                                                   aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                            <  4:319016/106‑1 (MQ=255)
                                                                                                                                                   aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                            >  3:319016/1‑106 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACAT  >  NC_000913/2691026‑2691320

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