Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,096,084 G→A 17.5% E160E (GAG→GAA hisB → fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,096,0840GA17.5% 237.9 / 38.6 97E160E (GAG→GAAhisBfused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (44/36);  new base A (10/7);  total (54/43)
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.84e-01

GCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGG  >  W3110S.gb/2096016‑2096154
                                                                    |                                                                      
gCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTTAACTGGCCAATGATTGGCGAGc                                                                       >  1:1508763/1‑70 (MQ=255)
gCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGc                                                                       <  1:3101950/70‑1 (MQ=255)
 cGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCaa                                                                     <  1:2855509/71‑1 (MQ=255)
   gAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACt                                                                   <  1:934243/71‑1 (MQ=255)
   gAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACt                                                                   <  1:3694052/71‑1 (MQ=255)
    aaaaCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTc                                                                  <  1:293769/71‑1 (MQ=255)
    aaaaCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTc                                                                  <  1:3626748/71‑1 (MQ=255)
            ggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCCCCAGACg                                                          >  1:3353941/1‑71 (MQ=255)
            ggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACg                                                          >  1:3846525/1‑71 (MQ=255)
            ggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACg                                                          >  1:1501721/1‑71 (MQ=255)
              cATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTg                                                        <  1:2938738/71‑1 (MQ=255)
              cATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTg                                                        <  1:2046374/71‑1 (MQ=255)
              cATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTg                                                        <  1:547055/71‑1 (MQ=255)
                ttACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGc                                                                       <  1:2704535/54‑1 (MQ=255)
                ttACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTc                                                                  <  1:1127720/59‑1 (MQ=255)
                ttACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTc                                                                  <  1:873724/59‑1 (MQ=255)
                 tACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAcc                                                     >  1:4017690/1‑71 (MQ=255)
                 tACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAcc                                                     >  1:1559865/1‑71 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCaa                                                                     >  1:3873949/1‑54 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGa                                                            >  1:4006272/1‑63 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGa                                                            >  1:452580/1‑63 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCg                                                    >  1:679025/1‑71 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCg                                                    >  1:2499935/1‑71 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCg                                                    >  1:2610890/1‑71 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCg                                                    >  1:3681184/1‑71 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCg                                                    >  1:3435335/1‑71 (MQ=255)
                  aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCg                                                    >  1:771564/1‑71 (MQ=255)
                  aCTGGATTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCa                                                              >  1:660422/1‑61 (MQ=255)
                   cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                   <  1:1286044/71‑1 (MQ=255)
                   cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                   <  1:866611/71‑1 (MQ=255)
                   cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                   <  1:3794796/71‑1 (MQ=255)
                   cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                   <  1:2537943/71‑1 (MQ=255)
                   cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                   <  1:3918839/71‑1 (MQ=255)
                   cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACAGt                                                   <  1:2229832/71‑1 (MQ=255)
                         tACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                             >  1:2719082/1‑71 (MQ=255)
                         tACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                             >  1:949044/1‑71 (MQ=255)
                            gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCAc                                          <  1:3273273/71‑1 (MQ=255)
                            gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCAc                                          <  1:2892866/71‑1 (MQ=255)
                             cTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         <  1:3292073/71‑1 (MQ=255)
                              tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                   >  1:1128969/1‑60 (MQ=255)
                              tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         >  1:2102778/1‑70 (MQ=255)
                              tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         >  1:1585634/1‑70 (MQ=255)
                                cGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGt                                        <  1:3337144/69‑1 (MQ=255)
                                   ccGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGc                                   >  1:2133557/1‑71 (MQ=255)
                                   ccGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGc                                   >  1:3285015/1‑71 (MQ=255)
                                     gcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                                 <  1:4065123/71‑1 (MQ=255)
                                     gcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                                 <  1:269000/71‑1 (MQ=255)
                                     gcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                                 <  1:3774575/71‑1 (MQ=255)
                                     gcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                                 <  1:3855146/71‑1 (MQ=255)
                                             cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                          >  1:3669482/1‑70 (MQ=255)
                                             cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                          >  1:1650420/1‑70 (MQ=255)
                                             cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                          >  1:59784/1‑70 (MQ=255)
                                               gAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                          >  1:559101/1‑68 (MQ=255)
                                                aaCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                               <  1:2804414/62‑1 (MQ=255)
                                                aaCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                               <  1:1481263/62‑1 (MQ=255)
                                                aaCTGGCCAATGATTGGAGAGCAACTCACCAGACGTGACCGTTACGCACACGTAGTGCGTAATACCAAaga                      <  1:159220/71‑1 (MQ=255)
                                                  cTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATa                             >  1:138061/1‑62 (MQ=255)
                                                  cTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATa                             >  1:2233017/1‑62 (MQ=255)
                                                    ggCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACg                  >  1:3133475/1‑71 (MQ=255)
                                                    ggCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACg                  >  1:128397/1‑71 (MQ=255)
                                                    ggCCAATGATTGGCGAGCAACTAACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACg                  >  1:1099058/1‑71 (MQ=255)
                                                      ccAATGATTGGCGAGCAACTCACCAGACGTGACCGtt                                                  <  1:1886898/37‑1 (MQ=255)
                                                      ccAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                      <  1:128264/65‑1 (MQ=255)
                                                        aaTGATTGGCGAGCAACTCACCAGACGTGACCGTTAc                                                >  1:3660851/1‑37 (MQ=255)
                                                        aaTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACg                  >  1:4022449/1‑67 (MQ=255)
                                                        aaTGATTGGCGAGCAACTCACCAGACGTGACAGTTAc                                                >  1:1862070/1‑37 (MQ=255)
                                                         aTGATTGGCGAGCAAGTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCa                <  1:4035262/68‑1 (MQ=255)
                                                         aTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTATTGCGt                                 <  1:1661626/51‑1 (MQ=255)
                                                         aTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCa                <  1:2452570/68‑1 (MQ=255)
                                                          tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         <  1:76054/42‑1 (MQ=255)
                                                          tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         <  1:753584/42‑1 (MQ=255)
                                                          tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         <  1:276108/42‑1 (MQ=255)
                                                          tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         <  1:1334278/42‑1 (MQ=255)
                                                          cgATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                         <  1:3983846/41‑1 (MQ=255)
                                                           gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt            >  1:605710/1‑70 (MQ=255)
                                                           gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt            >  1:2948666/1‑70 (MQ=255)
                                                           gATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt            >  1:2571198/1‑70 (MQ=255)
                                                           gATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt            >  1:704206/1‑70 (MQ=255)
                                                           gATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt            >  1:3972681/1‑70 (MQ=255)
                                                           gATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt            >  1:3403342/1‑70 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTa                                       >  1:1343146/1‑42 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAgaga                    >  1:3921756/1‑61 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAgaga                    >  1:1992085/1‑61 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTg           <  1:381081/70‑1 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa          >  1:1319771/1‑71 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa          >  1:1849906/1‑71 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa          >  1:2260572/1‑71 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa          >  1:2034704/1‑71 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa          >  1:3268712/1‑71 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa          >  1:3301453/1‑71 (MQ=255)
                                                            aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa          >  1:1941921/1‑71 (MQ=255)
                                                              tGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGAc         >  1:530640/1‑70 (MQ=255)
                                                              tGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACg        >  1:1713279/1‑71 (MQ=255)
                                                               ggCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGt       >  1:1455695/1‑71 (MQ=255)
                                                                 cGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc     <  1:2313776/71‑1 (MQ=255)
                                                                 cGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGAAGCAGATTGACGTTc     <  1:250768/71‑1 (MQ=255)
                                                                 cGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGAAGCAGATTGACGTTc     <  1:245969/71‑1 (MQ=255)
                                                                    gCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg  <  1:579718/71‑1 (MQ=255)
                                                                    gCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg  <  1:3063202/71‑1 (MQ=255)
                                                                    gCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg  <  1:1219482/71‑1 (MQ=255)
                                                                    |                                                                      
GCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGG  >  W3110S.gb/2096016‑2096154

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: