Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,266,174 T→C 15.1% K90K (AAA→AAG tufB ← protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,266,1740TC15.1% 167.6 / 19.5 66K90K (AAA→AAGtufBprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  ref base T (23/33);  new base C (4/6);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

GGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGG  >  minE/2266109‑2266222
                                                                 |                                                
gggCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTaa                                              <  1:1389953/70‑1 (MQ=35)
gggCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTaa                                              <  1:2551182/70‑1 (MQ=35)
gggCCGTCAGTCGCAGCAACTAACAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTaa                                              <  1:1137745/70‑1 (MQ=255)
  gCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACa                                            >  1:1250902/1‑70 (MQ=35)
    cGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         >  1:1671253/1‑71 (MQ=34)
                        aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGt                      >  1:1879085/1‑70 (MQ=11)
                        aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:4985370/71‑1 (MQ=11)
                        aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:268631/71‑1 (MQ=11)
                        aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:712045/71‑1 (MQ=11)
                        aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:297206/71‑1 (MQ=11)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:4657540/49‑1 (MQ=38)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:2950939/49‑1 (MQ=38)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:3536029/49‑1 (MQ=38)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:178601/49‑1 (MQ=38)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:2677146/64‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:2438177/64‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:3995979/64‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:801683/64‑1 (MQ=35)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:4300656/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:535194/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:2926817/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:2817966/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:3045157/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:1080182/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:3052658/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:3848834/1‑49 (MQ=38)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcc                               >  1:3194415/1‑57 (MQ=35)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                  >  1:1514593/1‑70 (MQ=35)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                 >  1:4584516/1‑71 (MQ=34)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                 >  1:3369030/1‑71 (MQ=34)
                            tCGCGCCGTCCATCTGCGCAGCACAGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                 >  1:2615751/1‑71 (MQ=21)
                                 ccGTCCATCTGCGCAGCACCGGTGATCATGTTTTCAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGc            >  1:2021762/1‑71 (MQ=21)
                                  cggccATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:1234672/68‑1 (MQ=21)
                                  cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:5083395/71‑1 (MQ=34)
                                  cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:3536838/71‑1 (MQ=34)
                                  cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:1024165/71‑1 (MQ=34)
                                  cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:2576480/71‑1 (MQ=34)
                                   gTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTTCCCCGGGCAGTCTACGTGTGCGt          <  1:116911/71‑1 (MQ=21)
                                   gTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGt          <  1:1005779/71‑1 (MQ=34)
                                   gTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGt          <  1:758690/71‑1 (MQ=34)
                                    tCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTa         <  1:4593377/71‑1 (MQ=34)
                                     ccATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGCGCGTa         >  1:2374542/1‑70 (MQ=2)
                                     ccATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGCGCGTa         >  1:2323996/1‑70 (MQ=2)
                                     ccATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGCGCGTAg        >  1:4104371/1‑71 (MQ=2)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCg                                   >  1:3382242/1‑43 (MQ=38)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCg                                   >  1:5064879/1‑43 (MQ=38)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCg                                   >  1:570476/1‑43 (MQ=38)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAg        <  1:2274659/70‑1 (MQ=35)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAg        <  1:1112507/70‑1 (MQ=35)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       >  1:4327106/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       >  1:2603834/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGAGCGTAGt       <  1:4400135/71‑1 (MQ=2)
                                      cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGAGCGTAGt       <  1:1730082/71‑1 (MQ=2)
                                        tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:88499/69‑1 (MQ=35)
                                        tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:1075632/69‑1 (MQ=35)
                                        tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:770768/69‑1 (MQ=35)
                                        tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:1465794/69‑1 (MQ=35)
                                        tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:4568808/69‑1 (MQ=35)
                                        tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:2374937/69‑1 (MQ=35)
                                        tCTGCGCAGCACCGGTGATCATGTGTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:738858/69‑1 (MQ=25)
                                         cTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:3748934/68‑1 (MQ=35)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTGTACGTGTGCGTAGTGACg   >  1:2382090/1‑71 (MQ=21)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACg   >  1:2796724/1‑71 (MQ=34)
                                            cgcAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACgg  <  1:1102974/70‑1 (MQ=31)
                                            cgcAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACgg  <  1:1077606/70‑1 (MQ=31)
                                            cgcAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACgg  <  1:3849002/70‑1 (MQ=31)
                                                                 |                                                
GGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGG  >  minE/2266109‑2266222

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