Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,721,254 A→C 23.4% intergenic (+124/‑229) atpC → / → glmU F1 sector of membrane‑bound ATP synthase, epsilon subunit/fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,721,2540AC23.4% 82.3 / 47.9 70intergenic (+124/‑229)atpC/glmUF1 sector of membrane‑bound ATP synthase, epsilon subunit/fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base C (5/11);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

GTGTGACCCGTCCTGAATAGCGTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTAATTTTCTG  >  W3110S.gb/3721187‑3721323
                                                                   |                                                                     
gtgtGACCCGTCCTGAATAGCGTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTAtt                                                                     >  1:760983/1‑70 (MQ=255)
 tgtgACCCGTCCTGAATAGCGTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCa                                                                   <  1:1495618/71‑1 (MQ=255)
          tCCTGAATAGCGTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGaa                                                          <  1:1067227/71‑1 (MQ=255)
            cTGAATAGCGTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGaaa                                                         >  1:3200877/1‑70 (MQ=255)
               aaTAGCGTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATaa                                                      <  1:3454509/70‑1 (MQ=255)
                    cGTTCACATAGATCCTGCTGATATAAAACCACCCTGTTTTCCTGTTTATTCATTGATCGAAATaa                                                      >  1:2147202/1‑65 (MQ=255)
                     gTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATaa                                                      <  1:2886120/64‑1 (MQ=255)
                     gTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCaaa                                               <  1:222809/71‑1 (MQ=255)
                     gTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCaaa                                               <  1:2735123/71‑1 (MQ=255)
                      ttCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCaaa                                               <  1:1083152/70‑1 (MQ=255)
                         acaTAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAAc                                            <  1:2499438/70‑1 (MQ=255)
                          caTAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTAtt                                                                     <  1:757309/44‑1 (MQ=255)
                            tAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATcc                                        <  1:2300016/71‑1 (MQ=255)
                            tAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATc                                         >  1:3243877/1‑70 (MQ=255)
                             aGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCa                                                                   >  1:3612018/1‑43 (MQ=255)
                               aTCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCa                                                 <  1:3251117/59‑1 (MQ=255)
                                 ccTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCt                                    >  1:105708/1‑70 (MQ=255)
                                  ctgctgATATAAAACCCCCCTGTTTTCCTGTTTAtt                                                                     <  1:3213723/36‑1 (MQ=255)
                                  ctgctgATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTg                                   <  1:268101/70‑1 (MQ=255)
                                     ctgATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGc                               >  1:1985888/1‑71 (MQ=255)
                                       gATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTg                                   >  1:1053758/1‑65 (MQ=255)
                                       gATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCt                              >  1:3066788/1‑70 (MQ=255)
                                       gATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCt                              >  1:1140953/1‑70 (MQ=255)
                                        atatAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCt                              >  1:2487464/1‑69 (MQ=255)
                                        atatAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTa                            <  1:1929079/71‑1 (MQ=255)
                                         tataAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTa                            <  1:3545570/70‑1 (MQ=255)
                                           taAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAAt                                                        <  1:568005/40‑1 (MQ=255)
                                           taAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCaaaa                                              <  1:853500/50‑1 (MQ=255)
                                           taAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTa                            <  1:1847006/68‑1 (MQ=255)
                                             aaaCCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATaa                                                      <  1:3157478/40‑1 (MQ=255)
                                             aaaCCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTaaa                          >  1:1655298/1‑68 (MQ=255)
                                             aaaCCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAAtt                        <  1:2113947/70‑1 (MQ=255)
                                             aaaCCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAAtt                        >  1:3630826/1‑70 (MQ=255)
                                              aaCCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTaa                      <  1:159033/71‑1 (MQ=255)
                                               aCCCCCCTGTTTTCCTGTTTATTCATTGATCGaaa                                                         <  1:3466849/35‑1 (MQ=255)
                                                  ccccTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCaaaa                                              <  1:3663074/43‑1 (MQ=255)
                                                     ttgTTTTCCTGTTTCTTtttcaatgatgtcagataaat                                                <  1:1048223‑M1/37‑22 (MQ=255)
                                                     ttgTTTTCCTGTTTCTTtttcaatgatgtcagataaat                                                <  1:2359972‑M1/37‑22 (MQ=255)
                                                     ttgTTTTCCTGTTTCTTtttcaatgatgtcagataaat                                                <  1:1167964‑M1/37‑22 (MQ=255)
                                                     ttgTTTTCCTGTTTCTTtttcaatgatgtcagataaat                                                <  1:727547‑M1/37‑22 (MQ=255)
                                                     ttgTTTTCCTGTTTCTTtttcaatgatgtcagataaa                                                 <  1:3181291‑M1/36‑21 (MQ=255)
                                                     ttgTTTTCCTGTTTCTTtttcaatgatgtcagataaa                                                 <  1:1507670‑M1/36‑21 (MQ=255)
                                                     ttgTTTTCCTGTTTCTTtttcaatgatgtcagataa                                                  <  1:2168559‑M1/35‑20 (MQ=255)
                                                     cTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACa                                           <  1:3557049/43‑1 (MQ=255)
                                                      tGTTTTCCTGTTTCTTtttcaatgatgtcagataaatt                                               >  1:1522072‑M1/1‑16 (MQ=255)
                                                      tGTTTTCCTGTTTCTTtttcaatgatgtcagataaat                                                <  1:1705872‑M1/37‑22 (MQ=255)
                                                      tGTTTTCCTGTTTCTTtttcaatgatgtcagataaa                                                 <  1:2079721‑M1/36‑21 (MQ=255)
                                                      tGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAAc                                            >  1:2814279/1‑41 (MQ=255)
                                                      tGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACAt                                          >  1:3259934/1‑43 (MQ=255)
                                                      tGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAgg                    <  1:1316659/65‑1 (MQ=255)
                                                      tGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAgg                    <  1:1557290/65‑1 (MQ=255)
                                                      tGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTa                  >  1:2732782/1‑67 (MQ=255)
                                                       gTTTTCCTGTTTGTTCATTGTTCGAAATAAGAGCa                                                 <  1:2644207/35‑1 (MQ=38)
                                                       gTTTTCCTGTTTCTTtttcaatgatgtcagataaatt                                               >  1:1301052‑M1/1‑15 (MQ=255)
                                                       gTTTTCCTGTTTCTTtttcaatgatgtcagataaatt                                               >  1:2456335‑M1/1‑15 (MQ=255)
                                                       gTTTTCCTGTTTCTTtttcaatgatgtcagataaat                                                >  1:1413530‑M1/1‑15 (MQ=255)
                                                       gTTTTCCTGTTTCTTtttcaatgatgtcagataaat                                                <  1:625902‑M1/36‑22 (MQ=255)
                                                       gTTTTCCTGTTTCTTtttcaatgatgtcagataaa                                                 <  1:2150284‑M1/35‑21 (MQ=255)
                                                        ttttCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGa                                  >  1:2593478/1‑49 (MQ=255)
                                                        ttttCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTaaa                          >  1:729760/1‑57 (MQ=255)
                                                        ttttCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCt            >  1:3081204/1‑71 (MQ=255)
                                                         tttCCTGTTTCTTtttcaatgatgtcagataaatt                                               >  1:3048882‑M1/1‑13 (MQ=255)
                                                          ttCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCt                                    >  1:3116730/1‑45 (MQ=255)
                                                          ttCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACg                                >  1:1383417/1‑49 (MQ=255)
                                                           tCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTa          <  1:2451832/70‑1 (MQ=255)
                                                            ccTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTa          <  1:1639586/69‑1 (MQ=255)
                                                                tttATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTAAtttt     >  1:2965364/1‑70 (MQ=255)
                                                                 ttATTCATTGATCGAAATAAGAGCAAAAACATCCACCTg                                   <  1:3203284/39‑1 (MQ=255)
                                                                  tATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTAAt        >  1:2549391/1‑65 (MQ=255)
                                                                  tATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTAATTTTCTg  >  1:1434404/1‑71 (MQ=255)
                                                                   aTTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTAATTTTCTg  >  1:3706550/1‑70 (MQ=255)
                                                                   |                                                                     
GTGTGACCCGTCCTGAATAGCGTTCACATAGATCCTGCTGATATAAAACCCCCCTGTTTTCCTGTTTATTCATTGATCGAAATAAGAGCAAAAACATCCACCTGACGCTTAAATTAAGGTACTGCCTTAATTTTCTG  >  W3110S.gb/3721187‑3721323

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: