Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 442,284 G→A 48.0% P493L (CCA→CTA)  speF ← ornithine decarboxylase isozyme, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE442,2840GA48.0% 11.2 / 144.4 100P493L (CCA→CTA) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base A (21/27);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 107 total aligned reads displayed.

TACACCGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTGGTACAAACTGG  >  minE/442217‑442352
                                                                   |                                                                    
tACACCGAACGCTTCATATTCGCCGTTACGTGCATCAATACATGGGGTTGTCAGCAACAGTTTGCATGGAt                                                                   >  1:1023830/1‑71 (MQ=255)
  cacCGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATc                                                                  >  1:1735089/1‑70 (MQ=255)
  cacCGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATc                                                                  >  1:3458637/1‑70 (MQ=255)
   acCGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATgg                                                                     <  1:1782622/66‑1 (MQ=255)
    ccGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATgg                                                                     <  1:3247504/65‑1 (MQ=255)
    ccGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATgg                                                                     <  1:1051623/65‑1 (MQ=255)
    ccGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATgg                                                                     <  1:589859/65‑1 (MQ=255)
     cGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACa                                                              <  1:1741274/71‑1 (MQ=255)
      gAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGa                                                                    >  1:797108/1‑64 (MQ=255)
      gAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGa                                                                    >  1:3868195/1‑64 (MQ=255)
      gAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACaa                                                             <  1:3001491/71‑1 (MQ=255)
       aaCGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACaa                                                             <  1:2151961/70‑1 (MQ=255)
        aCGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAg                                                           <  1:3821715/71‑1 (MQ=255)
        aCGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAg                                                           <  1:2606877/71‑1 (MQ=255)
              cATATTCGCCGTTACGTGCATCAATACCTGGGGTTTTCAGCAACAGTTTGCATAGATCGACAAAGTa                                                         <  1:2535895/67‑1 (MQ=255)
              cATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTAtt                                                       <  1:457797/69‑1 (MQ=255)
               atatTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGa                                                     <  1:3139390/70‑1 (MQ=255)
               atatTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGa                                                     <  1:4139131/70‑1 (MQ=255)
               atatTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGa                                                     <  1:2308571/70‑1 (MQ=255)
                tatTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGAt                                                    <  1:1786794/70‑1 (MQ=255)
                tatTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGa                                                                    >  1:1603546/1‑54 (MQ=255)
                  ttCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATtc                                                  >  1:1959646/1‑70 (MQ=255)
                  ttCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTa                                                         >  1:1166360/1‑63 (MQ=255)
                  ttCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTAtt                                                       <  1:1045554/65‑1 (MQ=255)
                  ttCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATtct                                                 >  1:2123020/1‑71 (MQ=255)
                  ttCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATtct                                                 >  1:2113235/1‑71 (MQ=255)
                   tCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATtctc                                                >  1:3355395/1‑71 (MQ=255)
                   tCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATtctc                                                >  1:417822/1‑71 (MQ=255)
                    cGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTctct                                               >  1:3238610/1‑71 (MQ=255)
                       cGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCa                                            <  1:1804806/71‑1 (MQ=255)
                         ttACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTg                                              <  1:3001273/67‑1 (MQ=255)
                           aCGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGa                                                     >  1:2736587/1‑58 (MQ=255)
                           aCGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGc                                        >  1:4120998/1‑71 (MQ=255)
                           aCGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACaaa                                                            >  1:734040/1‑51 (MQ=255)
                           aCGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTAtt                                                       >  1:198435/1‑56 (MQ=255)
                           aCGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTAtt                                                       >  1:709346/1‑56 (MQ=255)
                           aCGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTAtt                                                       >  1:504198/1‑56 (MQ=255)
                           aCGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAg                                         <  1:793479/70‑1 (MQ=255)
                            cGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGcc                                       <  1:3440367/71‑1 (MQ=255)
                            cGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGcc                                       >  1:3183842/1‑71 (MQ=255)
                            cGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGc                                        <  1:1530662/70‑1 (MQ=255)
                             gTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCt                                      <  1:1550399/71‑1 (MQ=255)
                              tGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCa                                            <  1:1806326/64‑1 (MQ=255)
                                cATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATa                                          >  1:1425455/1‑64 (MQ=255)
                                cATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGAt                                                                   >  1:3309632/1‑39 (MQ=255)
                                cATCAATACCTGGGGTTGTCAGCAACAGTTTGCATAGAt                                                                   >  1:712203/1‑39 (MQ=255)
                                  tCAATACCTGGGGTTGTCAGCAACAGTTTGCATgg                                                                     <  1:738009/35‑1 (MQ=255)
                                  tCAATACCTGGGGTTGTCAGCAACAGTTTGCATgg                                                                     <  1:1153156/35‑1 (MQ=255)
                                  tCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCaaa                                 <  1:346938/71‑1 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACaaa                                                            >  1:4006197/1‑42 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATa                                          >  1:1248668/1‑60 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCa                                   >  1:1910876/1‑67 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTa                                                         >  1:1175041/1‑45 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATtct                                                 <  1:4016282/53‑1 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCaa                                  >  1:711542/1‑68 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAg                               <  1:1353889/71‑1 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAg                               <  1:2674750/71‑1 (MQ=255)
                                    aaTACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAg                               <  1:804958/71‑1 (MQ=255)
                                     aTACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAg                               >  1:1212003/1‑70 (MQ=255)
                                      tACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTg                                                      <  1:3164618/46‑1 (MQ=255)
                                      tACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATa                                          <  1:512677/58‑1 (MQ=255)
                                      tACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATa                                          <  1:965101/58‑1 (MQ=255)
                                      tACCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGc                                        >  1:2472745/1‑60 (MQ=255)
                                       aCCTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAAt                            >  1:750987/1‑71 (MQ=255)
                                        ccTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAAt                            <  1:2370127/70‑1 (MQ=255)
                                        ccTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCa                                            >  1:2614872/1‑54 (MQ=255)
                                         cTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGt                                                          <  1:3252579/39‑1 (MQ=255)
                                         cTGGGGTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGc                          <  1:1001659/71‑1 (MQ=255)
                                           ggggTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGc                          >  1:1653349/1‑69 (MQ=255)
                                            gggTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGcc                         <  1:576322/69‑1 (MQ=255)
                                            gggTTGTCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGcc                         <  1:1055412/69‑1 (MQ=255)
                                               ttGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCa                                            <  1:1327066/47‑1 (MQ=255)
                                               ttGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCaa                                  <  1:420899/57‑1 (MQ=255)
                                               ttGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATg                           <  1:3664833/64‑1 (MQ=255)
                                               ttGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATg                           <  1:2612303/64‑1 (MQ=255)
                                                 gTCAGCAACAGTTTGCATGGATCGACAAAGTATTGAt                                                    >  1:2385461/1‑37 (MQ=255)
                                                 gTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATa                                          >  1:3938498/1‑47 (MQ=255)
                                                 gTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGaa                             >  1:1768808/1‑60 (MQ=255)
                                                  tCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATg                           <  1:1502988/61‑1 (MQ=255)
                                                  tCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAg                               <  1:2969673/57‑1 (MQ=255)
                                                  tCAGCAACAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTcc                 <  1:3590610/71‑1 (MQ=255)
                                                    aGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGaa                             >  1:1483513/1‑57 (MQ=255)
                                                    aGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTcccc               >  1:3357160/1‑71 (MQ=255)
                                                    aGCAACAGTTTGCATGGATCGACAAAGTATTGATTATCTGCATAGCCTTCAAAAGaa                             >  1:56127/1‑57 (MQ=255)
                                                      cAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTg             >  1:3169515/1‑71 (MQ=255)
                                                        aCAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCaa                                  <  1:2677714/48‑1 (MQ=255)
                                                        aCAGTTTGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGc                          <  1:1749888/56‑1 (MQ=255)
                                                            tttGCATGGATCGACAAAGTATTGATTCTCTGCATAGc                                        <  1:2293483/38‑1 (MQ=255)
                                                            tttGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTGGTACa        >  1:822310/1‑70 (MQ=255)
                                                             ttGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAg                       >  1:3803957/1‑54 (MQ=255)
                                                              tGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTc                  >  1:2435833/1‑58 (MQ=255)
                                                              tGCATAGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTGGTACAAAc     >  1:3299530/1‑71 (MQ=255)
                                                               gCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTGGTACAAAc     <  1:2685407/70‑1 (MQ=255)
                                                               gCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTGGTACAAAc     <  1:3596375/70‑1 (MQ=255)
                                                                cATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGaa                             >  1:23615/1‑45 (MQ=255)
                                                                 aTGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAgt                      <  1:1930063/51‑1 (MQ=255)
                                                                 aTGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTc                  <  1:1519808/55‑1 (MQ=255)
                                                                 aTGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTc                  <  1:1518202/55‑1 (MQ=255)
                                                                 aTGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTGGTACAAACTgg  <  1:3959533/71‑1 (MQ=255)
                                                                   ggATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTccc                >  1:3712129/1‑55 (MQ=255)
                                                                   |                                                                    
TACACCGAACGCTTCATATTCGCCGTTACGTGCATCAATACCTGGGGTTGTCAGCAACAGTTTGCATGGATCGACAAAGTATTGATTCTCTGCATAGCCTTCAAAAGAATGCCAGTGTTCCCCTGGTACAAACTGG  >  minE/442217‑442352

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: