Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,266,179 C→T 15.9% V89I (GTT→ATT) ‡ tufB ← protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,266,1790CT15.9% 166.2 / 19.3 63V89I (GTT→ATT) ‡tufBprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  ref base C (23/30);  new base T (4/6);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGG  >  Exported/2266111‑2266222
                                                                    |                                           
gCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACa                                            >  1:1250902/1‑70 (MQ=35)
  cGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         >  1:1671253/1‑71 (MQ=34)
                      aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGt                      >  1:1879085/1‑70 (MQ=11)
                      aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:297206/71‑1 (MQ=11)
                      aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:712045/71‑1 (MQ=11)
                      aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:268631/71‑1 (MQ=11)
                      aGGATCGCGCCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTc                     <  1:4985370/71‑1 (MQ=11)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:2950939/49‑1 (MQ=38)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:4657540/49‑1 (MQ=38)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:3536029/49‑1 (MQ=38)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAg                                         <  1:178601/49‑1 (MQ=38)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:801683/64‑1 (MQ=35)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:2677146/64‑1 (MQ=35)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:2438177/64‑1 (MQ=35)
                        gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                          <  1:3995979/64‑1 (MQ=35)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:4300656/1‑49 (MQ=38)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:3848834/1‑49 (MQ=38)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:535194/1‑49 (MQ=38)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:2926817/1‑49 (MQ=38)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:3052658/1‑49 (MQ=38)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:2817966/1‑49 (MQ=38)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:1080182/1‑49 (MQ=38)
                          tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                       >  1:3045157/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)
                                                                    |                                           
GCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGG  >  Exported/2266111‑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: