Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,647,451 T→C 16.2% intergenic (+131/‑142) asd → / → glgB aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding/1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,647,4510TC16.2% 171.1 / 19.2 104intergenic (+131/‑142)asd/glgBaspartate‑semialdehyde dehydrogenase, NAD(P)‑binding/1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (41/46);  new base C (7/10);  total (48/56)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 123 total aligned reads displayed.

CAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCT  >  minE/2647384‑2647508
                                                                   |                                                         
cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                          <  1:1997096/69‑1 (MQ=255)
cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                          <  1:578465/69‑1 (MQ=255)
cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                          <  1:1658248/69‑1 (MQ=255)
cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                          <  1:2430463/69‑1 (MQ=255)
cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                          <  1:2773068/69‑1 (MQ=255)
cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                          <  1:137789/69‑1 (MQ=255)
 aGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                          <  1:2287680/68‑1 (MQ=255)
    tGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTa                                                      <  1:2388023/69‑1 (MQ=255)
    tGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTa                                                      <  1:1919283/69‑1 (MQ=255)
        tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                    >  1:494705/1‑67 (MQ=255)
        tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                    >  1:741135/1‑67 (MQ=255)
        tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                    >  1:546999/1‑67 (MQ=255)
        tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAt                                                     >  1:2980491/1‑66 (MQ=255)
            gATTTACCGTGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                              <  1:63403/69‑1 (MQ=255)
            gATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                              <  1:1693502/69‑1 (MQ=255)
              tttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCt                                            >  1:1226141/1‑69 (MQ=255)
              tttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCt                                            >  1:2607754/1‑69 (MQ=255)
              tttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCt                                            >  1:1511177/1‑69 (MQ=255)
               ttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                           <  1:1068668/69‑1 (MQ=255)
                tACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                           <  1:477158/68‑1 (MQ=255)
                    gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTg                                      <  1:2502947/69‑1 (MQ=255)
                    gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTg                                      <  1:527612/69‑1 (MQ=255)
                    gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTg                                      <  1:118461/69‑1 (MQ=255)
                    gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTg                                      <  1:866998/69‑1 (MQ=255)
                              tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGa                                     >  1:2151902/1‑60 (MQ=255)
                              tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                            >  1:2035908/1‑69 (MQ=255)
                              tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                            >  1:1599283/1‑69 (MQ=255)
                              tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                            >  1:29352/1‑69 (MQ=255)
                              tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                            >  1:3122195/1‑69 (MQ=255)
                                      tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                      >  1:542727/1‑67 (MQ=255)
                                      tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                      >  1:2243395/1‑67 (MQ=255)
                                      tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                      >  1:1158496/1‑67 (MQ=255)
                                      tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                      >  1:3321864/1‑67 (MQ=255)
                                      tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                      >  1:3283615/1‑67 (MQ=255)
                                      tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                      >  1:1525475/1‑67 (MQ=255)
                                       ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagag                   >  1:2121939/1‑69 (MQ=255)
                                         cTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                              <  1:1386743/40‑1 (MQ=255)
                                            ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa               >  1:1518100/1‑68 (MQ=255)
                                            ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa               >  1:1244385/1‑68 (MQ=255)
                                            ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa               >  1:2720662/1‑68 (MQ=255)
                                            ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa               >  1:2485125/1‑68 (MQ=255)
                                            ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa               >  1:3233670/1‑68 (MQ=255)
                                             tCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                              <  1:401763/36‑1 (MQ=255)
                                              ctgtAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:1465274/4‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                           >  1:2318895/1‑38 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                           >  1:305415/1‑38 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                           >  1:223759/1‑38 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                           >  1:1476779/1‑38 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                           >  1:2184567/1‑38 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:461030/64‑1 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:685394/64‑1 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:1872992/1‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:2175501/64‑1 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:2978180/1‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:1391759/1‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:1218048/64‑1 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:524567/1‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:1906661/64‑1 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:1551166/1‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:1185947/64‑1 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:1681289/1‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 >  1:2259811/1‑64 (MQ=255)
                                              cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATg             <  1:286689/68‑1 (MQ=255)
                                                   aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca         <  1:408580/67‑1 (MQ=255)
                                                   aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca         <  1:1420721/67‑1 (MQ=255)
                                                   aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca         <  1:2870167/67‑1 (MQ=255)
                                                   aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca         <  1:161673/67‑1 (MQ=255)
                                                   aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca         <  1:624999/67‑1 (MQ=255)
                                                     ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag      <  1:8437/68‑1 (MQ=255)
                                                     ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag      <  1:817173/68‑1 (MQ=255)
                                                     ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag      <  1:981685/68‑1 (MQ=255)
                                                     ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag      <  1:1189603/68‑1 (MQ=255)
                                                     ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag      <  1:1440131/68‑1 (MQ=255)
                                                     ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag      <  1:2744820/68‑1 (MQ=255)
                                                     ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag      <  1:2750321/68‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:2553876/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:1400197/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:697333/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:1726701/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:1861793/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:3257593/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:2718412/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:3091221/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:3097602/56‑1 (MQ=255)
                                                      gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                 <  1:321531/56‑1 (MQ=255)
                                                       tGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgc   >  1:1320555/1‑69 (MQ=255)
                                                       tGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgc   >  1:1747538/1‑69 (MQ=255)
                                                            gAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca         <  1:731155/58‑1 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagaga                  >  1:641890/1‑45 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:2205955/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:562816/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:487769/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:401523/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:372336/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:2463595/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:362734/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:3330791/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:1391955/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:3209909/1‑61 (MQ=255)
                                                                acacGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct  >  1:3223247/1‑61 (MQ=255)
                                                                   |                                                         
CAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCT  >  minE/2647384‑2647508

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: