Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 123,736 A→G 11.4% K240K (AAA→AAG aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb123,7360AG11.4% 184.1 / 13.8 70K240K (AAA→AAGaceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base A (25/37);  new base G (3/5);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGATACCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATCAACTGTAACCTGCAG  >  W3110S.gb/123668‑123805
                                                                    |                                                                     
aaaGATACCTCTAAACAAACCGTTTACGCTTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCaaa                                                                       <  1:1013646/69‑1 (MQ=255)
aaaGATACCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAggg                                                                     >  1:3109039/1‑71 (MQ=255)
aaaGATACCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAg                                                                      >  1:244701/1‑70 (MQ=255)
 aaGATACCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGt                                                                    >  1:2778123/1‑71 (MQ=255)
 aaGATACCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGt                                                                    <  1:2081780/71‑1 (MQ=255)
  aGATACCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAGGGTg                                                                   >  1:1005488/1‑71 (MQ=255)
      aCCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGAt                                                               >  1:2052017/1‑71 (MQ=255)
      aCCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGAt                                                               >  1:4057726/1‑71 (MQ=255)
       ccTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATc                                                              <  1:2995223/71‑1 (MQ=255)
         tctAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATc                                                              <  1:470738/69‑1 (MQ=255)
          ctAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAGGGTGCGATCAc                                                            >  1:4328051/1‑70 (MQ=255)
           tAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCAc                                                            <  1:1546259/69‑1 (MQ=255)
           tAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCa                                                          >  1:1739516/1‑71 (MQ=255)
             aacaaacCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATc                                                        >  1:1086761/1‑71 (MQ=255)
               caaacCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAGGGTGCGATCACCATCGc                                                      <  1:1811893/71‑1 (MQ=255)
                aaacCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCt                                                     <  1:3153064/71‑1 (MQ=255)
                aaacCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCt                                                     <  1:3032568/71‑1 (MQ=255)
                 aacCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCt                                                     <  1:1503845/70‑1 (MQ=255)
                  acCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCAc                                                            >  1:1577792/1‑62 (MQ=255)
                  acCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCAc                                                            >  1:2507260/1‑62 (MQ=255)
                      tttACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGc                                                                  <  1:522985/52‑1 (MQ=255)
                      tttACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGc                                                                  <  1:3387661/52‑1 (MQ=255)
                      tttACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGc                                                      >  1:621040/1‑64 (MQ=255)
                        tACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTa                                                    >  1:4592906/1‑64 (MQ=255)
                        tACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTa                                                    >  1:4598140/1‑64 (MQ=255)
                        tACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTg                                              >  1:2786449/1‑70 (MQ=255)
                        tACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTg                                              >  1:3745838/1‑70 (MQ=255)
                        tACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTg                                              >  1:1156878/1‑70 (MQ=255)
                           gcgTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaaa                                          <  1:2190087/71‑1 (MQ=255)
                              ttCCTCGGTGAAGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaa                                           <  1:140809/67‑1 (MQ=255)
                               tCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTg                                      >  1:2327332/1‑71 (MQ=255)
                                 cTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGa                                    <  1:2279770/71‑1 (MQ=255)
                                 cTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGa                                    <  1:3961642/71‑1 (MQ=255)
                                 cTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGa                                    <  1:4411599/71‑1 (MQ=255)
                                  tCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTg                                      <  1:758976/68‑1 (MQ=255)
                                   cggtgacggtgaAATGGACGAACCGGAATCCAAGGGTGCGATCACCATCGCTACCCGTGaaaaa                                         <  1:3281382/64‑1 (MQ=255)
                                   cggtgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCAcc                                                           >  1:444078/1‑46 (MQ=255)
                                   cggtgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaaaa                                         >  1:3966040/1‑64 (MQ=255)
                                   cggtgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaaaa                                         >  1:3752146/1‑64 (MQ=255)
                                    ggtgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCAc                                                            >  1:574842/1‑44 (MQ=255)
                                    ggtgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTa                                                    >  1:3502890/1‑52 (MQ=255)
                                    ggtgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTa                                                    >  1:3973886/1‑52 (MQ=255)
                                    ggtgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATa                                  <  1:2950962/70‑1 (MQ=255)
                                      tgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAAcc                               <  1:2387152/71‑1 (MQ=255)
                                      tgacggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAAcc                               <  1:1334078/71‑1 (MQ=255)
                                        acggtgaAATGGACGAACCGGAATCCAAAGGTGCGATCa                                                             <  1:713460/39‑1 (MQ=255)
                                         cggtgaAATGGACGAACCGGAATCCAAGGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTgg                            <  1:1940013/71‑1 (MQ=255)
                                          ggtgaAATGGACGAACCGGAATCCAAGGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGt                           <  1:1293938/71‑1 (MQ=255)
                                           gtgaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATc                                                        <  1:2035888/41‑1 (MQ=255)
                                             gaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTa                                                    >  1:284402/1‑43 (MQ=255)
                                             gaAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCtt                        >  1:4048998/1‑71 (MQ=255)
                                              aaaTGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCtt                        <  1:4508824/70‑1 (MQ=255)
                                              aaaTGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCtt                        <  1:3535088/70‑1 (MQ=255)
                                                aTGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAAcc                               <  1:3885796/61‑1 (MQ=255)
                                                aTGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAAc                                >  1:236252/1‑60 (MQ=255)
                                                 tGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaa                                           <  1:1485584/48‑1 (MQ=255)
                                                 tGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaa                                           <  1:1836610/48‑1 (MQ=255)
                                                 tGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaa                                           <  1:4145782/48‑1 (MQ=255)
                                                  ggACGAACCGGAATCCAAGGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTa                   >  1:117240/1‑71 (MQ=255)
                                                  ggACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGtt                    <  1:4249742/70‑1 (MQ=255)
                                                   gACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaaa                                          <  1:62160/47‑1 (MQ=255)
                                                   gACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGaaaa                                          <  1:2270709/47‑1 (MQ=255)
                                                    aCGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATc                 <  1:2529011/71‑1 (MQ=255)
                                                     cGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATa                                  <  1:3217874/53‑1 (MQ=255)
                                                         ccGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATc                 <  1:4107183/66‑1 (MQ=255)
                                                          cGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATCAACTGt           <  1:1582987/71‑1 (MQ=255)
                                                           ggAATCCAAGGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATCAACTGTa          <  1:2452800/71‑1 (MQ=255)
                                                             aaTCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTg                             <  1:3002680/50‑1 (MQ=255)
                                                               tCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATCAACTGTAACCt      <  1:1224313/71‑1 (MQ=255)
                                                                ccAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATa                                  <  1:3710403/42‑1 (MQ=255)
                                                                ccAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTg                             <  1:1556551/47‑1 (MQ=255)
                                                                 cAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGt                     >  1:3544853/1‑54 (MQ=255)
                                                                   aaGGTGCGATCACCATCGCTACCCGTGAAAAACTg                                      <  1:4336009/35‑1 (MQ=255)
                                                                   aaGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATCAACTGTAACCTGCAg  >  1:3436771/1‑71 (MQ=255)
                                                                    |                                                                     
AAAGATACCTCTAAACAAACCGTTTACGCGTTCCTCGGTGACGGTGAAATGGACGAACCGGAATCCAAAGGTGCGATCACCATCGCTACCCGTGAAAAACTGGATAACCTGGTCTTCGTTATCAACTGTAACCTGCAG  >  W3110S.gb/123668‑123805

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: