Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,629,792 T→C 27.6% 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,7920TC27.6% 115.4 / 56.3 76D467G (GAC→GGC) guaAGMP synthetase
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base C (8/13);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

GTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCTGGCTGAC  >  W3110S.gb/2629722‑2629862
                                                                      |                                                                      
gTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGt                                                                        <  1:2738335/71‑1 (MQ=255)
gTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGCTTTCGACAGCACCCAGAGAGACAACCCAGt                                                                        >  1:4315257/1‑71 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCa                                                                      <  1:2867101/70‑1 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCAt                                                                     <  1:3527571/71‑1 (MQ=255)
   cGGCAGATGCGCCCAGTGTGCGGTCATAAAGCCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCAt                                                                     <  1:2841793/71‑1 (MQ=255)
       aGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACtt                                                                 <  1:1511436/71‑1 (MQ=255)
       aGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACt                                                                  >  1:1135949/1‑70 (MQ=255)
         aTGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTAc                                                               <  1:4174217/71‑1 (MQ=255)
          tGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                               >  1:2945542/1‑70 (MQ=255)
          tGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                               >  1:138968/1‑70 (MQ=255)
            cgcCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGAc                                                            >  1:2164192/1‑71 (MQ=255)
             gcCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGAcc                                                           <  1:153947/71‑1 (MQ=255)
               ccAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCAt                                                                     <  1:1269018/59‑1 (MQ=255)
               ccAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCAt                                                                     <  1:385507/59‑1 (MQ=255)
                  gtgtGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGc                                                      <  1:1883674/71‑1 (MQ=255)
                     tgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACt                                                                  >  1:51094/1‑56 (MQ=255)
                       cGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCAt                                                  <  1:3052416/70‑1 (MQ=255)
                       cGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCAt                                                  >  1:1089939/1‑70 (MQ=255)
                       cGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCAt                                                  <  1:1193715/70‑1 (MQ=255)
                         gTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATAc                                                                   >  1:2428000/1‑51 (MQ=255)
                         gTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTAc                                               >  1:1879649/1‑71 (MQ=255)
                          tCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACt                                                                  <  1:4187055/51‑1 (MQ=255)
                          tCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCAt                                                  <  1:2590486/67‑1 (MQ=255)
                           cATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTa                                                >  1:1966915/1‑68 (MQ=255)
                           cATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGc                                             >  1:3115664/1‑71 (MQ=255)
                            aTAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACt                                                                  <  1:3473253/49‑1 (MQ=255)
                            aTAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGc                                             >  1:2469799/1‑70 (MQ=255)
                               aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAAc                                         >  1:4090520/1‑71 (MQ=255)
                               aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAAc                                         >  1:4132649/1‑71 (MQ=255)
                               aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAAc                                         >  1:4138044/1‑71 (MQ=255)
                                 gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                                       <  1:3665922/71‑1 (MQ=255)
                                  tCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                                       >  1:574413/1‑70 (MQ=255)
                                  tCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGAcc                                                           <  1:3484535/50‑1 (MQ=255)
                                  tCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGAcc                                                           <  1:3807205/50‑1 (MQ=255)
                                  tCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCAt                                                         >  1:2621322/1‑52 (MQ=255)
                                   cGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                               >  1:3851775/1‑45 (MQ=255)
                                   cGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                               >  1:841780/1‑45 (MQ=255)
                                   cGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTAc                                                               >  1:3265287/1‑45 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                                  <  1:3102453/71‑1 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                                  <  1:97328/71‑1 (MQ=255)
                                      tGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                                  <  1:916742/71‑1 (MQ=255)
                                       ggTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTa                                 >  1:913363/1‑71 (MQ=255)
                                       ggTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTa                                 >  1:2389551/1‑71 (MQ=255)
                                        gTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                                  <  1:3179208/69‑1 (MQ=255)
                                          ttCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCAt                                                  >  1:2329795/1‑51 (MQ=255)
                                          ttCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGc                                             >  1:4264128/1‑56 (MQ=255)
                                                aGCACGCAGAGAGACAACCCAGTCATACTTACGATCATCGCCCAt                                                  <  1:4334498/45‑1 (MQ=255)
                                                 gCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGaa                       <  1:641004/71‑1 (MQ=255)
                                                 gCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGaa                       <  1:2471501/71‑1 (MQ=255)
                                                   aCGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAaca                     >  1:1562796/1‑71 (MQ=255)
                                                    cGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAacac                    <  1:2581379/71‑1 (MQ=255)
                                                    cGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAacac                    <  1:2376701/71‑1 (MQ=255)
                                                     gCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAaca                     >  1:55500/1‑69 (MQ=255)
                                                       agagagACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                                       <  1:579704/49‑1 (MQ=255)
                                                       agagagACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                                       <  1:1727893/49‑1 (MQ=255)
                                                        gagagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGaa                       <  1:93684/64‑1 (MQ=255)
                                                        gagagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg                <  1:3579251/71‑1 (MQ=255)
                                                        gagagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg                <  1:242127/71‑1 (MQ=255)
                                                        gagagaCAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACTGAACGTACCGGCAGGAacac                    <  1:38508/67‑1 (MQ=255)
                                                         agagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGaa                                     >  1:4525429/1‑49 (MQ=255)
                                                         agagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGa               <  1:1468373/71‑1 (MQ=255)
                                                          gagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTAc                                >  1:2792638/1‑53 (MQ=255)
                                                          gagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTAc                                >  1:2269893/1‑53 (MQ=255)
                                                          gagaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAaca                     >  1:1340829/1‑64 (MQ=255)
                                                           agaCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAAc             <  1:3723011/71‑1 (MQ=255)
                                                             aCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACTGAACGTACCGGCAGGAaca                     <  1:404214/61‑1 (MQ=255)
                                                             aCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAacac                    <  1:3384239/62‑1 (MQ=255)
                                                             aCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAaca                     <  1:449620/61‑1 (MQ=255)
                                                             aCAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAaca                     <  1:901057/61‑1 (MQ=255)
                                                             aCAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGaa                       >  1:3435116/1‑59 (MQ=255)
                                                             aCAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAaca                     <  1:3308536/61‑1 (MQ=255)
                                                              cAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACTcc          <  1:857711/71‑1 (MQ=255)
                                                                   cAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAAcaca                   >  1:1556534/1‑57 (MQ=255)
                                                                   cAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCTGGCt     >  1:1538906/1‑71 (MQ=255)
                                                                    aGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTAc                                >  1:2495943/1‑43 (MQ=255)
                                                                    aGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCg                              <  1:3996180/45‑1 (MQ=255)
                                                                     gTCATACTTACGACCATCGCCCATTACGCCAACGGAAc                                    >  1:480715/1‑38 (MQ=255)
                                                                     gTCATACTTACGACCATCGCCCATTACGCCAACGGAAc                                    >  1:2218138/1‑38 (MQ=255)
                                                                      tCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGc                            >  1:340575/1‑45 (MQ=255)
                                                                      tCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCTGGCTGAc  <  1:90140/71‑1 (MQ=255)
                                                                      tCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCTGGCTGAc  <  1:1803434/71‑1 (MQ=255)
                                                                      |                                                                      
GTACGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTGAACGCCTGGCTGAC  >  W3110S.gb/2629722‑2629862

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: