Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 407,009 G→A 40.6% I157I (ATC→ATT cobC ← predicted alpha‑ribazole‑5'‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE407,0090GA40.6% 51.7 / 79.3 113I157I (ATC→ATTcobCpredicted alpha‑ribazole‑5'‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base A (22/24);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 113 total aligned reads displayed.

ACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATT  >  minE/406951‑407070
                                                          |                                                             
acCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAg                                                       >  1:1726358/1‑67 (MQ=255)
acCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGa                                                      >  1:2766064/1‑68 (MQ=255)
acCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGa                                                      >  1:2333211/1‑68 (MQ=255)
acCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGa                                                      >  1:396371/1‑68 (MQ=255)
                     gCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAg                                                 <  1:292349/52‑1 (MQ=255)
                     gCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAg                                                 <  1:623454/52‑1 (MQ=255)
                     gCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAg                                                 <  1:549527/52‑1 (MQ=255)
                     gCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAg                                                 <  1:859851/52‑1 (MQ=255)
                     gCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAg                                                 <  1:2418242/52‑1 (MQ=255)
                     gCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACCCAg                                                 <  1:2071864/52‑1 (MQ=255)
                           tAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACt                           >  1:1371116/1‑68 (MQ=255)
                           tAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACt                           >  1:560795/1‑68 (MQ=255)
                           tAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACt                           >  1:519805/1‑68 (MQ=255)
                           tAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACt                           >  1:2510911/1‑68 (MQ=255)
                           tAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACt                           >  1:2646447/1‑68 (MQ=255)
                           tAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACt                           >  1:3385859/1‑68 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        >  1:242533/1‑37 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:3124142/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        >  1:2095119/1‑37 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:3310533/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:2720847/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:1769838/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:2861133/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:306620/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:1382688/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                        <  1:162836/37‑1 (MQ=255)
                             gCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTaa                         >  1:2765503/1‑68 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:399042/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1678678/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:2559321/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:2578149/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1271297/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:358646/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1258008/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:938751/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1135069/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:2796682/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:2912998/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1098811/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:3194017/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:2985157/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1041258/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:471910/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:829934/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1794615/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:1650070/65‑1 (MQ=255)
                                  ggcaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:2099704/65‑1 (MQ=255)
                                  ggaaggcaTGCCAATTAAACGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTaaaa                       <  1:836437/62‑1 (MQ=255)
                                        caTGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATa             >  1:2359559/1‑69 (MQ=255)
                                        caTGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATa             >  1:2455767/1‑69 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAata                     <  1:91222/57‑1 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAata                     <  1:1329119/57‑1 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAata                     <  1:2294701/57‑1 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAata                     <  1:755514/57‑1 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATATTg          >  1:2173334/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:3295990/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:2757516/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:1823879/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:1746463/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:2999710/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:1721295/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:3086933/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:1107516/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:630605/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:404918/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:1208036/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:2102005/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:1441811/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg          >  1:2639420/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTATGATAGTg          >  1:2632840/1‑68 (MQ=255)
                                            ccAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAAAATTCTGATAGTg          >  1:3188942/1‑68 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:3271966/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:633676/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:1054528/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:31021/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:1330008/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:133575/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:1478647/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:1599197/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:1608431/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:1671302/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:1806634/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:2287817/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:2363301/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:2376960/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:2508584/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:2701679/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:2895377/1‑66 (MQ=255)
                                                  aaaCGAGCAATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGata    >  1:3184150/1‑66 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTTCTGAAAtt  <  1:2954273/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:2369105/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:2673623/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:1870446/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:3401631/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:291201/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:1420647/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:1374682/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:1194929/68‑1 (MQ=255)
                                                    aCGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:951920/68‑1 (MQ=255)
                                                     cGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAAtt  <  1:2677646/67‑1 (MQ=255)
                                                          |                                                             
ACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATT  >  minE/406951‑407070

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