Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,210,827 T→A 19.7% intergenic (+17/+20) ytnA → / ← PP_1060 amino acid permease YtnA/glutamate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,210,8270TA19.7% 56.0 / 9.9 31intergenic (+17/+20)ytnA/PP_1060amino acid permease YtnA/glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (6/17);  new base A (4/2);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GGCGTACTTCGAGGATACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAACACATGCATCTGC  >  NC_002947/1210681‑1210973
                                                                                                                                                  |                                                                                                                                                  
ggCGTACTTCGAGGATACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCgtt                                                                                                                                                  >  2:90805/1‑149 (MQ=255)
    tACTTCGAGGATACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCgg                                                                                                                                              <  6:146631/149‑1 (MQ=255)
      cTTCGAGGATACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGcc                                                                                                                                            <  1:309728/149‑1 (MQ=255)
       ttCGAGGATACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTAGTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGGCGCGGTTCGGccc                                                                                                                                           >  2:296769/1‑149 (MQ=255)
          gAGGATACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCAAAGCGGCCCCac                                                                                                                                        <  5:21918/149‑2 (MQ=255)
              aTACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGtt                                                                                                                                    <  1:296769/149‑1 (MQ=255)
                   cgTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGc                                                                                                                                                >  4:185638/1‑132 (MQ=255)
                   cgTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGc                                                                                                                                                <  3:185638/132‑1 (MQ=255)
                         gCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGgccg                                                                                                                         <  3:196999/149‑1 (MQ=255)
                           ccTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCAAAGCGGcccc                                                                                                                                          <  7:60069/130‑1 (MQ=255)
                           ccTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCAAAGCGGcccc                                                                                                                                          >  8:60069/1‑130 (MQ=255)
                               tACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTAGAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACt                                                                                                                    <  4:46904/148‑1 (MQ=255)
                                      gCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCtgtg                                                                                                            <  5:350677/149‑1 (MQ=255)
                                                 cTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGaaa                                                                                                  <  7:91758/148‑1 (MQ=255)
                                                                                               aaGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGtt                                                   <  7:130777/149‑1 (MQ=255)
                                                                                                 gACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCACTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCa                                                 <  2:198010/149‑1 (MQ=255)
                                                                                                   cGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCAcc                                               <  5:161321/149‑1 (MQ=255)
                                                                                                                aCGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGTACTGGCCGTTCACCTCGCCGGTcagca                                  <  1:88432/149‑1 (MQ=255)
                                                                                                                      aCAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGGTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGTATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCa                             <  3:300679/148‑1 (MQ=255)
                                                                                                                             cctgacctgaCGTGGGGCCGCGTTTCGGCCCCACCGTTCTATAAGGCCCCCACATAGCTGTGCGGCGCGAAACAGTCCGCCCCGGCCCGTTGCCACATCCCCGAATAGAACAGGCCGATCAACAAGCCGGTCAGCAGCACACCAGGca                      >  4:213694/1‑147 (MQ=255)
                                                                                                                                 acctgacGTGGGGCCGCTTTGCGTCCACAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCTCCCGGGCCAGTTGCCATATCCTCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCaa                  <  3:266144/148‑1 (MQ=255)
                                                                                                                                   ctgacGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGaa               <  8:206924/149‑1 (MQ=255)
                                                                                                                                    aacgctGGGGCCGCAAAGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTc                              >  6:171250/6‑133 (MQ=255)
                                                                                                                                       gctGGGGCCGCAAAGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAAcaca           >  4:183029/3‑149 (MQ=255)
                                                                                                                                        gTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCtt                     >  4:201207/1‑138 (MQ=255)
                                                                                                                                        gTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCtt                     <  3:201207/138‑1 (MQ=255)
                                                                                                                                        gTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAACACAt          <  3:213694/149‑1 (MQ=255)
                                                                                                                                        ctGGGGCCGCAAAGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAACACAt          >  3:124439/2‑149 (MQ=255)
                                                                                                                                             gCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAACACATGCATc     >  7:85219/1‑149 (MQ=255)
                                                                                                                                             gCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAACACATGCATc     >  8:61026/1‑149 (MQ=255)
                                                                                                                                                gCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAACACATGCATCt    >  5:119127/1‑147 (MQ=255)
                                                                                                                                                gCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACATCGCCGGTCAGCAGCACACCTGGCTTCAAGAACACATGCATCTGc  >  4:86077/1‑149 (MQ=255)
                                                                                                                                                  tttGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCtt                     >  4:324832/1‑128 (MQ=255)
                                                                                                                                                  tttGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCtt                     <  3:324832/128‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGCGTACTTCGAGGATACGCGTGTGGCCCTGTACATCGGCCCGGCGTTCCTGGTGTTGCTGACAGTGCTGTATTACACGTTCCGCCTGGCGCCGAAAGACGCACAAGGTGTGACGAGTACAGCCTCCTGACCTGACGTGGGGCCGCTTTGCGGCCCCAGCGTTCTCTCAGGCCGCCACTTCGCTGTGCGGCTCGAAACTGTCCGCCCGGGCCAGTTGCCACATCCGCGAATAGAACTGGCCGTTCACCTCGCCGGTCAGCAGCTCACCTGGCTTCAAGAACACATGCATCTGC  >  NC_002947/1210681‑1210973

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: