Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 31,629 G→A 7.7% V271V (GTG→GTA carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb31,6290GA 7.7% 149.8 / 2.7 52V271V (GTG→GTAcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (20/28);  new base A (2/2);  total (22/30)
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.50e-01

CACTGTCGCGCCAGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGAAAAACGG  >  W3110S.gb/31560‑31694
                                                                     |                                                                 
cACTGTCGCGCCAGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGtgc                                                                  >  1:1950511/1‑71 (MQ=255)
 aCTGTCGCGCCAGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGtgct                                                                 >  1:402944/1‑71 (MQ=255)
 aCTGTCGCGCCAGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGtgct                                                                 >  1:3238376/1‑71 (MQ=255)
 aCTGTCGCGCCAGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGtgct                                                                 >  1:1905451/1‑71 (MQ=255)
 aCTGTCGCGCCAGCCCAAACGCTGAACGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCTGtgct                                                                 >  1:472458/1‑71 (MQ=255)
            aGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCg                                                      <  1:2095277/71‑1 (MQ=255)
            aGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCg                                                      <  1:1939528/71‑1 (MQ=255)
                    cGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGa                                               >  1:3797953/1‑70 (MQ=255)
                       tGACCGACAAAGAATACCAAATCATGCGTAACGCCTCGATGGCGGTACTGCGTGAAATCGGCGTTGAAAcc                                           >  1:3350095/1‑71 (MQ=255)
                       tGACCGACAAAGAATACCAAATCATGCGTAACGCCTCGATGGCGGTACTGCGTGAAATCGGCGTTGAAAc                                            >  1:815649/1‑70 (MQ=255)
                         aCCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAAt                                                        >  1:149175/1‑56 (MQ=255)
                         aCCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAAcc                                           <  1:2595103/69‑1 (MQ=255)
                         aCCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAAcc                                           <  1:3571040/69‑1 (MQ=255)
                         aCCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCgg                                         <  1:544484/71‑1 (MQ=255)
                         aCCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCgg                                         <  1:1164278/71‑1 (MQ=255)
                         aCCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCgg                                         <  1:910766/71‑1 (MQ=255)
                          ccGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGc                                                    <  1:1326645/59‑1 (MQ=255)
                          ccGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGc                                                    <  1:2417062/59‑1 (MQ=255)
                             aCAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGaa                                                          <  1:1025437/50‑1 (MQ=255)
                             aCAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGa                                               <  1:792218/61‑1 (MQ=255)
                             aCAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGa                                               <  1:1138442/61‑1 (MQ=255)
                             aCAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGa                                               <  1:2660392/61‑1 (MQ=255)
                              cAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGtt                                    >  1:3663819/1‑71 (MQ=255)
                                aaGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGaaa                                                         <  1:3192192/48‑1 (MQ=255)
                                aaGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTcc                                  <  1:2808902/71‑1 (MQ=255)
                                aaGAATATCAAATCATGCGTAACGCCTCGATGGCGGAGCTGCGTGAAATCGGGGt                                                  >  1:929839/1‑55 (MQ=255)
                                 aGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTcc                                  >  1:2221457/1‑70 (MQ=255)
                                  gAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCaa                                >  1:2968158/1‑71 (MQ=255)
                                     tatCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGaa                                                          <  1:2456956/42‑1 (MQ=255)
                                       tCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTcc                                  >  1:3451024/1‑64 (MQ=255)
                                       tCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTcc                                  >  1:1983631/1‑64 (MQ=255)
                                       tCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTcc                                  >  1:1439093/1‑64 (MQ=255)
                                       tCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTcc                                  >  1:3138702/1‑64 (MQ=255)
                                       tCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGtt                            >  1:2258090/1‑70 (MQ=255)
                                       tCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGtt                            <  1:120012/70‑1 (MQ=255)
                                         aaaTCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAg                         >  1:266629/1‑71 (MQ=255)
                                         aaaTCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAg                         >  1:1842774/1‑71 (MQ=255)
                                           aTCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAAc                                            <  1:3057209/50‑1 (MQ=255)
                                             cATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCg                                          >  1:3553024/1‑50 (MQ=255)
                                             cATGCGTAACGCCTCGATGGCGGTGCTGAGTGAAATCGGCGTTGAAACCg                                          >  1:3194210/1‑50 (MQ=255)
                                              aTGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGc                    <  1:1364428/71‑1 (MQ=255)
                                                  gTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTg                <  1:876262/71‑1 (MQ=255)
                                                  gTAACGCCTCGATGGCGGTACTGCGTGAAATCGGCGTTGAAACCggtgg                                      <  1:113222/49‑1 (MQ=255)
                                                  gTAACGCCTCGATGGCGGTACTGCGTGAAATCGGCGTTGAAACCggtgg                                      <  1:1452545/49‑1 (MQ=255)
                                                     aCGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTcc                                  >  1:874184/1‑50 (MQ=255)
                                                     aCGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAAc             >  1:715683/1‑71 (MQ=255)
                                                         cTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGtt                       <  1:356221/57‑1 (MQ=255)
                                                         cTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGtt                       <  1:1198391/57‑1 (MQ=255)
                                                         cTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCg          <  1:461432/70‑1 (MQ=255)
                                                           cGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTg                     >  1:1098582/1‑57 (MQ=255)
                                                           cGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGaa        >  1:2350763/1‑70 (MQ=255)
                                                            gATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGaaaa      <  1:225666/71‑1 (MQ=255)
                                                            gATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGaaa       <  1:633013/70‑1 (MQ=255)
                                                            gATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGaaa       <  1:2406438/70‑1 (MQ=255)
                                                               ggCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGAAAAACg   >  1:2913725/1‑71 (MQ=255)
                                                                gCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGAAAAACgg  <  1:2289503/71‑1 (MQ=255)
                                                                gCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGAAAAACgg  <  1:2702159/71‑1 (MQ=255)
                                                                     |                                                                 
CACTGTCGCGCCAGCCCAAACGCTGACCGACAAAGAATATCAAATCATGCGTAACGCCTCGATGGCGGTGCTGCGTGAAATCGGCGTTGAAACCGGTGGTTCCAACGTTCAGTTTGCGGTGAACCCGAAAAACGG  >  W3110S.gb/31560‑31694

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: