Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 857,256 A→C 10.9% D38E (GAT→GAG trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE857,2560AC10.9% 139.9 / 17.8 64D38E (GAT→GAGtrpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (24/33);  new base C (3/4);  total (27/37)
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

TGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTA  >  minE/857189‑857324
                                                                   |                                                                    
tGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATc                                                                    <  1:52274/70‑1 (MQ=255)
  gTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTGGTTGAACTGAGCCTGAAATTCAGGATCttt                                                                 <  1:1353749/71‑1 (MQ=255)
  gTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCttt                                                                 <  1:718769/71‑1 (MQ=255)
  gTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCtt                                                                  <  1:3134973/70‑1 (MQ=255)
     aGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTg                                                              >  1:2809623/1‑71 (MQ=255)
     aGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGCATCTTTTTg                                                              <  1:1697111/71‑1 (MQ=255)
       cgcgGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTgcg                                                            >  1:3380647/1‑71 (MQ=255)
       cgcgGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTgcg                                                            >  1:1386844/1‑71 (MQ=255)
        gcgGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTgcg                                                            <  1:3858878/70‑1 (MQ=255)
              ggACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCtttt                                                                >  1:3325947/1‑60 (MQ=255)
               gACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACa                                                    <  1:73897/71‑1 (MQ=255)
                       gCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGcttc                                            >  1:67241/1‑71 (MQ=255)
                       gCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGcttc                                            >  1:3604625/1‑71 (MQ=255)
                         aTAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGcttctt                                          <  1:2860671/71‑1 (MQ=255)
                          tAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTg                                                              <  1:3623178/50‑1 (MQ=255)
                            gTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCttcttc                                         <  1:2916301/69‑1 (MQ=255)
                            gTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCttcttc                                         >  1:3279430/1‑69 (MQ=255)
                            gTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCttcttc                                         >  1:13672/1‑69 (MQ=255)
                            gTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCttcttc                                         <  1:1316591/69‑1 (MQ=255)
                             tttttCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCt                                                                   <  1:1125253/42‑1 (MQ=255)
                             tttttCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCt                                                                   <  1:2901507/42‑1 (MQ=255)
                             tttttCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCt                                                                   <  1:263118/42‑1 (MQ=255)
                             tttttCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCt                                                                   <  1:667592/42‑1 (MQ=255)
                               tttCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCt                                    >  1:2377339/1‑71 (MQ=255)
                                      agGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCag                             <  1:1907357/71‑1 (MQ=255)
                                         tCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGc                          <  1:3680917/71‑1 (MQ=255)
                                            ttGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCag                             <  1:2433706/65‑1 (MQ=255)
                                            ttGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAgg                       >  1:274662/1‑71 (MQ=255)
                                            ttGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAgg                       >  1:1953785/1‑71 (MQ=255)
                                             tGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAgg                       >  1:830004/1‑70 (MQ=255)
                                              gAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGc                      >  1:2893877/1‑70 (MQ=255)
                                              gAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCa                     >  1:3625438/1‑71 (MQ=255)
                                              gAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCa                     >  1:3403320/1‑71 (MQ=255)
                                               aaCTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCAt                    <  1:2909932/71‑1 (MQ=255)
                                               aaCTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCAt                    <  1:1329657/71‑1 (MQ=255)
                                                aCTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCAt                    >  1:3189439/1‑70 (MQ=255)
                                                aCTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCAt                    >  1:390938/1‑70 (MQ=255)
                                                 cTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCa                         >  1:2810026/1‑64 (MQ=255)
                                                    aGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAg                 <  1:327363/69‑1 (MQ=255)
                                                     gCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGAt              <  1:2843103/71‑1 (MQ=255)
                                                      ccTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGAtt             <  1:2851977/71‑1 (MQ=255)
                                                      ccTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGAtt             <  1:614543/71‑1 (MQ=255)
                                                        tGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCAt                    <  1:2431834/62‑1 (MQ=255)
                                                        tGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCAt                    <  1:1920950/62‑1 (MQ=255)
                                                         gAATTTCAGGCTCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAg                 <  1:1252039/64‑1 (MQ=255)
                                                         gAATTTCAGGCTCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAg                 <  1:1819378/64‑1 (MQ=255)
                                                         gAATTTCAGGCTCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAg                 <  1:703896/64‑1 (MQ=255)
                                                         gAATTTCAGGCTCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAg                 <  1:2057816/64‑1 (MQ=255)
                                                         gAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTtg           <  1:1915203/70‑1 (MQ=255)
                                                         gAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCAAGCTGGCGCAGAGCAGGCATCAGGATTtg           <  1:3708969/70‑1 (MQ=255)
                                                          aaaTTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTtgtg         >  1:2678976/1‑71 (MQ=255)
                                                          aaaTTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTtgtg         >  1:3153576/1‑71 (MQ=255)
                                                          aaaTTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTtg           <  1:689077/69‑1 (MQ=255)
                                                           atttCAGGCTCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGa               >  1:3731702/3‑64 (MQ=255)
                                                           atttCAGGCTCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGa               >  1:2061563/3‑64 (MQ=255)
                                                           atttCAGGCTCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGa               >  1:335390/3‑64 (MQ=255)
                                                           aaTTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTtgtg         >  1:2417554/1‑70 (MQ=255)
                                                            aTTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTgg        <  1:3497246/70‑1 (MQ=255)
                                                            aTTCAGGATCTTTTTGCGCACTGACAAAAGCTTATTCCAGCTGGCGCaga                            >  1:2929733/1‑50 (MQ=255)
                                                             ttCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCa      >  1:3594557/1‑71 (MQ=255)
                                                               cAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTtgt          >  1:3542889/1‑65 (MQ=255)
                                                                aGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGt   <  1:2182044/71‑1 (MQ=255)
                                                                aGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGt   <  1:3223019/71‑1 (MQ=255)
                                                                 ggATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTa  <  1:2510068/71‑1 (MQ=255)
                                                                   |                                                                    
TGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTA  >  minE/857189‑857324

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: