Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,026,485 G→C 14.0% T161T (ACC→ACG) ‡ dppA‑II ← dipeptide ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,026,4850GC14.0% 108.6 / 10.5 43T161T (ACC→ACG) ‡dppA‑IIdipeptide ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base C (4/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

GAACACGAACGGCCCGGTACCGATAGGCTGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATGCGGTTGAAGGTG  >  NC_002947/1026338‑1026623
                                                                                                                                                   |                                                                                                                                           
gAACACGAACGGCCCGGTACCGATAGGCTGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAggg                                                                                                                                            <  1:7026/149‑1 (MQ=255)
gAACACGAACGGCCCGGTACCGATAGGCTGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAgg                                                                                                                                             <  1:202685/148‑1 (MQ=255)
gAACACGAACGGCCCGGTACCGATAGGCTGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAgg                                                                                                                                             <  7:191579/148‑1 (MQ=255)
         cGGCCCGGTACCGATAGGCTGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAAcacca                                                                                                                                   <  4:148368/149‑1 (MQ=255)
           gCCTGGTACCTATATGATGCTGATTGATGTCGCTGGGGCGTACACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCACCCTGAAcaccgtc                                                                                                                                 <  8:124877/149‑4 (MQ=37)
           gCCCGGTACCGATAGGCTGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCACCCTGAAcaccgtc                                                                                                                                 <  4:55309/149‑4 (MQ=255)
                      aTAGGCGGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGggtc                                                                                                                      <  6:14223/149‑1 (MQ=255)
                          gctgctGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAAGGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGggtcggtc                                                                                                                  >  8:40421/1‑149 (MQ=255)
                          gctgctGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGggtcggtc                                                                                                                  >  6:74230/1‑149 (MQ=255)
                              ctgATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGggt                                                                                                                       <  8:248053/140‑1 (MQ=255)
                              ctgATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGggt                                                                                                                       >  7:248053/1‑140 (MQ=255)
                               tgattgatGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGtcttc                                                                                                              <  6:111325/147‑1 (MQ=255)
                               tgattgatGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGtcttc                                                                                                              >  5:111325/1‑147 (MQ=255)
                                   tgatGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGAGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGaca                                                                                                        >  7:105566/1‑149 (MQ=255)
                                              gggCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATg                                                                                               <  2:287970/147‑1 (MQ=255)
                                                  gTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTg                                                                                         <  5:74230/149‑1 (MQ=255)
                                                  gTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTg                                                                                         <  6:262373/149‑1 (MQ=255)
                                                    ccACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTc                                                                                       >  2:19071/1‑149 (MQ=255)
                                                         ggCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGc                                                                                  >  4:136288/1‑149 (MQ=255)
                                                                  ctgctCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGt                                                                         <  7:221548/149‑1 (MQ=255)
                                                                     ctCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTg                                                                       >  1:48955/1‑148 (MQ=255)
                                                                        ggCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGaa                                                                     <  3:108207/147‑1 (MQ=255)
                                                                         gCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATa                                                                  <  3:106200/149‑1 (MQ=255)
                                                                          cGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAg                                                                 <  7:133257/149‑1 (MQ=255)
                                                                                    ggAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCgg                                                       >  4:123429/1‑149 (MQ=255)
                                                                                      agagGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGtt                                                                        <  6:164326/130‑1 (MQ=255)
                                                                                      agagGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGtt                                                                        >  5:164326/1‑130 (MQ=255)
                                                                                        agGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGCGCGATGTTCTTGTCCAGGCGCATGCCGTTGAAATAGGGAAGCTCggtagg                                                   >  4:98578/1‑149 (MQ=255)
                                                                                            tGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCAAGGCCAATGCCGTAGAAATAGGGAAACTCggtaggggag                                               >  8:55719/1‑149 (MQ=255)
                                                                                                ggCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCgg                                                       >  1:123913/1‑137 (MQ=255)
                                                                                                ggCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCgg                                                       <  2:123913/137‑1 (MQ=255)
                                                                                                 gCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGcc                                                                                 <  4:262049/110‑1 (MQ=255)
                                                                                                 gCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGcc                                                                                 >  3:262049/1‑110 (MQ=255)
                                                                                                     aGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCAGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACgg                                      <  8:105566/149‑1 (MQ=255)
                                                                                                          aTGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAAc                                   <  3:123429/147‑1 (MQ=255)
                                                                                                                aGGTTCTGGATGAATGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGAcc                           <  7:40421/149‑1 (MQ=255)
                                                                                                                    tCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTg                       >  1:26265/1‑149 (MQ=255)
                                                                                                                     cTGGATGAAGGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTg                       <  7:195641/148‑1 (MQ=255)
                                                                                                                        gATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGa                   <  2:48955/149‑1 (MQ=255)
                                                                                                                             aCGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATg              <  3:136288/149‑1 (MQ=255)
                                                                                                                                   cATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATGCGGtt         >  5:22552/1‑148 (MQ=255)
                                                                                                                                     accgtggTGTTCACCCTGAACACCACGGTCAATGGGTCGGTCAGTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATGCGGTTg        >  1:214429/6‑148 (MQ=37)
                                                                                                                                       cgtggTGTTCACCCTGAACACCACGGTCAAAGGGTCCGTCTTCTCGACACCGGCGATCTTCTTGTCCAGGCCCATAACCTTGAAATAGGGAAACTCGGTAGGCTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATCCGGTAGAAgg    >  6:247483/4‑149 (MQ=255)
                                                                                                                                        gtggTGTTCACCCTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATGCGGTTGAAGGt   >  5:50848/3‑149 (MQ=255)
                                                                                                                                         tggTGTTCACCCTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATGCGGTTGAAGGt   >  5:203989/2‑148 (MQ=255)
                                                                                                                                         cgggcttCAGGGTGAACACCACGGTCAATGTGTCGGTCTTCTCGACACGGGCGATGTGCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATGCGGTTGAAGGTg  <  4:16651/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GAACACGAACGGCCCGGTACCGATAGGCTGCTGATTGATGTCGCTGGGGCGTCCACTGGCCAGTAGCTGCTCGGCGTATTCAGCGGAGAGGATGGAGGCGAAGCTCATGGCCAGGTTCTGGATGAACGCGGCATCGACGGTGTTCAGGGTGAACACCACGGTCAATGGGTCGGTCTTCTCGACACGGGCGATGTTCTTGTCCAGGCCCATGCCGTTGAAATAGGGAAACTCGGTAGGGTAGGCCTGACGGAACCGATGACCCTTGTCGAGCATGCGGTTGAAGGTG  >  NC_002947/1026338‑1026623

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: