Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,987,171 G→A 38.9% G12G (GGC→GGT araG ← fused L‑arabinose transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,987,1710GA38.9% 42.0 / 58.0 54G12G (GGC→GGTaraGfused L‑arabinose transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (22/11);  new base A (7/14);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 2.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

AACCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACAC  >  W3110S.gb/1987102‑1987237
                                                                     |                                                                  
aaCCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGc                                                                   >  1:5045028/1‑71 (MQ=255)
aaCCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGc                                                                   <  1:5057999/71‑1 (MQ=255)
         ggCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATg                                                          <  1:3104073/71‑1 (MQ=255)
          gCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAAt                                                           >  1:1876065/1‑69 (MQ=255)
                  gTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACg                                                 <  1:4185835/71‑1 (MQ=255)
                   tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGTAAACGTTTTACCGATGCCGCGAAATGAGAGATa                                                   <  1:4389947/68‑1 (MQ=255)
                   tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGaa                                                             <  1:4189541/58‑1 (MQ=255)
                     aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACgg                                                >  1:3840712/1‑69 (MQ=255)
                            atatCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTg                                       <  1:4165907/71‑1 (MQ=255)
                              atCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGtt                                     >  1:2124216/1‑71 (MQ=255)
                                 cGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATgaga                                                       >  1:3333434/1‑50 (MQ=255)
                                 cGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTg                                    >  1:1975091/1‑69 (MQ=255)
                                 cGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCa                                  >  1:1764931/1‑71 (MQ=255)
                                 cGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCa                                  >  1:3136645/1‑71 (MQ=255)
                                  gTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCAt                                 <  1:1676583/71‑1 (MQ=255)
                                  gTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCAt                                 <  1:3770383/71‑1 (MQ=255)
                                  gTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCAt                                 <  1:1462949/71‑1 (MQ=255)
                                   tCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATgagaga                                                     <  1:1359202/50‑1 (MQ=255)
                                   tCAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAAt                                                           <  1:1182185/44‑1 (MQ=255)
                                    cAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATgagaga                                                     <  1:522683/49‑1 (MQ=255)
                                    cAGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAAt                                                           <  1:1924112/43‑1 (MQ=255)
                                     aGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAAc                              >  1:1889302/1‑71 (MQ=255)
                                     aGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGAGAAATGAGAGATACGGGGTAGACTGTTGCATAAc                              >  1:4722623/1‑71 (MQ=255)
                                     aGCGCCTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAAc                              <  1:2621097/71‑1 (MQ=255)
                                          cTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGa                                          >  1:4682559/1‑54 (MQ=255)
                                          cTTAACGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGtctc                         >  1:66097/1‑71 (MQ=255)
                                           ttAACGCCGGGAAACGTTTTACCGATGCCGCGAAATgaga                                                       >  1:1227761/1‑40 (MQ=255)
                                           ttAACGCCGGGAAACGTTTTACCGATACCGCGAAATgagaga                                                     <  1:4989634/42‑1 (MQ=255)
                                            tAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGtt                                     <  1:4105335/57‑1 (MQ=255)
                                            tAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCAt                                 >  1:2533345/1‑61 (MQ=255)
                                            tAACGCCGGGAAACGTTTTACCGATACCGCGAAATgagaga                                                     <  1:70651/41‑1 (MQ=255)
                                            tAACGCCGGGAAACGTTTTACCGATACCGCGAAATgagaga                                                     <  1:1247949/41‑1 (MQ=255)
                                             aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGa                                          >  1:982767/1‑51 (MQ=255)
                                             aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGt                      >  1:3592022/1‑71 (MQ=255)
                                             aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGt                      >  1:3150035/1‑71 (MQ=255)
                                             aaCGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATa                                                   >  1:806751/1‑42 (MQ=255)
                                              aCGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGt                      <  1:691716/70‑1 (MQ=255)
                                              aCGCCGGGAAACGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTg                     <  1:2316720/71‑1 (MQ=255)
                                                 ccGGGAAACGTTTTACCGATGCCGCGAAATgagaga                                                     >  1:4626659/1‑36 (MQ=255)
                                                 ccGGGAAACGTTTTACCGATGCCGCGAAATgagaga                                                     >  1:3133045/1‑36 (MQ=255)
                                                      aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTg             <  1:4620676/71‑1 (MQ=255)
                                                      aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTACGTGaa                   >  1:2889192/1‑65 (MQ=255)
                                                         cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCAc               >  1:4007380/1‑66 (MQ=255)
                                                         cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCAc               >  1:163495/1‑66 (MQ=255)
                                                         cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCAc               >  1:183344/1‑66 (MQ=255)
                                                         cGTTTTACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTc          >  1:1618657/1‑71 (MQ=255)
                                                            tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGtt       >  1:5043832/1‑71 (MQ=255)
                                                            tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGtt       >  1:4323567/1‑71 (MQ=255)
                                                             ttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCg         <  1:466228/68‑1 (MQ=255)
                                                             ttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCg         <  1:353706/68‑1 (MQ=255)
                                                             ttACCGATGCCGAGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGtt       <  1:4070558/70‑1 (MQ=255)
                                                             ttACCGATACCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGtt       >  1:4310897/1‑70 (MQ=255)
                                                                 cGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCacac  >  1:2033262/1‑71 (MQ=255)
                                                                  gATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGtt       <  1:4537507/65‑1 (MQ=255)
                                                                     gCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTg                     <  1:2821195/48‑1 (MQ=255)
                                                                     |                                                                  
AACCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACAC  >  W3110S.gb/1987102‑1987237

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