Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,629,792 T→C 10.7% D467G (GAC→GGC)  guaA ← GMP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,629,7920TC10.7% 139.9 / 9.9 56D467G (GAC→GGC) guaAGMP synthetase
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (3/3);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCT  >  W3110S.gb/2629722‑2629855
                                                                      |                                                               
gTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGt                                                                 >  1:1029826/1‑71 (MQ=255)
gTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCAACAGCACGCAGAGAGACAACCCAGt                                                                 <  1:566083/71‑1 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCATAGAGACAACCCAGTCAt                                                              <  1:1910709/71‑1 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCa                                                               <  1:2358521/70‑1 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCa                                                               <  1:2384682/70‑1 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCAt                                                              <  1:1296509/71‑1 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCAt                                                              <  1:1459969/71‑1 (MQ=255)
      cAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACt                                                           <  1:1934143/71‑1 (MQ=255)
          tGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                        >  1:1219996/1‑70 (MQ=255)
          tGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                        >  1:1587725/1‑70 (MQ=255)
          tGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACg                                                       <  1:195023/71‑1 (MQ=255)
             gcCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGAc                                                     <  1:251197/70‑1 (MQ=255)
                  gtgtGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGc                                               >  1:302755/1‑71 (MQ=255)
                   tgtgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGcc                                              <  1:1533913/71‑1 (MQ=255)
                     tgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCa                                                               >  1:372403/1‑52 (MQ=255)
                     tgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCa                                            >  1:2384629/1‑71 (MQ=255)
                     tgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCa                                            >  1:890690/1‑71 (MQ=255)
                     tgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCa                                            >  1:2120690/1‑71 (MQ=255)
                       cGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACt                                                           >  1:180757/1‑54 (MQ=255)
                           cATAAAGTCGATGGTTTCGACAGCAAGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGc                                      >  1:2650012/1‑71 (MQ=255)
                            aTAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGc                                      >  1:1925432/1‑70 (MQ=255)
                              aaaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCa                                                               >  1:2393013/1‑43 (MQ=255)
                              aaaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCaa                                   >  1:1817389/1‑71 (MQ=255)
                               aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCa                                    >  1:2538284/1‑69 (MQ=255)
                               aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAAc                                  >  1:577127/1‑71 (MQ=255)
                               aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAAc                                  >  1:1096617/1‑71 (MQ=255)
                                 gTCTATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCaa                                   >  1:1483260/1‑68 (MQ=255)
                                 gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                                <  1:2375073/71‑1 (MQ=255)
                                     aTGGTTTCGACAGCACGCAGAGAGACAACCCAGCCAt                                                              <  1:2371265/37‑1 (MQ=255)
                                     aTGGTTTCGACAGCACGCAGAGAGACAACCCAGCCAt                                                              <  1:1182263/37‑1 (MQ=255)
                                     aTGGTTTCGACAGCACGC‑‑AGAGACAACCCAGTCATACTTACGACCATAGcc                                              <  1:2706168/51‑1 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACg                            >  1:1970999/1‑70 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                           <  1:2147340/71‑1 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                           <  1:2713278/71‑1 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                           <  1:1421548/71‑1 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                           <  1:2076111/71‑1 (MQ=255)
                                       ggTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                           <  1:1905270/70‑1 (MQ=255)
                                       ggTTTCGACAGCACGCAGAGAGACAACCCAGCCATAc                                                            >  1:494316/1‑37 (MQ=255)
                                       ggTTTCGACAGCACGCAGAGAGACAACCCAGCCATAc                                                            >  1:25655/1‑37 (MQ=255)
                                         tttCGACAGCACGCAGAGAGACAACCCAGTCATACTTACg                                                       <  1:699409/40‑1 (MQ=255)
                                          ttCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                        >  1:2393091/1‑38 (MQ=255)
                                          ttCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                        >  1:270/1‑38 (MQ=255)
                                          ttCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                        >  1:2812047/1‑38 (MQ=255)
                                                aGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGa                               <  1:1873393/57‑1 (MQ=255)
                                                aGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAgg                  >  1:2711474/1‑70 (MQ=255)
                                                   aCGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAaca              <  1:2575862/71‑1 (MQ=255)
                                                    cGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAgcg                 >  1:1786896/1‑65 (MQ=255)
                                                    cGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAacac             <  1:441775/71‑1 (MQ=255)
                                                        gagagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGt          <  1:2040206/70‑1 (MQ=255)
                                                          gagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTAc                         >  1:2735638/1‑53 (MQ=255)
                                                          gagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTAc                         >  1:2440732/1‑53 (MQ=255)
                                                          gagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGaa       >  1:1835511/1‑71 (MQ=255)
                                                           agaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAAc      <  1:916121/71‑1 (MQ=255)
                                                            gaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGaa                >  1:473375/1‑60 (MQ=255)
                                                                aCCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCt  <  1:201558/70‑1 (MQ=255)
                                                                     gTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg         <  1:2700008/58‑1 (MQ=255)
                                                                     gTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg         <  1:1909998/58‑1 (MQ=255)
                                                                     gTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg         <  1:575210/58‑1 (MQ=255)
                                                                      |                                                               
GTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCT  >  W3110S.gb/2629722‑2629855

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: ‑