Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,464 A→T 46.4% V126V (GTA→GTT)  tufA → protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,169,4640A→T46.4% 99.5 0.0 ‑0.064V126V (GTA→GTT) tufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (16/13):  ref base (19/15):  total (35/29)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.23e-01

TCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTC‑CTGAACAAATGCGACATGGTTGATGACGAAGAGCTGCTGGAACTGGTTGAAATGGAAGTTCGTGAACTTCT  >  NC_007779/4169364‑4169559
                                                                                                    |                                                                                                
tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                                                                                  >  1:10526026/1‑101 (MQ=32)
tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                                                                                  >  1:11799994/1‑101 (MQ=33)
tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                                                                                  >  1:12171896/1‑101 (MQ=33)
tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGCAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCCGGTa                                                                                                  >  1:7374746/1‑101 (MQ=12)
 cACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAg                                                                                                 >  1:8658149/1‑101 (MQ=32)
  aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                                                                                >  1:3829040/1‑101 (MQ=2)
  aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                                                                                <  1:3939415/101‑1 (MQ=2)
  aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                                                                                >  1:12897506/1‑101 (MQ=2)
  aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                                                                                <  1:1285988/101‑1 (MQ=2)
  aCCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCTGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                                                                                >  1:7421410/1‑101 (MQ=17)
  aCCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                                                                                <  1:8386249/101‑1 (MQ=32)
  aCCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                                                                                <  1:6471301/101‑1 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                                                                               >  1:4585808/1‑101 (MQ=2)
   ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                                                                               <  1:15157608/101‑1 (MQ=2)
   ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                                                                               <  1:14903809/101‑1 (MQ=2)
   ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                                                                               >  1:7780257/1‑101 (MQ=2)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               <  1:6885824/101‑1 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               >  1:3467075/1‑101 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               >  1:3297787/1‑101 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               >  1:12316217/1‑101 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               >  1:8170322/1‑101 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               <  1:5543245/101‑1 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               <  1:14044363/101‑1 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               <  1:13513714/101‑1 (MQ=32)
   ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                                                                               >  1:667615/1‑101 (MQ=32)
    cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                                                                                              >  1:899465/1‑101 (MQ=2)
    cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                                                                                              >  1:11472295/1‑101 (MQ=2)
     ggTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGt                                                                                             >  1:8086794/1‑101 (MQ=2)
     ggTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGt                                                                                             <  1:12543802/101‑1 (MQ=2)
     ggTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGt                                                                                             <  1:11185072/101‑1 (MQ=2)
     ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                                                                                             >  1:628947/1‑101 (MQ=32)
     ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                                                                                             >  1:1670189/1‑101 (MQ=32)
     ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                                                                                             >  1:13246194/1‑101 (MQ=32)
     ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                                                                                             >  1:2551651/1‑101 (MQ=32)
     ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                                                                                             <  1:7613203/101‑1 (MQ=32)
      gTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGtt                                                                                            <  1:15141921/101‑1 (MQ=2)
      gTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGtt                                                                                            >  1:11348692/1‑101 (MQ=2)
      gTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGtt                                                                                            >  1:281759/1‑101 (MQ=32)
       tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                                                                                           >  1:9212435/1‑101 (MQ=2)
       tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                                                                                           >  1:8761536/1‑101 (MQ=2)
       tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                                                                                           >  1:12252783/1‑101 (MQ=2)
       tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                                                                                           >  1:6892589/1‑101 (MQ=32)
       tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                                                                                           <  1:5431863/101‑1 (MQ=32)
       tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                                                                                           >  1:12794790/1‑101 (MQ=32)
       tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                                                                                           <  1:14624617/101‑1 (MQ=32)
        gctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTcc                                                                                          >  1:10814200/1‑101 (MQ=2)
         ctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCg                                                                                         >  1:14048703/1‑101 (MQ=2)
         ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                                                                                         <  1:2307968/101‑1 (MQ=32)
         ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCCGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                                                                                         >  1:4886678/1‑101 (MQ=17)
          tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                                                                                        >  1:1294060/1‑101 (MQ=32)
          tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                                                                                        >  1:6206477/1‑101 (MQ=32)
          tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                                                                                        <  1:5329266/101‑1 (MQ=32)
          tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                                                                                        <  1:8430094/101‑1 (MQ=32)
           gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                                                                                       <  1:12993368/101‑1 (MQ=2)
           gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                                                                                       >  1:15047943/1‑101 (MQ=2)
           gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                                                                                       <  1:15197629/101‑1 (MQ=2)
           gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                                                                                       <  1:3164264/101‑1 (MQ=2)
           gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                                                                                       <  1:7959954/101‑1 (MQ=32)
           gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                                                                                       >  1:92986/1‑101 (MQ=32)
           gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                                                                                       >  1:1140915/1‑101 (MQ=32)
            ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                                                                                      <  1:672885/101‑1 (MQ=31)
            ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                                                                                      >  1:6671999/1‑101 (MQ=31)
             tCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTAca                                                                                     >  1:12836509/1‑101 (MQ=2)
             tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTCGGCGTTCCGTTca                                                                                     <  1:10342330/101‑1 (MQ=2)
             tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAca                                                                                     <  1:3416541/101‑1 (MQ=11)
                                 tGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGACAGGTTGGCGTTCCGTACATCCTCGTTGCA‑CTGAACaag                                                                <  1:4336976/101‑2 (MQ=0)
                                                                              cttccgatctgGCTCGCCAGGTTGGCGTTCCGTACATCCTCGTTACA‑CTGAACAAGTGCGACATGGTTAACGATGAAGAGATCATCGAACTCGTCGAGATg                   <  1:209922/93‑1 (MQ=0)
                                                                              acgttctcctgGCTCGCCAGGTTGGCGTTCCGTACATCCTCGTTGCA‑CTGAACAAGTGCGACATGGTTGACGATGAAGAGATCATCGAACTCGTCGAGATg                   <  1:3594107/93‑1 (MQ=0)
                                                                                 ttctcctgGCTCGCCAGGTTGGCGTTCCGTACATCCTCGTTGCA‑CTGAACAAGTGCGACATGGTTGACGATGAAGAGATCATCGAACTCGTCGAGATGGag                <  1:7131727/96‑2 (MQ=0)
                                                                                            cgccAGGTTGGCGTTCCGTACATCCTCGTTGCA‑CTGAACAAGTGCGACATGGTTGACGATGAAGAGATCATCGAACTCGTCGAGATGGAGATTCGCGAGCt     >  1:5390705/4‑101 (MQ=0)
                                                                                            cgccAGGTTGGCGTTCCGTACATCCTCGT‑TGCACTGAACAAGTGCGACATGGTTGACGATGAAGAGATCATCGAACTCGTCGAGATGGAGATTCGCGAGCt     >  1:7159497/4‑101 (MQ=0)
                                                                                               cAGGTTGGCGTTCCGTACATCCTCGT‑TGCACTGAACAAGTGCGACATGGTTGACGATGAAGAGATCATCGAACTCGTCGAGATGGAGATTCGCGAGCTGCt  >  1:6462960/1‑101 (MQ=0)
                                                                                                    |                                                                                                
TCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTC‑CTGAACAAATGCGACATGGTTGATGACGAAGAGCTGCTGGAACTGGTTGAAATGGAAGTTCGTGAACTTCT  >  NC_007779/4169364‑4169559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.