Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,358 A→C 12.8% I22L (ATC→CTC) ‡ setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,3580AC12.8% 184.9 / 24.8 86I22L (ATC→CTC) ‡setBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base A (45/30);  new base C (2/9);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 1.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGA  >  minE/1395290‑1395428
                                                                    |                                                                      
ccAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATc                                                                      <  1:894684/71‑1 (MQ=255)
 cAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCg                                                                     <  1:2851816/71‑1 (MQ=255)
 cAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCg                                                                     <  1:3157539/71‑1 (MQ=255)
     atgcatAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGc                                                                 <  1:5144516/71‑1 (MQ=255)
        catAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCtt                                                              >  1:2210088/1‑71 (MQ=255)
        catAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCt                                                               <  1:5021084/70‑1 (MQ=255)
          tAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATc                                                                      <  1:2138289/61‑1 (MQ=255)
          tAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTtc                                                            >  1:5070636/1‑71 (MQ=255)
          tAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTtc                                                            >  1:685511/1‑71 (MQ=255)
           aaCTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTtc                                                            >  1:4036094/1‑70 (MQ=255)
            aCTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTtct                                                           >  1:2071719/1‑70 (MQ=255)
                ccccGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTtctc                                                          <  1:228408/67‑1 (MQ=255)
                 cccGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGtt                                                                   >  1:3771341/1‑57 (MQ=255)
                 cccGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                     <  1:2579060/71‑1 (MQ=255)
                  ccGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGcc                                                                <  1:5264894/59‑1 (MQ=255)
                  ccGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAc                                                        >  1:2706281/1‑67 (MQ=255)
                  ccGCAGTCTCCAGCGCGAAATCGTATGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                     >  1:635835/1‑70 (MQ=255)
                   cGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                   >  1:5225024/1‑71 (MQ=255)
                    gCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                  >  1:5306107/1‑71 (MQ=255)
                    gCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                  >  1:660005/1‑71 (MQ=255)
                    gCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGGTTTTAATCGTTGCCTTTCTCACCGGTAt                                                  >  1:1389234/1‑71 (MQ=255)
                       gTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATc                                                                      <  1:931724/48‑1 (MQ=255)
                        tctcCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCg                                                                     <  1:1807780/48‑1 (MQ=255)
                         ctcCAGCGCGAAATCGTTTGACCGGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCgg                                             >  1:3612387/1‑71 (MQ=255)
                          tcCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCg                                                      <  1:968773/61‑1 (MQ=255)
                          tcCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCggg                                            >  1:4405991/1‑71 (MQ=255)
                            cAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGt                                                                    >  1:2444831/1‑45 (MQ=255)
                            cAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                          >  1:4176297/1‑71 (MQ=255)
                            cAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                          >  1:729124/1‑71 (MQ=255)
                              gcgcGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGtt                                                                   <  1:2222784/44‑1 (MQ=255)
                              gcgcGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGtt                                                                   <  1:2995409/44‑1 (MQ=255)
                              gcgcGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCa                                                         >  1:4068343/1‑54 (MQ=255)
                              gcgcGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCa                                                         >  1:2890089/1‑54 (MQ=255)
                              gcgcGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCa                                                         >  1:1899379/1‑54 (MQ=255)
                              gcgcGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCa                                                         >  1:861092/1‑54 (MQ=255)
                                 cgAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                     >  1:1719305/1‑71 (MQ=255)
                                 cgAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                     >  1:1707370/1‑71 (MQ=255)
                                 cgAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                     >  1:7046/1‑71 (MQ=255)
                                  gAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGt                                                    <  1:4810592/55‑1 (MQ=255)
                                  gAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGt                                                    <  1:2350916/55‑1 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCttt                                                             <  1:1821564/45‑1 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTtct                                                           <  1:616655/47‑1 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCa                                                         >  1:1568340/1‑49 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCa                                                         >  1:4271485/1‑49 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAcc                                                       <  1:1375379/51‑1 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAcc                                                       <  1:1220346/51‑1 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAcc                                                       <  1:2395392/51‑1 (MQ=255)
                                   aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGt                                                    <  1:3022663/54‑1 (MQ=255)
                                    aaTCGTTTGACCTGACCTCGAGGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                  >  1:5423496/1‑55 (MQ=255)
                                    aaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                  >  1:1540246/1‑55 (MQ=255)
                                     aTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                     <  1:1561261/51‑1 (MQ=255)
                                     aTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                     <  1:902423/51‑1 (MQ=255)
                                     aTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                     <  1:2479554/51‑1 (MQ=255)
                                       cGTTTGCCCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAc                                                        >  1:4283757/1‑46 (MQ=255)
                                       cGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAc                                                        >  1:1054151/1‑46 (MQ=255)
                                        gTTTGACCTGACCTCGACGGCGTTTTTAATCGTTg                                                                  >  1:210461/1‑35 (MQ=255)
                                        gTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                              >  1:2179918/1‑71 (MQ=255)
                                         tttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTtctc                                                          <  1:2904769/42‑1 (MQ=255)
                                          ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                   >  1:4508323/1‑48 (MQ=255)
                                          ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                   >  1:459019/1‑48 (MQ=255)
                                          ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                   >  1:3335461/1‑48 (MQ=255)
                                          ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAcccc                             >  1:1240193/1‑70 (MQ=255)
                                                 gacctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                          <  1:3104764/50‑1 (MQ=255)
                                                     tCGGCGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAcccc                             <  1:5507853/59‑1 (MQ=255)
                                                     tCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCaa                                  <  1:3094563/54‑1 (MQ=255)
                                                      ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCTCTCTTCAGACTCCTACCCTAAGTa                  <  1:1912440/66‑1 (MQ=255)
                                                      ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                  <  1:5461045/66‑1 (MQ=255)
                                                      ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                  <  1:5606809/66‑1 (MQ=255)
                                                      ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                  <  1:961379/66‑1 (MQ=255)
                                                      cGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                          >  1:4075428/1‑45 (MQ=255)
                                                      cGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                              >  1:570409/1‑57 (MQ=255)
                                                       gACGGCGTTTATAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTAttt               >  1:345318/1‑71 (MQ=255)
                                                         cGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTAtattc             <  1:1329297/71‑2 (MQ=255)
                                                           gcatttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCt           <  1:5433845/68‑1 (MQ=255)
                                                           gcatttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCt           <  1:433957/68‑1 (MQ=255)
                                                             gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa         >  1:3885576/1‑71 (MQ=255)
                                                             atttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCtcg         >  1:3109744/2‑69 (MQ=255)
                                                              tttttACTTGTCGCCTTTCTGACGGGTATTGCGGgcg                                          <  1:1611621/37‑1 (MQ=12)
                                                              tttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCa                                   <  1:4844851/44‑1 (MQ=16)
                                                              tttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                           <  1:4913978/52‑1 (MQ=255)
                                                               ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc       >  1:2495962/1‑71 (MQ=255)
                                                               ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc       >  1:5049856/1‑71 (MQ=255)
                                                                tttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCtaa                     >  1:3363584/1‑54 (MQ=255)
                                                                tttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                     >  1:4477010/1‑40 (MQ=255)
                                                                tttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc       >  1:4926033/1‑70 (MQ=255)
                                                                 ttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGa     >  1:4581485/1‑71 (MQ=255)
                                                                    aTCGTTGCCTTTCTCACCGGTATTGCGGGTGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATg   <  1:2188262/70‑1 (MQ=255)
                                                                    aTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGa  <  1:4649244/71‑1 (MQ=255)
                                                                    |                                                                      
CCAGCATGCATAACTCCCCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGA  >  minE/1395290‑1395428

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: